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PREFACE 


On  26-27  June  1984  the  Second  Tri-Service  Cloud  Modeling 
Workshop  was  held  at  the  Naval  Surface  Weapons  Center,  White 
Oak ,  MD . 

The  purpose  of  the  meeting  was  the  exchange  of  information 
between  users  of  cloud  data  information,  representing  the  devel¬ 
opers  of  electrooptical  surveillance  (and  other)  systems,  and 
the  producers  of  cloud  information  in  the  DoD  weather  community. 

This  document  summarizes  the  results  of  the  workshop  and 
related  work.  Part  I  consists  of  the  DoD  letter  establishing 
the  Tri-Service  Cloud  Modeling  Effort,  the  letter  of  invitation 
and  program  of  the  meeting,  and  a  brief  summary  report.  Part  II 
consists  of  the  briefing  charts  and  supplementary  material  - 
classified  material  is  included  as  a  supplementary  volume.  Part 
III  consists  of  the  list  of  invitees  and  Part  IV  of  the  list  of 
attendees,  while  Part  V  consists  of  supplementary  comments 
received  after  the  meeting,  which  serve  to  expand  on  the  summary 
report.  Finally,  Part  VI  summarizes  a  related  meeting  held  at 
IDA  on  24  January  1984  on  the  use  of  airborne  lidars  to  collect 
meteorological  data. 

The  preparation  of  this  document  was  supported  by  DARPA/DEO 
under  IDA  Task  A-88.  The  material  contained  herein  is  presented 
for  the  use  and  convenience  of  interested  parties.  It  has  not 
been  evaluated,  analyzed,  or  subjected  to  IDA  review. 


SUMMARY 


The  Second  Tri-Service  Cloud  Modeling  Workshop  was  held  at 
the  Naval  Surface  Weapons  Center,  White  Oak,  Maryland,  on  26-27 
June  1984. 

Clouds  affect  DoD  systems — high-energy  laser  weapons 
systems,  optical  surveillance  systems  and  some  others — in  a 
wide  variety  of  ways.  These  cloud  effects  are  at  times  not 
well  understood  and  even  more  frequently  not  well  described. 
Indeed,  there  is  a  major  gap  between  the  existence  of  cloud 
data  and  its  utilization  for  the  realistic  evaluation  of  DoD 
weapons  systems  performance.  To  quote  from  the  DoD  letter 
establishing  the  Tri-Service  CMW  (see  p.  3  for  the  full  letter): 

"The  realistic  evaluation  of  DoD  weapon  system  performance 
is  a  critical  element  of  the  multimillion/billion  dollar  acqui¬ 
sition  decisions.  The  weather  limitations  on  the  performance 
must  play  a  key  role  in  the  evaluation  process. 

"We  (DoD  meteorologists)  have  worked  very  hard  to  obtain 
recognition  of  the  weather  effects  by  the  operational  evaluation 
community;  however,  now  that  we  have  achieved  this  recognition, 
we  are  failing  to  produce  even  the  most  basic  binary  cloud  data 
and  methodologies  needed  to  support  the  evaluation  programs." 

In  view  of  this  recognized  deficiency,  the  purpose  of  the 
Tri-Service  Cloud  Modeling  Workshop  is  the  exchange  of  informa¬ 
tion  between  users  of  information  on  clouds,  representing  the 
developers  of  a  variety  of  DoD  systems,  and  the  producers  of 
cloud  information  in  the  weather  community. 

The  work  reported  herein  falls  into  several  broad  categor¬ 
ies,  namely: 

•  Data  bases  and  their  limitations — see  pp.  67,  249,  517, 
537,  569,  599,  627 
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•  Cloud  models — see  pp .  17,  45,  119,  323,  377 

•  Cloud  description  (small  and  large  scale) — see  pp. 

145,  201,  311,  353,  417,  471,  501,  613,  C-1* 

•  System  needs--see  pp,  83,  135,  685,  C-22* 

•  System  applications — see  pp.  99,  199 
The  following  general  comments  can  be  made: 

1.  The  meeting  provided  a  forum  for  systems  users 
(however  defined)  to  learn  what  kinds  of  cloud 
data  exist  and  how  related  problems  are  being 
addressed,  for  systems  analysts  to  report  on 
how  they  use  cloud  data,  and  for  the  people 
who  generate  the  data  to  see  what  kind  of 
questions  are  being  asked. 

2.  There  was  an  exchange  of  information  with  non- 
DoD  agencies,  with  presentations  made  by  NASA 
and  NOAA  (see  pp .  537,  569). 

A  brief  report  by  the  steering  committee  which  discusses 
some  outstanding  technical  problems  is  given  on  p.  11.  The 
generic  problem  of  the  lack  of  adequate  interaction  between 
cloud  users  and  cloud  developers  which  led  to  the  establishment 
of  the  Tri-Service  Cloud  Modeling  Workshop  is  still  in  evidence; 
it  can  be  mitigated  by  actions  on  the  part  of  the  steering 
committee,  which  must  actively  seek  a  balance  in  terms  of 
stronger  user  representation,  and  also  serve  as  a  focus  on 
cloud  data  bases  and  methodologies  for  users.  To  achieve  these 
actions  requires  explicit  support. 


^See  classified  supplementary  volume. 


CONTENTS 


PREFACE 

iii 

SUMMARY 

V 

PART  I  . 

Dr.  Edith  W.  Martin's  Letter 
Letter  of  Invitation  to  the  Workshop 
Administrative  Letter 
Agenda 


Steering  Committee's  Report  on  the  Workshop  1 

PART  II.  BRIEFING  CHARTS  AND  SUPPLEMENTARY  MATERIALS  15 

CLOUDS  PROGRAM — A  BRIEF  PERSONAL 

HISTORY/PURPOSE  -  Lt  Col  Vernon  Bliss  17 

OVERVIEW  OF  METHODOLOGIES  AND  CAVEATS  FOR 
MODELING  CFLOS,  CFFOV  AND  CFI  - 

Donald  D.  Grantham  45 

CLOUD  DATA  BASES:  OVERVIEW  AND  CAVEATS  - 

J .  Bunting  67 

DARPA'S  INTEREST  IN  CLOUD  DATA  - 

David  Zimmerman  83 

EVALUATION  OF  THE  EFFECTS  OF  CLOUDS  ON  SBL 

APPLICATIONS  -  Ronald  J.  Nelson  99 

CLOUD  MODELING  EFFORTS  AT  THE  NAVAL 

ENVIRONMENTAL  PREDICTION  RESEARCH  FACILITY  - 
Paul  M.  Tag  119 

TRIDENT  TARGETING  AND  TEST  FLIGHT  ANALYSES  - 

Susan  Masters  135 

SYNTHETIC  CLOUD  SCENES  FOR  IRST  SENSOR  DESIGN 
APPLICATIONS  -  Alex  T.  Maksymowicz, 

J oh nH. Allen  145 


INFRARED  CLOUD  BACKGROUNDS  AND  SENSOR 

PERFORMANCE  -  W.J.  Tropf,  B.J.  Sandford, 
J.H.  Schummers,  J.  schroeder 

vi  i 


199 


ro  lo  ^ 


THE  NAVY  BACKGROUND  MEASUREMENT  AND  ANALYSIS 

PROGRAM  (BMAP)--FY  1984  -  Bernard  V.  Kessler  249 

AN  EMPIRICAL  CLOUD  MODEL  FOR  THE  SWIR  AND  LWIR 

PASS  BANDS  -  Mike  Scarborough,  Bob  Pilgrim  309 

ASSESSMENT  OF  THE  LOWTRAN  6  CIRRUS  MODEL  - 

John  Hornstein  323 

PROFILES  OF  OPTICAL  EXTINCTION  COEFFICIENTS 
CALCULATED  FROM  DROPLET  SPECTRA  OBSERVED  IN 
MARINE  STRATUS  CLOUD  LAYERS  - 


V.  Ray  Noonkester  351 

INTERACTIVE  SOFTWARE  FOR  CALCULATING  CLOUD- 

FREE  INTERVALS  -  Lt  Col  Vernon  Bliss  377 

3-D  CLOUD  SIMULATION  USING  THE  SAWTOOTH 

MODEL  -  Irving  I.  Gringorten  397 

CORRELATION  FUNCTIONS  FOR  CLOUD  COVER  AND 

CEILING  -  Ralph  Shapiro  417 

CLOUD  VISIBILITY  MODELING  OF  JOINT  MESOSCALE 
PROBABILITY  IN  CENTRAL  EUROPE  - 

Oscar  M.  Essenwanger  471 

A  DISCUSSION  OF  THE  SYNTHETIC  IR  CLOUD  SCENES 
OF  PHOTON  RESEARCH  ASSOCIATES  - 

Edward  J.  Stone  501 

REAL-TIME  NEPHANALYSIS  (RTNEPH)  - 

Lt  Col  William  M.  Cox  517 

CLOUD  DATA  BASES  FROM  THE  INTERNATIONAL 

SATELLITE  CLOUD  CLIMATOLOGY  PROJECT  (ISCCP)  - 
Herbert  Jacobowitz  537 

MULTI-YEAR  GLOBAL  CLOUD  DATA  SET  FROM  NIMBUS-7 

Satellite  -  P.  Hwang,  L.  Stowe,  P.K.  Bhartia  569 

VIEWING  ANG1.E  BIAS  IN  CLOUD  CLIMATOLOGIES  - 

J.  Bunting,  G.  Gustafson,  J.  Arck  599 


DIGITAL  IMAGING  AND  ANALYSIS  USING  CLOUD 
CLIMATOLOGY  DATA  (FROM  SATELLITES)  - 
T.  Vonder  Haar,  T.  Brubaker 


613 


A  NEW  DATA  BASE  OF  CLOUD  VARIABLES  FOR 
ALTITUDES  UP  TO  10,000  FEET  AGL  AND  THE 
IMPLICATIONS  FOR  LOW  ALTITUDE  AIRCRAFT  ICING 
Richard  K  .  Jeck 

PART  III.  LIST  OF  INVITEES 

PART  IV.  LIST  OF  ATTENDEES 

PART  V.  SUPPLEMENTARY  COMMENTS 

PART  IV.  ON  THE  USE  OF  AIRBORNE  LIDARS  TO  COLLECT 

METEOROLOGICAL  DATA  (WINDS,  DENSITIES, 
AEROSOLS)  DURING  MISSILE  REENTRY  OVER 
BROAD  OCEAN  AREAS 


OFFICE  OF  THE  UNDER  SECRETARY  OF  DEFENSE 
WASHiNCt'bN  DC  20)01 


10  August  1982 


MOORANDUM  FOR  DIRECTOR  OF  LABORATORIES,  HEADQUARTERS  AIR  POROE 
SYSTEMS  OOMRND 

SUBJECT:  Cloud  Models  for  Realistic  Eveduation  of  DoO  Weapon  System 
Perfotmanoe 


■Riis  is  to  request  seme  assistance  in  providing  planning  guidance  for  a  tri- 
Service  initiative  to  improve  meteorological  support  for  the  DoD  performance 
evrduation  activities. 

The  realistic  evaluation  of  DoD  weapon  system  perfonwice  is  a  critical  element 
of  the  maltimillion/billion  dollar  acquisition  decisions  which  must  often  be 
made  by  DoD.  The  weather  limitations  on  the  performance  of  weapon/sensor  systems 
must  play  a  key  role  in  the  evaluation  process.  Utility/engagement  models  and 
war  gaming  models  all  need  realistic  weather  models  to  handle  the  historical/cllma- 
toLogical  data  needed  to  insure  realistic  and  relevant  results.  Therefore, 
the  meteorological  oonriunity  must  be  prepared  to  support  the  operational  evedu- 
ation,  acquisition,  studies  and  analysis,  and  war  gamng  oemnunities  with  the 
historical  data  bases,  methodologies  and  ccmpwter  programs  which  can  efficiently 
interact  with  the  system  perfonnance  models  to  produce  the  needed  realistic 
output. 

We  have  been  working  to  adequately  address  the  clear  air  atmospheric  transmission 
problems  with  sane  success,  but  still  have  major  difficulties  handling  cloud 
data  effectively  and  efficiently,  especially  in  support  of  the  evaluation  of 
systems  needing  a  cloud-free  line-of-sight  (CFLOS)  or  cloud-free  field-of-view 
(CFFCV)  for  successful  operations.  We  have  vorked  very  hard  to  obtain  recog¬ 
nition  of  the  weather  effects  by  the  operational  evaluation  camunity;  however, 
rcM  that  we  have  achieved  this  recognition,  we  are  failing  to  produce  even  the 
most  basic  binary  cloud  data  and  methodologies  needed  to  support  the  eveduation 
programs.  We  need  to  pull  together  now  and  resolve  seme  of  the  recurring  issues 
in  cloud  modeling  (temporal  and  spatiaQ.  distribution  not  microphysical  modeling) . 
This  includes  identifying  the  complementary  efforts  by  various  organizations 
needed  to  achieve  the  goal  of  a  DoD  accepted  cloud  modeling  methodology  and 
data  beise  to  support  DoD  performance  evaluation  and  war  gaming  activities.  In 
addressing  systems  like  Trident,  Pershing,  or  MX;  Copperhead,  GBU-15  or  Harpoon; 
space  based  laser,  laser  communications  or  HEL  and  all  the  follow-on  systems, 
the  meteorological  oonnunity  roust  be  prep^ared  to  support  the  performanoe  evalu¬ 
ation  models  supporting  the  acquisition  of  these  systens. 

In  view  of  the  above,  I  would  like  the  Air  Force  Geophysics  Laboratory  to  take 
the  lead  and  work  with  the  Navad  &ivironmental  Prediction  Rese2irch  Facility 
and  the  Engineer  Topographical  Laboratories  to  form  a  tri-Service  plcinning 
group  to  1)  identify  the  key  issues  to  be  resolved,  2)  outline  the  work  needed, 

3)  discuss  methods  for  achieving  the  goal,  4)  identify  agencies  best  able  to 
provide  solutions,  and  5)  draft  a  charter  which  would  cover  these  key  issues 
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and  provide  for  continuing  group  activities  until  the  tasks  are  ocmpleted. 

While  related  to  the  current  DARPA  sponsored  CFTjOS  project,  the  DoO-wide  effort 
I  cun  initiating  is  intended  to  be  broader  and  of  longer  duration  to  insure 
advances  are  obtained  to  support  many  types  of  DoO  systems  affected  by  clouds. 

I  would  very  much  appreciate  your  support  and  assistance  in  achieving  this 
goal  of  a  DoD  accepted  cloud  modeling  methodology  and  data  base  to  support  the 
performance  evaluation  and  war  gaming  activities  of  [)oD.  Ihe  initial  planning 
meeting  should  be  held  within  60  days. 

_ 

Edith  W.  Martin 

j  Deputy  Under  Secretary  of  Defense 
>r**^'^-fbr  Research  and  Engineering 
/j  (Research  and  Advanced  Technlogy) 

oc;  Dr.  McClatchey  ^ 

LtCol  Pfeffer 
Col.  Ramsay 
LtCol  Christensen 
Dr.  Bauer 
Maj.  Ttxnlinson 
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;;j°  AFGI/LTT 

.EJECT  Second  Annual  Tri-Service  Cloud  Madeling  workshop 


12  Mar  84 


TO  Distribution  List 
Dear  Colleague 

1 .  This  is  an  announcement  and  a  call  for  papers  for  the  Second  Annual 
Tri-Service  Cloud  Modeling  Workshop,  which  will  be  held  at  the  Naval  Surface 
Weapons  Center,  White  Oak,  MD  on  26-28  June  1984. 

2.  The  first  workshop,  held  at  IDA  in  June  1983,  was  attended  by  about  80 
representatives  of  the  three  Services,  Universities,  and  Industry.  The 
presentations  from  that  workshop  are  contained  in  IDA  Record  Document  D-16 
(Vol.  I  Unclassified,  Vol.  II  Secret)  and  is  available  to  Government  Agencies 
by  request  throLgh  Dr.  E.  Bauer,  IDA  (Institute  for  Defense  Analyses, 

1801  N.  Beauregard  St.,  Alexandria,  VA  22311),  and  to  others  by  request 
through  DARPA. 

3.  The  workshop  to  be  held  26-28  June  1984  will  emphasize  Government  needs 
Lor  cloud  modeling  and  will  examine  those  techniques,  methodologies,  and  data 
bases  that  are  available,  or  under  developnent,  which  can  be  used  in  evaluating 
cloud  inpacts  on  DoD  systems.  We  are  therefore  calling  for  presentations  in 
the  following  subject  areas: 

a.  Cloud  data  requirements  for  high-energy  laser  weapons  systems, 
surveillance,  IRST,  ccnminication  systems,  etc. 

b.  Cloud  modeling  techniques  including  simulation,  3nall-scale  structure, 
IR  backgrounds,  cloud  edge  effects,  cloud-free  line  of  sight,  etc. 

c.  Cloud  data  bases  which  may  be  available  or  applicable  to  the  CLOUDS 
program,  such  as  global  data  bases  (3D/Kr  NEPH,  International  Satellite  Cloud 
Climatology  Program,  etc.),  and  regional  or  climatic  data  bases. 

4.  Please  submit  unclassified  abstracts  in  sufficient  detail  (200-400  words) 
to  provide  reviewers  with  a  basis  for  conparison  and  acceptance.  Indicate  in 
your  abstract  the  classification  of  the  presentation,  time  requirements  and 
the  need  for  special  projection  facilities.  Please  send  abstracts  to  Donald 
D.  Grantham,  AFGL/LYT,  Hanscom  AFB,  MA  01731  by  20  ^ril  1984.  Authors  will 
be  notified  of  the  presentations  by  18  May  1984. 
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5.  The  overall  classification  of  this  meeting  will  be  Secret/No  Foreign 
Nationals.  Clearances  should  be  sent  to  Visitor  Control,  Security  Office, 
Institute  for  Defense  Analyses,  1801  N.  Beauregard  St.,  Alexandria,  VA  22311 
(Tel.  (703)  845-2290  or  AV289-2063)  by  15  June  1984.  The  visit  request  should 
specify  that  the  purpose  is  to  attend  the  2nd  Tri-Service  Cloud  Modeling 
Workshop . 

6.  A  Record  Document  will  be  prepared  and  distributed  as  soon  as  possible 
after  the  meeting.  Speakers  are  asked  to  bring  reproducible  copies  of  their 
vtgraphs  and  a  title/abstract  page  to  the  meeting. 

7.  Arrangements  have  been  made  to  provide  time  and  facilities  for  Government 
Technical  Direction  meetings  and  reviews  before  and  after  the  Workshop.  If 
your  agency  would  like  to  take  advantage  of  these  arrangements,  please  contact 
me  (AV478-2982  or  (617)  861-2982)  . 

Sincerely 


DCNALD  D.  GRANTHAM 

Chief,  Tropospheric  Structure  Branch 
Atmospheric  Sciences  Division 


DEPARTMENT  OF  THE  AIR  FORCE 
AIR  FORCE  GEOPHYSICS  LABORATORY  (AFSC) 
HANSCOM  AIR  FORCE  BASE.  MASSACHUSETTS  01731 


ntPiY  TO 

ATTsof  LYT/(ol7)  861-2982/AV  A78-2982  31  May  1984 

SUBJECT  Second  Tri-Service  Cloud  Modeling  Workshop,  26-28  June  1984 


TO  Distribution  (Ref  Atch  to  AFGL  Ltr  12  flar  84) 

1.  aased  on  the  enthusiastic  response  to  our  12  March  1984  call  for  papers 
for  the  subject  Workshop,  we  have  developed  the  enclosed  tentative  agenda. 

In  general,  presentation  times  will  be  20  minutes  including  discussion.  The 
few  exceptions  to  this  time  limit  have  been  notified  individually.  The  Work¬ 
shop  will  be  held  at  the  Naval  Surface  Weapons  Center  in  the  Auditorium  of 
the  Administration  Building  (see  NSWC  enclosure)  beginning  at  0830  each  day. 

2.  Tne  overall  classification  of  the  Workshop  will  be  Secret/No  Foreign 
Nationals.  Security  clearances  are  being  handled  through  IDA.  If  you  have 
not  yet  submitted  your  clearance,  please  send  it  to  Visitor  Control,  Security 
Office,  Institute  for  Defense  Analyses,  1801  N.  Beauregard  St.,  Alexandria, 

VA  22311  (Telephone  (703)  843-2183)  by  13  June  1984,  noting  that  it  is  for 
attendance  at  the  Tri-Service  Cloud  Workshop  at  NSWC.  If  you  have  a  current 
security  clearance  at  IDA,  please  either  write  or  phone  the  Security  Office  to 
notify  them  that  you  will  be  attending  the  Tri-Service  Cloud  Modeling  Workshop. 

3.  A  record  document  will  be  prepared  and  distributed  as  soon  as  possible 
after  the  meeting.  Speakers  are  asked  to  bring  reproducible  copies  of  their 
vugraphs  and  a  title/abstract  page  to  the  Workshop. 

4.  NSWC  White  Oak  does  not  have  military  quarters  available.  Those  traveling 
on  military  orders  may  pick  up  a  Certificate  of  Non-Availability  of  Quarters 
from  the  Travel  Office  located  in  the  Administration  Building.  A  list  of  motel 
accommodations  is  included  on  the  enclosed  map  showing  their  locations  with 
respect  to  NWSC  White  Oak.  Their  will  be  no  bus  transportation  provided 
between  motels  and  NSWC;  therefore,  car  rentals  may  be  desirable. 

3.  Coffee  and  doughnuts  will  be  available  each  morning.  Lunches  will  be 
available  in  the  Administration  Building  cafeteria.  Plans  are  being  made  for 
an  informal  (beer  and  soda)  party  from  1700-1830  on  20  June. 

0.  The  Cloud  Modeling  Workshop  will  conclude  on  Wednesday,  27  June  with  a 
Wrap-up  Session  at  1300-1600.  Immediately  following  the  Wrap-up,  Dr.  Ernie 
Bauer,  IDA,  will  introduce  an  extra  topic  on  Nuclear  Winter  that  may  be  of 
interest  to  many  of  the  workshop  attendees.  A  Topical  Review  on  this  subject 
will  be  held  the  next  morning,  0830-1200  Thursday,  28  June  1984. 


3  Atchs 

DONALD  D.  GRANTHaM  1.  Agenda 

Chief,  Tropospheric  Structure  Branch  2.  NSWC  Admin  Bldg 

Atmospheric  Sciences  Division  3.  NSWC  Area  Map 
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SECOND  TRI- SERVICE  CLOUD  HODELING  WORKSHOP 
26-28  June  1984 

Navel  Surface  Weapons  Center 
Adalnlatration  Building  Auditorium 
White  Oak,  Silver  Spring,  Maryland 

AGENDA 

Tuesday,  26  June  1984 

Session  I  -  Workshop  Overviews,  0830-1000  -  Dr.  Bemaupd  Kessler,  NSWC, 

Sessicjn  Chairman 

NSWC  Welcome  -  Capt.  J.  J.  Brannan,  USN,  Officer  in  Charge,  White  Oak 
CLOUDS  Program  -  Lt.  Col.  Vernon  Bliss,  AFWL 

Cloud  Model  Methodologies,  Overview  and  Caveats  -  Donald  D.  Grantham,  AFGL 

Cloud  Data  Bases,  Overview  and  Caveats  -  James  Bunting,  AFGL 

Session  II  -  Cloud  Model  Applications,  1020-1230  -  Dr.  Paul  Twitchell,  Naval 
Air  Systems  Oomuutd,  Session  Chairman 

•  "DARPA’s  Interest  in  Cloud  Data,"  David  Zinrerraan,  PRA 

*  "Evaluation  of  the  Effect  of  Clouds  on  SBL  Applications," 

Ronald  J.  Nelson,  SAl 

**  "Effect  of  Cloud  Clutter  on  Target  Detection,"  James  L.  Griggs,  Jr.,  SAl 

"Current  Cloud  Modeling  Efforts  at  the  Naval  Environmental  Prediction 
Research  Facility,  Paul  M.  Tag,  NEPRF 

"  "Trident  Targeting  and  Test  Flight  Analyses,"  Susan  Masters,  NSWC 

**  "Cloud  Effects  on  Surveillance  from  Space  and  High  Altitudes,"  Ernest 
Bauer  and  Hans  Wolf hard,  IDA 

Session  III  -  Cloud  Models  for  IRST  Applications,  1330-1700  -  Dr.  Elton  kvara,  ASL, 
Session  Chaiman 

"Synthetic  Cloud  Scenes  for  IRST  Sensor  Design  Applications,"  Alex  T. 
Maksymowlcz,  Lockheed,  and  John  H.  Allen,  SRI 

"IR  Radiative  Properties  of  Cirrus  and  Cumulus  Cloud  and  IRST  System 
Application,"  G.  Gal  and  T.  Wlnarske,  Lockheed 

"Infrared  Cloud  Backgrounds  and  Sensor  Performance,"  W.J.  Tropf, 

A.N.  Vareck,  Johns  Hopkins  U.,  B.P.  Sandford,  J.H.  Schummers,  AFGL, 
and  J.  Schroeder,  ONTAK 

"The  Navy  Background  Measurement  and  Analysis  Program  (BMAP),  FY84," 

Bernard  V.  Kessler,  NSWC 

"An  Empirical  Cloud  Model  for  SWIR  and  LWIR  Pass  Bands,"  Mike  Scarborough 
and  Robert  A.  Pilgrim,  Teledyne  Brown 

*Paper  classified  SECRET.  Briefing  charts  Included  herein  are  UNCLASSIFIED. 

**Paper  and  briefing  charts  classified  SECRET.  These  papers  are  presented 
in  Volume  11. 
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Wednesday,  27  June  1964 

Session  IV  ~  Cloud  Model  Methodologies.  0830-1100  -  Dr.  Roland  Nagel,  NEPRF, 

Session  Chairman 

'Assessment  of  the  L0WTRAN6  Cirrus  Model,'  John  Hornsteln,  NRL 

'Profiles  of  Optical  Extinction  Coefficients  Calculated  from  Droplet 
Spectra  Observed  In  Marine  Stratus  Cloud  Layers,'  V.  Ray  Noonkester,  MOSC 

'Interactive  Software  for  Calculating  Cloud-Free  Intervals,'  Lt.  Col.  Vernon 
Bliss,  AFWL 

'3-D  Cloud  Simulation  using  the  Sawtooth  Model,'  Irving  I.  Grlngorten,  AFGL 

'Correlation  Functions  for  Cloud  Cover  and  Celling,'  Ralph  Shapiro,  SASC 

'Cloud  and  Visibility  Modeling  of  Joint  Mesoscale  Probability  in  Central 
Europe,'  Oskar  H.  Essenwanger,  US  Army  Missile  Command 

'The  PRA  Synthetic  Cloud  Scenes,'  Edward  J.  Stone,  NRL 

Session  V  -  Cloud  Data  Bases.  1100-1230/1330-1445  -  L/C  Ed  Tonlinson,  AVC, 

Session  Chaiman 

"The  Real-Time  Nephanalysls  (RTNEPH),"  Ma j .  S.  Cox,  AFGWC 

"Cloud  Data  Bases  from  the  International  Satellite  Cloud  Climatology  Project 
(ISCCP),"  Herbert  Jacobowltz,  NOAA 

'Multi-Year  Global  Cloud  Data  Set  From  NIMBUS-?  Satellite,  P.H.  Hwang,  NASA, 
L.L.  Stowe,  NOAA,  and  P.K.  Bhartla,  SASC 

'Viewing  Angle  Bias  In  Cloud  Climatologies,"  James  Bunting,  AFGL 

"Digital  Imaging  and  Analysis  from  Cloud  Climatology  Data,"  T.  Vonder 
Haar  and  T.  Brubaker,  CSU 

"A  New  Data  Base  of  Cloud  Variables  at  Subfreezing  Temperatures,"  Richard 
K.  Jeck,  NRL 

Session  VI  -  Workshop  Wrap-Up,  1500-1630  -  Panel:  Dr.  B.  Kessler,  Dr.  E.  Bauer, 

L/C  V.  Bliss,  Dr.  Elton  Avara 

Thursday,  28  June  1984 

TOPICAL  REVIEW  -  "Cloud  Physical  Aspects  of  Nuclear  Winter"  -  Dr.  Leon  Wittwer,  DNA 

Chairman 

"The  Smoke  Plume  Rise  and  Scavenging  Problems,"  Ernest  Bauer,  IDA* ** 


*Paper  classified  SECRET.  Briefing  charts  included  herein  are  UNCLASSIFIED. 

**Coples  of  the  briefing  are  available  on  request  from  IDA. 


Subject:  Report  on  the  Second  Tri-Service  Cloud 

Modeling  Workshop,  26-27  June  1984 

From:  Steering  Committee:  D.  Grantham  (Chairman), 

E.  Avara,  E.  Bauer,  P.  Twitchell 

This  is  a  maturing  specialty  meeting,  as  can  be  seen  in 
the  historical  context: 

Oct.  1981  Initial  Space-Based  Laser  (SBL)  meeting. 

Presentations  by  systems  analysts  and  meteorologists, 
but  n_o  bottom  1  ine . 

June  1982  SBL-  first  presentation  of  Zero-Order  Assessment 
(see  IDA  B-6*,  Feb.  83).  Things  are  beginning  to 
jell.  Col.  Try  sees  the  need  for  a  coordination 
process  which  led  to  Dr.  Edith  Martin's  Aug.  82 
letter.  Copy  attached. 

Nov.  1982  Tri-Service  Cloud  Modeling  Workshop  planning 

meeting  at  AFGL,  which  is  Lead  Lab.  (see  D-16**, 
p  .  9) . 

June  1983  First  Tri-Service  Meeting  at  IDA.  152  invitees, 

72  attendees.  Successful,  but  too  big  for  IDA- 
next  meeting  at  NSWC  (B.  Kessler,  host). 

June  1984  Second  Tri-Service  Meeting  at  NSWC.  300  invitees, 

95  attendees.  Third  meeting  will  be  held  at  AFGL 
in  May  1985,  back-to-back  with  Atmospheric  Trans¬ 
mission  Conference. 

This  has  been  a  successful  program  initiation  by  DARPA 
(STO  and  DEO),  OUSDRE,  and  IDA.  The  attendance  at  the  meetings, 
quality  of  presentations,  discussion  and  overall  interest  in  the 
topic  shows  that  it  fills  a  need,  but  it  is  clearly  time  to 
integrate  the  meeting  with  the  Atmospheric  Transmission  Conference 
While  funding  in  the  CLOUDS  area  is  still  inadequate,  some 
funding  is  becoming  available  as  the  utility  of  the  work  is 
getting  recognized. 

*IDA  Annotated  Briefing  B-6,  "Cloud  Effects  on  Space-based 
Laser  Weapons  Systems,"  E,  Bauer  et  al.,  February  1983. 

**IDA  Record  Document  D-16,  "Presentations  at  the  Tri-Service 
Cloud  Modeling  Workshop  held  at  IDA,  June  1983,"  E.  Bauer, 
Editor,  July  1983. 


The  biggest  deficiency  of  the  June  1984  meeting  was  clearly 
the  lack  of  user  presentations  listing  needs,  and  this  was 
recognized  during  the  wrapup  session.  We  call  for  a  survey  of 
users'  needs  and  a  marketing  effort  overseen  by  the  steering 
panel . 

(Recall  that  the  CLOUDS  program  concept  was  designed  to 
start  with  users*  needs,  which  call  for  tailored  ou t pu t  to  use 
d e s igna t ed  approaches  to  pr oc es s/u t il ize  relevant  meteorological 
data .  This  concept  of  a  customer-driven  program  as  developed  by 
the  original  steering  group  clearly  remains  va'lid  and  needs 
implementation)  . 

Technical  issues  addressed  in  the  discussion  include  the 
following : 

1.  RV  issues-  if  it  is  a  problem,  it  should  be  significant  for 
both  USAF  and  USN 

2.  Interaction  between  military  and  civilian  agencies  (e.g.,  ISCCP) 
should  be  encouraged,  and  available  data  should  be  advertized  for 
all  relevant  users. 

3.  Important  objects  of  the  meeting  include  both  critical  review 
of  existing  work  (e.g.,  CFLOS,  cloud  radiation  models/data,  cirrus 
models/data  (e.g.,  LOWTRAN  6  Cirrus))  and  the  avoidance  of 
duplication  (e.g.,  in  model  development) 

4.  As  now  constituted,  the  steering  group  represents  scientists/ 
developers  rather  than  users  who  have  requirements  and  funding. 

This  is  a  fundamental  weakness  which  must  be  rectified  if  the 
Tri-Service  Workshop  is  to  be  successful. 

5.  Cirrus  is  clearly  an  area  which  is  ripe  for  development. 

6.  Small-Scale  Structure  needs  emphasis: 

(a)  what  is  needed? 

(b)  how  should  data  be  analyzed /proc es sed ? 

(PSD  vs  third  differences,  Sobel  transform,  etc.) 

(c)  what  are  the  data?  (how  sharp  are  edges  and  structure 

at  different  wavelengths?  -  see,  e.g.,  Z.  G.  Sztankay's 
letter  in  Part  V. 

7.  Validation  of  models  used  by  systems  analysts:  the  meteoro¬ 
logical  community  will  always  be  unhappy  with  the  data  and  models 
used,  but  will  get  system  development  funds  only  if  they  can  show 
that  the  difference  is  critical.  Typ ically.  Lab  Director's  funds 
and  other  6.1  and  6.2  money  will  have  to  be  used  to  address 
validation-  a  vague  feeling  of  unhappiness  won't  extract  system 
development  funds. 


A  topic  that  most  Workshop  attendees  agreed  would  be  an 
asset  to  the  user  community  was  the  establishment  of  a  Clearing 
House  which  would  be  a  tri-service  focal  point  for  cloud  data 
bases  and  methodologies.  This  group  could  provide  background 
information  for  "newcomers"  and  guide  users  to  the  data  bases 
and  techniques  most  applicable  to  their  requirements.  This 
would  include  such  items  as  caveats,  computer  formats,  etc. 

The  tri-service  committee  will  initiate  action  on  this  item 
during  the  next  year. 
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Figure  21.20 

(a)  Frequency  function  for  arbitrary  discrete  distribution  of  Z.  (b)  Cumulative  distnbuiioa 
function  corresponding  to  (a).  It  illustrates  choice  of  a  random  number  drawn  from  the  arbitrary 
distribution  of  Z. 


CUMULATIVE  CLOUD  COVER  DISTRIBUTION  FOR  WHITE  SANDSCOOMST) 


Tenths  of  Cloud  Cover 


length  (km) 

Figure  4.  Exponential  clear  run* 


j  =  mean  length 

’  =  probability  CFLOS 

2  =  probability  no  CFLOS 

P  +  Q  =  1 

MEAN  CLEAR  LENGTH 

y  (p)  =8P+1.2 

p  (P)  =  (P/Q) (8Q  +  1.2)  .5<P<1.0 

MEAN  CLOUDY  LENGTH 

P(Q)=8Q+1.2 

U  (0)  =  (Q/P)  (8P  +  1.2)  .5<Q<_1.0 


INTERVAL  LENGTH  as  a  FUNCTION  OF  PCFLOS 


PCFLOS 


INTERVAL  from  SOLUTION  OF  CUM»JLATIVE  EXPONENTIAL  DISTRIBUTION 
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Probab! I 


CATIONS  OF  DATA  SET  (U) 
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LUND-SHANKLIN  CFLOS  fflPEL 


I  lIND-SHANKLIN  CFLOS  MODEL 


ITATION  AVERAGE  PROBABILITIES  OF  CLFOS  DETERMINED 
AIRCRAR  OBSERVATIONS  AT  A  DEPRESSION  ANGLE  OF  6i 


CLOUD  FREE  FIELDS  OF  VIEW  fPDEL 


LUND,  ET  AL.,  JAM,  VOL  19,  1980 


probability  that  a  10»,  50®,  90®,  130®  AND  170®  FOV  AT 

ZENITH  KILL  HAVE  A  CLOUD-FREE  FRACTION  WHEN  THE  OBSERVER 
REPORTS  5/10  SKY  COVER 
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NEEDS  BETTER  VERTICAL  CORRELATION  OF  CLOUD  STRUCTURE 


MALICK,  ET  AL,  SPIE  VOL  195,  1979 


observed  total  sky  cover,  in  tenths. 


SRI  MODEL  FOR  CFI 


CALIBRATED  ANALYTICAL  MODELING  OF  CLOUD-FREE  INTERVALS 


AFWL  CLOUDS  1(SCENE  GENERATOR)  AND  CLOUD  IMPACT  MODEL  ^(CIM) 


NELSON,  R.J.  AND  LEWIS,  M- ,  1979:  AFWL  TR  79-108 


AFWL  CLOUDS  1(SCENE  GENERATOR)  AND  CLOUD  IMPACT  MODEL  2(CIM)  -  Continued 


CLOUD  DATA  BASES:  OVERVIEW  AND  CAVEATS 


J.  Bunting 

Air  Force  Geophysics  Laboratory 


CLOUD  DATA  BASES:  OVERVIEW  AND  CAVEATS 


26  JUNE  1984 


SOURCES  OF  CLOUD  DATA 


CLOUD  OBSERVATIONS  FROM  WEATHER  STATIONS 


HUMAN  INTERPRETATION 


WEATHER  STATION  SOURCES 


CLOUD  AND  CFL08  ATLASES 


RADIOSONDE  OBSERVATIONS 


rOO  FEW  REPORTS  FOR  QOOD  CLOUD  MAPS 


AIRCRAFT  OBSERVATIONS 


NO  STANDARDS  FOR  OBSERVINQ  OR  ARCHIVINQ 


SATELLITE  CLOUD  CLIMATOLOGIES 


ISCCP 


BEWARE  OF 


EFFICIENT  DATA  HANOLINO  FOR  LARQE  DATA  BASES 


HI -CAMP,  TEAL  RUBY) 


AVD  CLOUD  CLUTTER  DATA  BASE 


CLOUD  CLUTTER  ISSUES 


AVD  CLOUD  SCENE  SIMULATION 


Cloud  Scene  Realizations 
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AIRCRAFT  FLYTHROUGH  TRACKS 


LWIR  FLYTHROUGH  CONTRASTS 
HIGH  CLOUD  SCENE 


DARPA  CLOUD  INTERESTS 
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Evaluation  of  the  Effect  of  Clouds  on  SBL  Applications 

Ronald  J.  Nelson 
Science  Applications,  Inc. 

Abstract 


Conventional  wisdom  has  long  held  that  the  atmosphere  and  its 
phenomena,  and  clouds  in  particular,  could  act  to  reduce  the  potential 
utility  of  space  based  laser  (SB,.)  weapon  systems  in  defensive  or  offensive 
roles  against  endo- atmospheric  targets.  In  1980  a  study  was  commissioned 
with  the  objective  of  determining  the  various  ways  in  which  clouds  could 
effect  SBL  applications  against  long  range  bombers.  The  results  showed 
that  the  problem  was  complicated  but  that  the  magnitude  of  the  problem 
could  probably  be  scoped  using  some  reasonable  -  and,  more  importantly, 
traceable  -  assumptions.  This  paper  presents  the  status  of  a  follow  on 
to  the  1980  study. 

Individual  clouds,  like  the  atmosphere  itself,  scatter  and  absorb 
radiation  passing  through  them.  However,  clouds  act  almost  as  discontinuities 
in  an  otherwise  relatively  homogeneous  medium  so  that,  for  many  practical 
purposes,  clouds  can  be  treated  as  if  they  were  infinite  attenuators.  In 
this  sense,  it  is  sufficient  to  know  whether  or  not  a  cloud  is  in  the  line 
of  sight  between  a  potential  target  and  an  SBL.  However,  clouds  are  also 
sources  of  radiation.  The  radiation  coming  from  individual  cloud  elements 
is  in  sharp  contrast  to  that  coming  from  areas  immediately  adjacent  to  the 
boundary  of  an  individual  element.  Basically,  the  background  radiance  is 
blocked  by  the  cloud  element  (assuming  infinite  attenuation)  which  makes 
its  own,  generally  lower  intensity,  contribution  to  the  total  energy  impinging 
on  the  SBL  sensor  aperture. 

The  principal  objective  of  this  effort  is  the  development  of  a 
quantitative,  first-order  estimate  of  the  effect  of  clouds  on  the  potential 
utility  of  a  space  based  laser  weapon  system  used  against  a  specific  set  of 
airborne  and  ground  based  targets  taking  into  account  the  effect  of  clouds 
on  the  overall  weapon  system,  i.  e.,  on  both  the  laser  beam  and  the  sensor. 
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The  presentation  will  provide  an  overview  of  the  logic  used  to 
accomplish  the  ojbective,  including  a  review  of  the  logic,  analysis 
methodology  and  data  analysis  accomplished  to  date.  This  is  a  companion 
paper  to  the  one  being  presented  by  Dr.  J.  L.  Griggs. 

The  presentation  is  classified  SECRET.*  Time  requested  is  20  minutes 
not  including  time  for  questions  and  answers. 


*0nly  the  UNCLASSIFIED  briefing  charts  presented  at  the  workshop  are  included  here. 
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OBJECTIVE 


Science  Applications.  Inc 


PROGRAM  FLOW  -  TI-83-13 


AIRBORNE/MARITIME 
GROUhO/CORT I  RENTAL 


REGION  Of-  INTEREST 


CONCEPTUAL  LASER  EFFECTIVENESS  MODEL 


KILL  ASSESSMENT 


Science  Applications,  Inc 


CLOUD  INFORMATION  SOURCES 


Science  Applications.  Inc 


DATA  SELECTION  PROCESS 


CFS  SPACE  AS  A  FUNCTION  OF  HEIGHT  CATEGORY 
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FIRST-ORDER  ASSESSMENTS 


115 


Science  Applications,  Inc 


GUARANTEED  SUCCESS"  ~  CLEAR  SKIES 


Science  Applications.  Inc 


GUARANTEED  FAILURE"  -  SINGLE  LAYER  OVERCAST 
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CLOUD  MODELING  EFFORTS  AT  THE  NAVAL 
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ABSTRACT 

Deterministic  synthetic  cloud  scenes  in  the  8-12  micrometer  band 
have  been  generated  for  the  side-looking  viewing  geometry 
characteristic  of  infrared  search  and  track  (IRST)  scenarios. 
LANDSAT  data  were  used  as  the  point  of  departure  for  generating 
three-dimensional  cloud  morphology  and  liquid  water  content 
distributions.  The  integrated  liquid  water  path  was  computed  for 
each  pixel  in  a  256  x  1024  pixel  scene  by  projecting  the  line  of 
sight  through  the  three-dimensional  cloud  distribution.  The 
liquid  water  path  , in  turn,  determined  the  value  of  the  cloud 
radiance,  using  a  model  derived  from  a  detailed  multiple 
scattering  treatment  of  radiative  transfer  in  clouds. 

A  five-frame  temporal  sequence  was  generated  for  each  of  two 
scenes,  a  high  cirrus  layer  and  a  cumulus  layer.  Each  sequence 
corresponds  to  specific  sensor  motion,  stare  point,  and 
interframe  interval.  Each  frame  represents  a  different  projected 
view  of  the  three-dimensional  cloud  distribution.  Thus,  the 
sequence  displays  differential  motion  effects  between  frames, 
wherein  displacement  of  features  depends  on  their  distance  from 
the  sensor.  This  allows  realistic  simulation  of  sensor  motion 
effects  on  signal  processing  operations. 

Descriptions  are  given  of  tne  viewing  scenarios,  sampling 
parameters,  and  the  statistical  properties  of  the 
three-dimensional  cloud  and  radiance  distributions  for  each  of 
the  two  scenes.  To  illustrate  the  application  of  the  synthetic 
scenes  to  IRST  sensor  design  efforts,  pictorial  examples  are 
presented  of  these  scenes,  together  with  their  focal  plane  images 
during  various  stages  of  a  signal  processing  simulation. 
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SYNTHETIC  CLOUD  SCENES  FOR 
IRST  SENSOR  DESIGN 
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IR  Radiative  Properties  of  Cirrus  and 
Cumulus  Clouds  and  IRST  System  Application 

G.  Gal  and  T.  Winarske 
Lockheed  Palo  Alto  Research  Laboratory 

A  serious  problem  In  the  IR  radiometric  measurements  from  an  alrboume 
platform  Is  the  Interpretation  of  the  target  signal  through  the  atmosphere 
and  the  Interference  of  clouds  particularly  the  globally  distributed  cirrus 
layers  or  local  cumulus  clouds.  For  various  system  applications  such  as 
HI-CAMP,  IRST  deterministic  cloud  radiance  distribution  must  be  generated. 

Our  goal  was  to  assemble  a  computer  model  to  calculate  cloud  spectral  radiance 
and  atmospheric  transmittance  based  on  realistic  cloud  physical  parameters  and 
various  Irradlance  models  such  as:  direct  emission  Including  multiple  scat¬ 
tering  effects,  scattering  of  diffuse  earthshine,  upwelllng  and  downwelling 
alrshlne,  scattered  sun  and  solar  radiance. 

Furthermore  a  first  order  attempt  has  been  made  to  obtain  an  estimate  of 
the  effects  of  a  nearby  cloud  (mutual  coupling)  and  cloud  edge  radiance  values. 
Evaluation  of  these  contributions  require  the  solution  of  the  radiative  trans¬ 
fer  equation  for  a  coupled  sky-cloud-solar-earth  system  for  the  particular 
geometric  configuration.  Neither  physical  input  nor  the  computational  metho¬ 
dology  Is  currently  available.  Our  approach  was  to  break  down  the  problem  to 
the  solution  of  well  defined  model  problem  to  obtain  cloud  Internal  radiation 
field  for  given  cloud  model  parameters  and  obtain  some  quantitative  estimate 
of  the  relative  Importance  of  the  various  irradlance  contributions. 

A  plane  parallel  atmosphere  with  horizontal  extent  Is  treated  as  being 
made  up  of  stratified  layers.  The  cloud  layer  Is  located  at  an  arbitrary 
altitude  divided  into  small  sublayers.  Clear  air  extent  below  and  above  the 
cloud.  The  cloud  layer  consists  of  clear  atmospheric  molecules  and  arbitrary 
droplet  size  distribution  described  by  modified  gamma  distribution  or  similar 
continuous  function.  Altitude  dependent  meterologlcal  and  mlcrophyslcal  are 
input  data.  These  are  air  pressure,  temperature  dew  temperature,  droplet 
concentration,  etc.  Radiative  transport  equation  for  a  monochromatic  radiation 


through  thus  Inhomogeneous  scattering  and  absorbing  layers  which  assumed  to 
be  In  local  thermodynamic  equilibrium  Is 

dI,(T;0,(t.) 

W  - -  a)^(T)J^(T;0,()») 

Where  I  and  J  Is  the  spectral  Intensity  and  source  term  respectively;  u.Ct) 

A  A  A 

single  scattering  albedo;  u  Is  the  cosine  of  the  zenith  angle  0;  41  Is  the 
azimuth  angle.  The  optical  depth  t  Is  directly  obtained  from  rigorous  Mle 
theory  of  scattering  from  the  assumed  cloud  particle  size  distribution. 

Similarly a  the  cloud  averaged  scattering  phase  function  Is  obtained  as  a  sum 
of  the  Rayleigh  (molecular)  and  Mle  (particulate)  components.  Utilizing  the 
Lengendre  and  Fourier  series  expansion  of  the  phase  function  the  solution  of 
the  radiative  transfer  equation  Is  obtained  by  a  numerical  method  outlined  by 
Dave*Gazdag  with  modification  Introduced  by  Low. 

For  our  Cloud  Clutter  Model  Development  program  our  Interest  has  been 
focused  In  the  8  £  X  (pm)  ^  12  spectral  band.  Spectral  up/down-welllng 
radiance  values  are  obtained  as  a  function  of  Nadir  angle  for  an  Azimuthal 
symmetric  case  (sun  has  been  omitted  for  this  particular  spectral  region)  for 
a  cirrus,  cumulus  and  cumulus  edge  with  or  without  atmospheric  layer  above  or 
below  the  cloud  layer. 

This  model  solution  Is  the  point  of  departure  for  constructing  engineering 
models  for  different  system  applications,  such  as  the  shallow  angle  viewing 
IRST  study.  Particular  details  of  the  cloud  scenes,  both  physical  and  morpho¬ 
logical,  are  taken  Into  account  by  this  simplified  engineering  model.  Estimates 
have  been  obtained  for  radiance  values  for  multilayer  clouds,  cumulus  edge 
effects,  as  well  as  for  the  different  source  Irradlance  contribution  to  the 
total  radiance  values.  This  Is  particularly  Important  for  very  shallow  viewing 
angles  where  the  LOS  does  not  Intercept  the  earth  surface. 

For  a  general  data  analysis  and  background  study  this  computer  model  • 
may  be  applicable  upon  completltlon  of  the  numerical  solution  for  other 
spectral  regions.  It  would  be  especially  Important  to  compare  the  prediction 
of  cumulus  edge  radiance  values  with  recently  obtained  experimental  data. 


RADIANCE  MODELING 


ENGINEERING  MODEL;  DETERMINISTIC,  PHYSICS-BASED  ALGORITHMS  FOR 

GENERATING  SYNTHETIC  CLOUD  SCENES 


SOURCES  OF  RADIANCE  IN  THE  ATMOSPHERE 


SURFACE  RADIANCE 
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DUE  TO  EARTH  CURVATURE.  DIRECT  EARTH  RADIANCE  MAY  BE  OMITTED 
FOR  SCENARIO  I  AND  II  MODELING 


SPECTRAL  RADIANCE  DUE  TO 
ATMOSPHERE  PLUS  EARTH 


ANCLE  FROM  NADIR  (deg) 


OPTICAL  DEPTH  OF  A  HOMOGENEOUS  CLOUD  LAYER 
(U.S.  STANDARD  ATMOSPHERE) 
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SPECTRAL  RADIANCE  AT  H  =  20  km 
FOR  CIRRUS  LAYER 


ANCLE  FROM  NADIR  (deg) 
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CIRRUS  LAYER  EMISSIVITY  AT  A  =  10.5  urn 
AS  A  FUNCTION  OF  CLOUD  DEPTH 
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ANGULAR  VARIATION  OF  THE  UPWELLING 
HOMOGENEOUS  STRATOCUMULUS  CLOUD  LAYER 
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ISOLATED  CUMULUS  CLOUD  LAYER 
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SPECTRAL  RADIANCE  FOR  CUMULUS  CLOUD  LAYER 


SPECTRAL  RADIANCE  AT  H  •  10  km  FOR 
ATMOSPHERE  PLUS  CUMULUS  LAYER 
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OPTICALLY  THICK  CUMULUS 
(NOT  YET  STUDIED) 


NADIR  ANCLE  9  (deg) 


FROM  THE  ITERATIVE  EXACT  SOLUTION 
(MULSCT)  CODE 


WAVELENGTH 


COMPUTATIONAL  ’ENGINEERING” 
MODEL  FOR  fs  sA(um)  ^121  BAND 


AVERAGED  MIE  COEFFICIENTS 


DIFFUSE  CLOUD  RADIANCE  AS  A  FUNCTION 
OF  NORMALIZED  LIQUID  WATER  PATH 
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COMPLETE  NUMERICAL  SOLUTION  FOR  OTHER  SPECTRAL  REGION  (VISIBLE. 
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Abstract 


The  performance  of  an  Infrared  sensor  may  be  greatly  affected  by 
cloud  clutter.  The  Johns  Hopkins  University  Applied  Physics  Laboratory,  as 
part  of  the  Naval  Sea  Systems  Command's  Wide  Area  Guidance  and  Control 
Program,  Is  Investigating  Infrared  backgrounds  and  missile  seeker  signal 
processing.  Clouds  are  a  potentially  significant  clutter  source  for  both 
horizontal  and  down-viewing  sensors*  A  library  of  representative  cloud 
imagery  has  been  assembled  from  data  taken  by  the  Airborne  Measurements 
Branch,  Infrared  Technology  Division  of  the  Air  Force  Geophysics  Laboratory. 
Typical  scenes  include  horizon  scans  with  a  foreground  cloud  deck,  cloud 
formations  (e.g.  cirrus,  alto-cumulus,  cumulo-nimbus,  cirro-stratus,  etc.)  and 
terrain  background  partially  obscured  by  clouds. 

These  data  have  been  analyzed  to  characterize  clouds  In  terms  of 
their  spatial  properties  such  as  size,  intensity,  contrast  and  edge 
sharpness.  This  includes: 

o  Scene  statistics  -  the  mean,  variance,  or  higher  moments,  are  used 
to  determine  the  variation  within  a  scene; 
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o  Power  spectral  density  -  the  frequency  distribution  of  clutter 
power  (spatial  correlation  length  can  be  found  from  the  Fourier 
transform  of  the  power  spectral  density); 

o  Image  processing  and  enhancement  techniques  -  properties  such  as 
slopes,  local  averages,  textures  analysis,  and  edge  detection  are 
examined  to  determine  scene  content. 

o  Bandpass  filtering  analogous  to  that  used  by  scanning  Infrared 
seekers; 

This  paper  will  expand  on  the  results  presented  at  the  1983  Trl~ 
Services  Infrared  Background  Symposium  (MITRE  October  1983).  New  methods  to 
characterize  cloud  Imagery  and  to  predict  sensor  performance  against  clutter 
backgrounds  will  be  presented  and  discussed. 

1.0  Introduction 

Future  Naval  air  defense  missiles  will  Include  an  infrared  guidance 
mode.  Applications  of  Infrared  missile  guidance  In  this  mission  Include: 

1.  avoiding  ECH,  2.  Improving  hlgh-altltude  missile  dynamics,  and 
3.  countering  low  target  radar  cross  section.  The  Naval  Sea  Systems  Command 
Wide-Area  Guidance  and  Control  Exploratory  Development  Program  (conducted  by 
The  Johns  Hopkins  University  Applied  Physics  Laboratory)  has  recently 
Initiated  a  study  of  the  capability  of  advanced  highly  sensitive  Infrared 
seekers  to  acquire  targets  under  high  clutter  conditions.  Long  range 
detection  of  air  vehicle  targets  at  both  high  and  low  altitudes  Is  a  goal  of 
this  work.  The  effort  Includes  both  the  development  of  signal  processing 
techniques  and  the  evaluation  of  signal  processing  against  realistic  Infrared 
backgrounds.  This  paper  Is  the  second  report  of  Initial  results  from  this 


Although  terrain  backgrounds  have  been  extensively  studied  over  the 
past  several  years,  the  effects  of  atmospheric  clouds  on  sensor  performance  Is 
an  area  that  has  not  been  fully  Investigated.  Clouds  formed  by  atmospheric 
water  vapor  and  Ice  can  reduce  target  visibility,  ranging  from  partial 
obscuration  to  total  occultatlon,  and  consequently  reduce  sensor 
performance.  Of  equal  Importance  are  limitations  arising  from  vle«rlng  a 
target  against  a  cloud  background  or  a  mixture  of  clouds,  sky  and  terrain. 

This  paper  discusses  these  limitations  and  presents  analysis  which 
characterize  clouds  In  terms  of  their  spatial  properties.  The  data  used  for 
this  work  were  acquired  by  the  Airborne  Measurements  Branch  of  the  Air  Force 
Geophysics  Laboratory  (AFGL). 

2.0  AFGL  Cloud  Data 

The  Air  Force  Geophysics  Laboratory  operates  an  instrumented  NKC-135A 
aircraft  to  acquire  target  and  background  signature  data.  This  system 
Includes  a  number  of  spectral,  radiometric,  and  spatial  instruments  capable  of 
making  state-of-the-art  Inflight  measurements.  Among  these  is  a  calibrated 
Infrared  mapper  (FLIR)  which  Is  used  to  obtain  absolute  radiance  Imagery  in 
the  8-14  pm  region. 

Since  1978  the  AFGL  has  used  the  FLIR  mapper  to  acquire  signature 
data  of  air  vehicles,  ground  targets,  and  background  scenes.  From  the  wide 
variety  of  available  data  a  digital  tape  of  31  calibrated  scenes  was  prepared 
for  analysis.  These  data  Include  cloud  deck  with  horizon,  cloud  formations, 
and  terrain  backgrounds  partially  obscured  by  clouds.  Table  1  summarizes  the 
basic  characteristics  of  the  scenes  in  terms  of  the  mean  and  standard 
deviation  of  radiance  and  an  approximate  correlation  length. 
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Table  1 


AFGL  Scene  Parameters 


AFGL 

Scene 

Scene  Radiance  (W/cm^/sr)* 

Correlation 

Image 

Designation 

IFOV 

Mean  Deviation 

101 

I0504A 

0.5  mrad 

2.88  X  10“^ 

1.62  X  10"^ 

102 

1051 OA 

3.44  X  10“J 

1.06  X  10“J 

103 

I0804A 

3.25  X  10“J 

3.26  X  10"t 

104 

I0804B 

3.82  X  10“5 

3.46  X  10“^ 

105 

lOlOlA 

3.01  X  10“3 

7.04  X  10"5 

106 

lOlOlB 

2.89  X  lO'l 

1.21  X  10"^ 

107 

lOlOlC 

3.23  X  10~3 

2.37  X  10“^ 

108 

I0103A 

1.74  X  10“3 

2.18  X  10“J 

109 

101 03B 

2.80  X  10“3 

1.11  X  10"* 

110 

10103C 

2.54  X  10“^ 

8.84  X  10"^ 

Ill 

I0901Z 

2.25  X  10“3 

1.23  X  10“* 

112 

I0902A 

2.27  X  10“J 

1.45  X  10"* 

113 

101 02Z 

1.11  X  10~^ 

1.94  X  10"* 

114 

I0104A 

1.54  X  10~3 

2.14  X  10"* 

115 

B0104B 

1.87  X  10“3 

8.08  X  10"5* 

116 

I1201A 

1.55  *  '0“J 

3.02  X  10"* 

117 

I1201B 

2.41  X  10“3 

1.35  X  10"* 

118 

11201C 

3.13  X  IQTI 

2.21  X  10"* 

119 

I1201D 

2.10  X  10“^ 

2.52  X  10"* 

120 

I1201E 

2.41  X  10*J 

8.95  X  10"5 

121 

I2003A 

0.1  mrad 

4.01  X  10”3 

1.61  X  10"* 

1122 

I2003B 

3.98  X  10"3 

1.74  X  10"* 

123 

I2003C 

4.02  X  10”3 

2.05  X  10"* 

124 

I2003D 

3.94  X  10“^ 

1.49  X  10"* 

125 

I2003E 

4.09  X  10"3 

1.64  X  10"* 

126 

I2003F 

4.13  X  10“3 

1.33  X  10"* 

127 

12003G 

4.06  X  10"3 

1.56  X  10"* 

128 

I2003H 

4.28  X  10“^ 

1.34  X  10"* 

129 

I2003I 

3.94  X  10"3 

1.22  X  10"* 

130 

I2503A 

0. 5  mrad 

2.02  X  10~3 

1.86  X  10"* 

131 

I2503B 

1.81  X  10~3 

3.20  X  10"* 

Scene 


Length  (mrad) 

18.5 
12.0 
20.0 

24.5 

6.8 

6.5 

16.0 

18.5 
14.0 

7.3 

12.0 

4.4 

21.7 
19.0 

22.5 
25.0 

13.5 

18.5 

21.0 

17.5 

13.7 

3.5 

4.6 
3.2 
4.1 
3.4 
4.0 

3.6 
3.0 

20.0 

12.0 


*  7.84-13.16  vm  spectral  bandpass  except  Image  115  mhlch  is  oncallbrated. 
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3.0  Cloud  Edge  Characterization 


Cloud  edges  may  be  a  significant  source  of  false  alarms.  Both  the 
extent  and  steepness  of  the  edges  will  contribute  to  the  seeker  tracking 
problem.  Our  previous  paper  (Reference  1)  used  a  series  of  linear  directional 
masks  to  examine  edge  properties  which  showed  that  these  operators  are  not 
adequate  for  delineating  edge  properties.  This  paper  will  use  a  3x3  nonlinear 
Sobel  (Reference  2)  enhancement  operator  to  Investigate  cloud  edge  properties 
for  the  four  images  (Figures  1-4  of  Reference  1).  The  temperature  contours 
for  Images  123  and  124  (Figures  3  and  4  of  Reference  1)  are  given  for 
reference  In  Figures  1  and  2. 

A  Sobel  operator  is  of  the  form: 

T’(1,J)  -  (x^+y^)^/^ 

where 

X  -  (c+2f+i)  -  (a+2d+g) 

T  -  (a+2W-c)  -  (g+2h+i) 

for  a  3x3  neighborhood  of  pixel  T(l,j): 

a  b  c 

d  T(i,j)  f 
g  h  1 

and  the  letters  a,b,  etc.  represent  the  pixel  values  (e.g. ,  radiance)  at  the 
respective  positions  relative  to  pixel  T(i,J).  The  results  of  applying  a 
Sobel  operator  to  Figures  1  and  2  (Images  123  and  124)  are  shown  In  Figures  3 
and  4.  The  range  in  values  for  the  two  images  was  from  slightly  greater  than 
zero  to  24  (Image  123)  and  32  (Image  124).  However  the  mean  and  distribution 
of  values,  as  shown  In  Table  2,  were  similar  for  both  Images. 
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Laukel  marvuino 


Sobel  valiM 


Table  2.  Distribution  of  Sobel  Values 


Value  of  Sobel 
Operator 

Image 

123 

Image 

124 

Scene 

Content 

0-2 

28Z 

29Z 

Intra-cloud 

2-4 

44 

38 

Structure 

4-6 

15 

17 

6-10 

7 

9 

10-20 

6 

7 

Cloud 

>20 

«1 

«1 

Boundaries 

Mean 

4.0 

4.0 

Figures  3a  and  4a  give  the  contours  of  the  tvio  Images  generated  by  the  Sobel 
operator  for  values  2  through  4.  In  each  of  these  figures  all  values  less 
than  the  lower  threshold  are  Indicated  by  the  absence  of  shading  while  those 
above  the  upper  threshold  are  Indicated  by  a  high  density  broken  diagonal 
shading.  Figures  3b  and  3c  for  Image  123,  and  4b  and  4c  for  Image  124  are 
plots  of  the  data  in  Figures  3  and  4  with  different  threshold  values.  These 
plots  more  clearly  show  the  structure  within  the  cloud  formations  as  well  as 
the  magnitude  and  sharpness  of  the  edges.  Both  Figures  3  and  4  show  a  large 
amount  of  edge  content  both  within  and  between  the  cloud  formations.  Within  a 
formation  the  values  are  small  (generally  less  than  4)  and  the  gradients  are 
not  steep.  Of  Interest  Is  that  these  values  are  comparable  to  those  within 
the  terrain  (see  the  upper  left  of  Image  123  and  the  upper  central  portion  of 
Image  124. 


The  steepest  gradients  are,  not  surprisingly,  seen  either  between  the 
cloud  formations  or  between  the  clouds  and  terrain,  as  Is  seen  in  Image  123. 
The  appearance  of  a  thicker  contour  line  in  Figures  3  and  4  for  this  region  is 
the  result  of  several  closely  spaced  contour  levels.  The  values  for  the  edge 
betmen  the  cloud  formations  range  from  4  to  12  while  the  strong  edge  at  the 
top  of  the  scenes  has  values  between  4  and  20. 


4.0  Seeker  Acquisition  Performance  Using  Power  Spectral  Density  Data 


In  a  previous  paper  (Reference  1)  we  described  AFGL  cloud  scenes  with 
one-dimensional  power  spectral  densities  averaged  over  all  scan  lines.  Figure 
3  shows  the  power  spectra  for  two  typical  scenes,  one  sampled  at  0.5  mrad 
intervals  and  the  other  at  0.1  mrad  Intervals.  The  k  (k  -  spatial  frequency 
In  wavenumbers)  dependence  of  clutter  power  is  typical  of  many  backgrounds. 

The  high  concentration  of  noise  at  low  frequencies  indicates  a  preponderance 
of  large  features  In  the  scene.  An  Infrared  seeker  searching  for  point 
targets  exploits  this  frequency  rolloff  of  clutter  power  using  a  processor 
(adaptive  filter)  that  blocks  the  predominately  lower-frequency  scene  noise. 

The  relatively  flat  portion  of  the  power  spectral  density  is  believed 
to  be  the  FLIR  Internal  noise  floor  (see  reference  1).  When  an  afocal 
telescope  Is  used  to  Increase  FLIR  resolution,  the  backgrourd  power  spectrum 
is  unchanged  while  the  Internal  noise  power  spectrum  Is  lower  and  spread  In 
spatial  wavenumber  (l.e.,  the  total  Internal  noise  power  Is  unchanged). 

Hence,  the  higher  resolution  scenes  give  clutter  data  at  higher  spatial 
frequencies. 

In  Reference  1  %re  determined  the  required  target  irradlance  for 
detection  of  point  clutter  In  various  cloud  scenes  using  the  following 
assumptions: 

1)  The  target  signal  Is  an  impulse,  l.e.,  the  target  Fourier 
transform  Is  flat  (to  the  background  cutoff  frequency)  with  a 
magnitude  of  E^Ax  where  Is  the  target  irradlance  and  Ax  Is  the 
spatial  sample  Interval. 


2)  A  fixed  slgnal-to-nolse  ratio  Is  used  as  the  threshold,  l.e., 
(S/N)^  ■  50  (after  filtering)  Is  used  as  an  arbrltary  threshold 
for  detection  and  false  alarm  rate.  (Note  that  a  scene  with 
Gaussian  background  statistics  and  49152  pixels  would  require  a 
threshold  of  approximately  16  to  Insure  a  false  alarm  probability 
of  less  than  0.1). 


Fig.  5  Scene  Power  Spectra. 


These  assumptions  were  used  in  Reference  1  to  determine  the  required 
target  irradiance  for  the  cases  of  an  ideal  bandpass  filter  (passing  all 
frequencies  for  k  to  kmax,  kmax  ■  l/2Ax)  and  a  matched  filter.  Furthermore, 
FLIR  noise-limited  acquisition  range  was  estimated  using  the  flat  portion  of 
the  power  spectra. 

Matched  Filter  Results 


Table  3  summarizes  the  required  target  Irradiance  using  a  matched 
filter  for  30  cloud  scenes.  Three  results  are  given:  (1)  required  target 
irradiance  for  the  scene  data  (background  plus  FLIR  noise),  (2)  required 
target  irradiance  for  FLIR  noise  only  (white  noise  assumed),  and  (3)  required 
target  Irradiance  against  the  background  only  (l.e.,  the  background  power 
spectral  density  is  extrapolated  beyond  the  FLIR  noise  floor).  The  first  two 
results  are  an  extension  of  those  given  in  Reference  1;  the  third  is  new. 

Comparison  of  required  target  irradiance  for  the  scene  data  and  the 
FLIR  noise  level  shows  that  the  FLIR  noise  significantly  affects  the  usefulness  of 
the  scene  data  to  determine  detection  performance  in  backgrounds.  In  order  to 
get  an  estimate  of  the  effect  of  background  alone  on  detection,  scene  data 

power  spectral  density  was  used  for  low  wavenumbers  and  extrapolated  with  a  k~ 

O 

function  to  higher  spatial  frequency.  The  validity  of  this  extrapolation  to  1 
cycle/arad  (adequate  for  examination  of  0.5  arad  scenes)  is  confirmed  by 
examination  of  power  spectral  density  of  high  resolution  scenes  (see  Figure 
5).  This  artlfical  removal  of  FLIR  noise  lowers  the  required  target 
irradiance  for  detection  by  a  factor  of  2. 

Comparison  of  0.5  mrad  scene  results  with  those  from  0.1  arad 
resolution  scenes  shows  the  value  of  high  spatial  resolution  in  clutter 
suppression.  Presuming  a  k  dependence  of  background  clutter,  the  matched 
filter  performance  will  improve  (required  target  irradiance  decrease)  as  the 
square  of  the  sensor  resolution.  Unfortunately,  diffraction  Halts  the 
resolution  achievable  by  small  aperture  seekers.  The  diffraction  limit  of  a 
two-inch  aperture  is  approximately  0.5  arad  at  10  pa  . 
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Table  3 

Required  Target  Irradlance  for  Material 
Filter  Detection  in  Several  AFGL  Cloud  Scenes 


Required  Target  Irradlance 

(10“^^W/cm2) 

Scene 

Background 

FLIR 

Background 

Image 

IFOV  Description 

Plus  Noise 

Noise  Only 

Only 

101 

0.5  mrad 

Down  Viewing  Cloud 

3.98 

2.96 

1.81 

102 

Down  Viewing  &  Water 

3.46 

2.67 

1.53 

103 

Down  Viewing  &  Terrain 

8.36 

6.46 

4.80 

104 

Down  Viewing  &  Terrain 

8.67 

6.82 

4.67 

105 

Clouds  &  Horizon 

3.53 

2.73 

1.73 

106 

Cloud  Deck 

4.97 

3.68 

2.95 

107 

Cloud  Structure 

5.53 

4.22 

3.70 

108 

Clouds  &  Horizon 

9.81 

8.14 

2.94 

109 

Monument  Valley 

4.17 

3.21 

2.34 

110 

Monument  Valley 

4.09 

3.12 

2.23 

Ill 

Close  Clouds 

3.39 

2.71 

1.50 

112 

Close  Clouds 

3. 38 

2.72 

1.63 

113 

Clouds  &  Horizon 

5.29 

4.36 

2.62 

114 

Clouds  &  Horizon 

5.75 

4.65 

3.01 

116 

Horizon  &  Clouds 

6.20 

4.62 

3.26 

117 

Clouds  &  Terrain 

5.20 

4.35 

2.07 

116 

Clouds  &  Horizon 

5.95 

4.57 

3.02 

119 

Clouds  &  Horizon 

5.64 

4.34 

3.36 

120 

Down  Viewing  Clouds 

3.32 

2.69 

1.40 

131 

Clouds  Below  Horizon 

4.54 

3.55 

3.23 

132 

Down  Viewing  Clouds 

4.96 

3.86 

3.70 

21  scene 

0.5  mrad 

5.23 

4.12 

2.74 

Average 

121 

0. 1  mrad 

Structured  Clouds 

.147 

.110 

.070 

122 

ft  ft 

.152 

.107 

.083 

123 

ft  ft 

.175 

.131 

.106 

124 

ft  ft 

.150 

.122 

.072 

125 

ft  ft 

.171 

.134 

.084 

126 

••  M 

.154 

.118 

.077 

127 

M  •« 

.165 

.125 

.084 

128 

■f  •• 

.160 

.127 

.077 

129 

ft  ft 

.144 

.114 

.063 

9  scene 

0. 1  mrad 

.157 

.121 

.080 

average 


One  way  of  testing  the  validity  of  this  extrapolated  data  technique 
is  to  use  the  0.1  mrad  resolution  scenes  with  an  assumed  0.5  mrad  resolution 


sensor  (l.e.,  only  use  the  background  power  spectrum  out  to  1.0  cycle/mrad 
spatial  frequency.)  These  results  are  shown  In  Table  4.  This  9  scene  average 
of  2.54  X  10  W/cm  for  the  required  target  Irradlance  compares  closely  to 
the  value  of  2.74  x  10  M/cisr  obtained  from  an  average  of  21  extrapolated 

scenes  with  0.5  mrad  resolution  (see  Table  3).  Hence  the  extrapolation 
appears  valid. 

Using  characterized  scenes  (l.e.,  power  spectral  densities)  for 
performance  analysis  gives  a  good  answer  for  ultimate  seeker  performance 
provided  the  assumptions  used  are  valid.  Two  principal  assumptions  likely  to 
be  violated  are  that  the  scene  Is  (1)  stationary  and  (2)  obeys  gausslan 
statistics.  The  stationary  criterion  requires  the  scene-average  statistics  to 
be  representative  of  the  entire  scene:  the  spatial  correlation  between  scene 
elements  Is  dependent  only  on  the  spatial  separatioi.  "'tween  the  points  and 
not  In  position  In  the  scene.  Power  spectral  density  cnaracterlzatlon  of  the 
scene  Includes  only  a  measure  of  scene  variance,  not  higher  moments.  Complete 
scene  statistics  are  needed  to  calculate  the  threshold-to-scene  noise  ratio 
required  for  detection  (l.e.,  set  the  false  alarm  rate). 


5.0  Seeker  Acquisition  Performance  Using  Butterworth  Filters 

Required  target  Irradlance  can  be  more  accurately  calculated  by 
directly  filtering  the  scene.  This  approach  Is  straightforward,  gives  false 
alarm  statistics,  but  unfortunately  can  not  separate  FLIR  noise  from  scene 
contributions  In  the  data. 


215 


Table  4 

Required  Target  Irradlance  for  Detection  %d.th 
0. 5  mrad  Resolution  Seeker  Using  0. 1 
mrad  Resolution  AFGL  Cloud  Scenes 


Description 

■red  Clouds 


Required  Target 
Irradlance  W, 


2.48 


mnr 


Filtering  was  performed  digitally,  on  a  llne-by-llne  basis,  using 
Butterworth  high-pass  filters  of  various  orders.  Two  filters  were  of 
particular  interest;  a  first-order  filter,  which  approximates  a  prewhltenlng, 
or  matched  filter;  and  a  sixth-order  filter,  which  has  an  extremely  sharp 
cuton,  and  represents,  approximately,  a  rectangular  filter.  Figure  6  shots 
the  spectral  character  of  these  filters,  and  demonstrates  their  effect  on  a 
typical  image  power  spectrum.  Figure  7  compares  an  image  after  high-pass 
filtering  with  a  sixth-order  filter,  with  the  same  image  after  convolution 
with  the  template 


[-1  2  -1]  . 

Both  images  were  contrast  enhanced  to  roughly  equate  their  means  and  variances 
and  thresholded  at  the  same  level.  The  convolved  image  appears  very  similar 
to  the  digitally  filtered  image,  indicating  the  sharp  high-pass  character  of 
the  sixth-order  filter  and  the  similarity  of  these  methods  in  Isolating  scene 
features. 


Statistical  analysis  of  images  after  filtering  leads  to  a 
determination  of  the  target  Irradlance  required  for  a  constant  false  alarm 
rate.  In  the  case  of  AFGL  Images,  which  contain  only  49,152  pixels,  the  false 
alarm  rate  is  fixed  at  zero,  with  the  detection  threshold  set  Just  above  the 
maximum  pixel  value  in  the  filtered  scene.  As  will  be  shown,  consideration  of 
both  positive  and  negative  thresholds  is  Important,  especially  when  high  order 
filters  are  applied.  A  negative  threshold  is  set  at  Just  below  the  minimum 
pixel  value  of  the  filtered  scene  to  preclude  false  alarms. 

Table  5  presents  required  target  irradlances  for  five  different 
images,  chosen  for  their  variety  of  background  characteristics.  These 
Irradlances  are  computed  with  the  assumption  that  a  positive  detection 
threshold  is  set.  With  a  positive  threshold,  the  first-order  Butterworth 
filter  does  not  require  as  high  a  target  irradlance  for  detection  as  a  sixth- 
order  filter  at  the  same  cuton  frequency.  With  a  negative  threshold  is 
considered,  this  is  not  true.  Table  6  gives  required  target  irradlances 

217 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAuAfL  Maryland 


Itt  order 
("Whitening") 

Butterwofth 

/ 

filter 

? 

/ 

A 

/ 

/ 

Itt  order 
("Whitened") 


Spatul  frequency 


Original  image 
(image  5) 

power  spectral  density 


Spatial  frequency 


Filter  transfer 
function  (power) 


Spatial  frequency 
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Fig.  6  Digital  Filter  Characteristics. 


a.  AFGL  image  5  convolved 


b.  AFGL  image  5  filtered 


F  ig.  7  Image  Processing  Comparison. 


Table  5 


Filter  Performance  Using  Positive  Threshold,  FAR  ■  0 


1st  Order  6th  Order 


Image 

(T/o) 

E^,req(10“^^W/cm^) 

(T/o) 

E^,req{10“^^ 

101 

38.2 

3.85 

28.6 

6.87 

105 

186.1 

7.97 

135.9 

12.77 

106 

75.4 

7.76 

55.5 

12.17 

108 

36.3 

9.40 

17.3 

11.72 

109 

33.3 

4.06 

29.8 

6.63 

Table  6 


Filter  Performance  Using  Negative  Threshold,  FAR  ■  0 


Ist  Order  6th  Order 


Image 

(T/o) 

E^,req(10"^Vcm2) 

(T/o) 

Ej,req(10"^^W/cm^) 

101 

28.1 

6.51 

32.3 

3.43 

105 

103.2 

15.40 

62.0 

4.75 

106 

44.6 

11.06 

67.8 

6.91 

108 

18.5 

18.32 

17.3 

6.86 

109 

22.6 

6.05 

29.8 

3.71 
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assuming  a  negative  threshold  for  detection*  Two  phenomena  contribute  to  this 
effect.  The  first  can  be  seen  If  an  examination  of  the  Impulse  responses  of 
these  filters  Is  made.  Figure  8  Is  a. plot  of  the  response  of  f Irst-and-sixth- 
order  Butterworth  high-pass  filters  to  an  impulse  of  height  100  digital 
counts.  The  filter  cuton  In  both  cases  was  0.25  arbitrary  frequency  units, 
with  the  maximum  frequency  being  1  unit.  The  extremely  sharp  ne^atlve  peak  In 
the  response  of  the  higher  order  filter  gives  rise  to  a  low  Impulse  absolute 
amplitude  loss  after  filtering,  one  comparable  to  the  first-order  filter 
loss. 


The  second  phenomenon  Involves  the  sensor  Internal  noise,  which  Is 
almost  certainly  gausslan  In  nature,  while  the  lower  frequency  background  Is 
definitely  non-gausslan.  As  the  sixth-order  filter  passes  less  of  the  non- 
gausslan  background  than  the  first  order,  the  sixth-order  filtered  scene 
becomes  more  gausslan  and  Its  variance  Is  lower  than  the  first  or  any  lower 
order  filtered  scene.  This  lower  variance  coupled  with  a  low  filter  loss, 
leads  to  improved  performance  on  the  part  of  the  sixth-order  filter. 

In  comparison  with  calculations  earlier  In  this  paper  of 
theoretically  determined  required  target  irradlances.  Table  7  presents  the 
Irradlances  In  Table  5,  given  an  output  slgnal-to-background  noise  ratio  of  50 
In  power.  The  average  required  target  Irradlance  of  these  scenes  Is  slightly 
higher  than  the  level  for  an  ideal  matched  filter,  for  both  the  first  and 
sixth-order  cases.  Here  too,  of  course,  the  sixth-order  filter  requires  a 
lower  target  Irradlance  for  a  zero  false  alarm  rate  than  the  first-order 
filter  with  the  same  cuton  frequency. 

Required  target  Irradlances  have  been  calculated  for  first-  and  sixth- 
order  filters  spanning  a  range  of  cuton  frequencies.  The  result  Is  plotted  in 
Figure  9,  for  AFGL  image  105,  against  the  ideal  bandpass  filter  required 
irradlance  (see  Reference  1).  Required  target  irradlances  for  the  ideal 
filter  case  are  lower  than  those  for  a  realizable  filter,  until  higher  cuton 
frequencies  are  reached.  Here,  the  sixth-order  filter,  due  to  its  finite 
slope  in  the  stopband  region,  improves  slightly  in  its  performance  over  the 
ideal  filter. 
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Fig.  8  Filter  Impulse  Response. 


Filter  Performance  Assuming  an  Output  SNR  of  50 
Ej.,req(10"^^W/cm^)  @  SNR  -  50 


Image 

1st 

Order 

6th 

Order 

+  Thresh 

-  Thresh 

+  Thresh 

-  Thresh 

101 

5.09 

6.51 

9.08 

4.27 

105 

4.84 

15.40 

7.75 

4.26 

106 

7.36 

11.06 

10.45 

8.05 

108 

15.97 

18.32 

19.92 

11.66 

109 

5.26 

6.05 

8.59 

4.81 

Average 

7.70 

11.47 

11.16 

6.61 

Average  of  ideal 

matched  filter  case, 

same  images  5.29 


Average  of  ideal 
Bandpass  filter  case 
same  images 


5.48 


2 


Butterworth  digital  filters  of  first-  and  sixth-order  appear  to  be 
fairly  successful  as  approximations  to  matched,  or  prewhltenlng  filters,  and 
ideal  bandpass  filters,  respectively.  Because  of  their  real  nature  and  the 
non-gaussian  and  non-stationary  properties  of  scenes,  these  filters  do  not 
perform  as  well  as  matched  ideal  filters  in  suppressing  clutter,  nor  do  they 
relate  to  each  other  as  in  the  ideal  case.  The  sixth-order  filter  performance 
(assuming  a  negative  threshold  is  set)  is  better  than  the  first-order  filter 
performance,  whereas  theoretically  the  prewhitening  filter  outperforms  higher- 
order  filters. 


6.0  Summary  and  Conclusions 

A  data  base  of  cloud  Imagery  for  missile  seeker  performance 
evaluation  has  been  established  and  characterized.  Included  in  the 
characterization  is  an  analysis  of  sharpness  of  cloud  edges.  This  work  has 
shown  that  the  boundary  between  clouds  and  terrain  have  very  large  gradient 
values.  There  is  also  a  great  deal  of  edge  structure  within  cloud 
formations.  The  magnitude  and  steepness  of  Intra-cloud  edges  are,  however, 
small  compared  to  the  edges  between  formations  and  equivalent  to  those  within 
terrain  observed  at  long  range. 

Cloud  data  are  used  to  estimate  the  minimum  detectable  target 
Irradiance  of  future,  hlgh-sensltlvlty ,  scanning  missile  seekers.  Both  power 
spectral  density  characterizations  and  directly-filtered  raw  data  have  been 
used  to  determine  the  impact  of  scene  structure  on  detection.  Comparison  of 
these  two  techniques  shot^d: 

o  For  a  given  threshold,  and  filter  type  directly-filtered  scenes 
showed  higher  clutter  compared  to  using  processed  data 

o  Directly-filtered  scenes  required  higher  thresholds  than  that 
expected  assuming  Gaussian  statistics  to  achieve  a  given  false 
alarm  rate.  There  »ras  a  large  scene-to-scene  variation  in  the 
threshold  required  %rlth  directly  filtered  scenes. 
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Matched  filters  give  the  best  clutter  suppression  when  the  scenes  are 
characterized  by  their  power  spectral  density.  Ideal  high  pass  filters 
approach  the  performance  of  a  match  filter,  but  are  never  better.  Using 
directly-filtered  images  the  reverse  is  true:  sharp  cuton  bandpass  filters 
perform  better  than  'whitening'  filters.  This  result  will  be  investigated 
further . 

Finally  correlation  between  steep  gradients  in  scenes  found  by  edge 
detection  and  high  scene  clutter  after  direct  bandpass  filtering  has  been 
made.  This  connection  will  be  further  investigated  to  determine  the  utility 
of  image  characterization  in  making  detection  estimates  in  clutter. 

This  work  will  continue  with  further  investigation  of  filter 
performance  on  an  expanded  data  base.  Additional  discrimination  techniques, 
such  as  time-domain  pulse-shape  analysis,  will  be  explored. 
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over  matched  filters. 


Infrared  Cloud  Backgrounds  and  Sensor  Performance 


Useability  of  our  data; 

Sensor  sensitivity  requirements. 


INFRARED  BACKGROUND  PROCESSING  INVESTIGATION 


THE  NAVY  BACKGROUND  MEASUREMENT  AND 
ANALYSIS  PROGRAM  (BMAP) — FY  1984 

Bernard  V.  Kessler 
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Scannac 

Collecting  Area  (3.9-4.&U) 

(8-11) 

80  cn? 

Optical  Transaission  (3.9-4.87) 

85%  ^ 

Optical  Transaission  (8-11) 

80% 

Detector  Collecting  Efficiency (3.9-4.8U) 

^  (0-11) 

80% 

Total  Optical  Efficiency  (3.9-4.8U} 

68% 

Total  Optical  Efficiency  (8-11) 

64% 

Detector  Area  (3.9-4.8U) 

(8-11) 

•  •  *  *  *• 

2.58  X  10”®  CO? 

Detector  D*jvv  (1  Hz  -  1  kHz  P)  f3.9-4.8W) 

2  X  10^^  cm  Hz 

Detector  D*av  <1  Hz-1  kHz  P)  (7.6-11.3W) 

1/2 

1.5  X  lOlO  cm  Hz/W 

Electronic  Bandwidth  (3.9-4. 8U)  ' 

(7.6-11.3H) 

Hz  -  1  kHz 

ff.S 

Hunber  of  channels  single  color 

16 

2  color  (per  color) 

8 

Elevation  F.Q.V. 

/ 

5.33  mr  or  0.31  deg 

Azimuth  P.O.V. 

49  ar  or  2.81  deg 

Instantaneous  P.O.V. 

0.33  X  0.33  mr 

HEX  3.9  -  4.8y 

1.5  X  10”^^  watts/ca^ 

HEX  7.6  -  11. 6y 

2  X  10“^®  watts/ca^ 

Scanning  frequency 

2  scans  per  second 

Scanning  speed  during  data  acquisition 

36®/sec 

TABLE  1  (Cont.) 


2,4,5,10  or  continuous 


75  in/sec 


NRC 


"NR 


MONTAUK  CLOUD  DATA 
NAV100  FRAME  16 

9  HOURS  43  MINUTES  56.8  SECONDS 
CHANNEL  4 


SAMPLE  NUMBER 


Objective 

-  Assess  Raytheon  Sensor/Data  Recorder  Performance 


nrC 


POWER:  900  WATT  (sTART-Up) 
400  WATT  (continuous) 


r 


Low  PRESSURE  N2  BOTTLE;  UNITHON  POWER  CONVERTER 


nrc 


ANCILLARY  DATA 


PROGRESS  VERSUS  FLIGHT  OBJECTIVES 
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BMAP  Montauk  Measurements  Summary  A.  Hirschhai 
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A  VARIETY  OF  GOOD  DATA  OBTAINED  WITH  EXTENSIVE  ANNOTATION*  GoOD  PROCEDURAL 
EXPERIENCE*  REPORT  IN  PUBLICATION* 


MONTAUK  aOUD  DATA 
MAV100 

HOURS  43  MINUTES  50.3  SECONDS 


SAMPLE  NUMBER 


SAMPLE  NUMBER 


IMPROVED  CALIBRATION  PROCEDURE 


o 

"■  % 

•-x 

s: 

CD 

h- 

< 

■  =5.- 

•— « 

CSS 

ua 

CO 

UJ 

I 

SD 

O- 

t>/a 

o 

CO 

csS 

UI 

CO 

o 

UJ 

CD 

*— • 

• 

oa. 

OS 

sc 

1— 

CD 

zn 

H* 

-c 

CD 

o 

css 

to 

tss 

UJ 

to 

UJ 

o 

-  O- 

■  UJ 

— 

UJ 

u. 

.  -  v  CD 

-J 

CJ 

••  tJ 

U- 

z 

o 

to 

ess 

«: 

UJ 

e 

»>x 

as 

CSS 

o 

SIS 

CD 

to 

Ui 

■sr 

H- 

ixt  UJ 


=3  C3 
DC 

I—  >- 
on  ^ 
2:  3s 
o 


on  on 

=:  H-  o 
•  ui 
ca 


UJ  ►- 
oe:  < 


fc  t— 

PQ  o. 


>-  ^ 


•-x  UJ 
>■  <-> 
•c:  QS 
U  =3 
C3 

-J  c/3 
•< 

ss  • 

cc  u. 

UJ  UJ 
t—  css 


<c 

>•  >- 
UJ  UJ 
-J  _J  to 
UJ  UJ  S£ 


3:  ss:  UJ 
ss  t-> 
UJ  «gc  s 
ca  cu  o 
■a:  o 

s:  >■  css 
n  o_ 
to 

ss  >-  o 


cs)  n  c-> 
s:  szs  CO 
^  <c  . 

CD  CD  I 
s:  SS.  X 


CO  <  o 
C^  UJ 
•  to 


UJ  J—  »— 
H—  CO  O 
UJ  I  S 


c  as 


os  es  Q£ 
<C  <  CD 
UJ  UJ  u_ 
pc;  PS 

to 

f=J  PS  I— 
UJ  UJ  CD 

h-  CO  ea- 

<  t—t 

os  >  cas 
exa  o  -sc 

^  OS  UJ 

_J  u.  ss 
c  s;  •— 


CD  to  to 

■-4  «a: 


<  a  h- 

ESS  »-i 

UJ 
CD  SZ 

as  a  • 
•  ca  UJ 
cs  z 

<C  UJ  UJ 

>—  es  cj 

CO  <  CO 


1 

CD 

u_ 

pa 

1 

CD  >»«  _J 

<C 

css 

UJ 

1 

CO 

o 

CD 

»-x 

to 

1— 

Z 

u. 

OS 

ca 

CD 

»— < 

UJ 

o 

to 

UJ 

OS 

t— 

cs_ 

»— 

<c 

-J 

.D 

UJ 

u* 

SS3 

CSS 

c=s 

to 

u. 

u. 

H* 

is 

V-.v," 

ii 


I 


>■-.- 


!• 

►T'" 


*  *•»* 

t\'i\ 


RMAP  DATA  SHEET 

FIqutauk  Point  Llghx 


nATF/TlME  (LQCAlJ. 

8/15/85  0934 
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<;rFMF  PESCRrPTIOrt 

Broken  clouds  near  solar  azimuth^  varied  illumination 


fsFMsoR  Data.  ..  '  .  '  . 

Boresight  (tripod). setting  az 
Internal  ref  source  temp  10Q°F 


TELESCOPE  AMBIENT  TEMP  74.3°F 
nifiTTAt  Tape  Data  (scanner) 


Stop 


LWIR 

lS0q3440  _ 

+5  SEC 

'8/8 
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+20  sEcr 

1  WIR 

. "'3520 
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ALT 
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+10  SEC 

Analog  Data,  (scanner) 
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Channels  8/8 
Counter  - _ 
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JL20L 


EL  6.S* 


Reel  Ho.  HAV  T07 


Counter 


,  ✓ 

V* 


t  miA  mjfmtCdKik  i*C«i 


r.ODEL  PREDICTS  BACKGROUND  CLOUD  CLUTTER  RADIANCE  AND  TEXTURE 

1-  •  MODEL  INDICATES  CLOUDS  TO  BE  MEASURED;  AN  AID  IH  SITE 

SELECTION  .  ' 

I  '  *  * 

2-  MODEL  DEFINES  SUPPORTING  EO  MET  DATA  EASE  REQUIRED  TO 
RECOVER  INHERENT  CLOUD  STRUCTURE 


3-  MODEL  RESULTS  COMPARED  TO  MEASUREi)  DATA 

-  REFINE  MODEL 

-  IMPROVE  MEASUREMENT  PROCESS 


VALIDATED  MODEL  CAN  GEMLRAIE  sVNlHtllC  CLOUD  CLUTTER  DATA 
BASE 

5-  KOl^a/DATA  AGREEMENT  TELLS  US  VJHEN  VrE  ARE' FINISHED; 
DISCREPANCIES  INDICATE  FUTURE  FORK 


...  .  ■  ■  BMP 

#;  PBnfiliflH  rHAR(/lCTEIilZATION)  ■  .  Ij.AY  BB 


—More  characteri7ation  hetri.cs  will  re  added 


f1D-fll52  735  PRESENTATIONS  AT  THE  TRI-SERVICE  CLOUD  MODELING 
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ANALVSES  ALEXANDRIA  VA  E  BAUER  AUG  84  IDA-H-9-VOL-1 
I  UNCLASSIFIED  IDA/HQ-84-28971  HDA9e3-84-C-0e31  F/G  4/2 


MICROCOPY  RESOLUTION  TEST  CHART 
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4.  SENSOR  TECHNICAL  REQUIRB'lENTS  SUMMARY  -yj 


Table  I 

4.1  Dual  Window  IR  Operation 

DUAL  WINDOW  IR  OPERATION 

.  IFOV  (AZ  X  ELj  MR)  V 
.  TFOV  (aZ  X  el) 

ACQUIRE  VIDEO  DURING  BOTH 
..  DIRECTIONS  V 

SCAN  RATE  ; 

FRAMES/SECCONE  FRAME  IS  ONE 
BACK  FORTH  SCAN) 

SENSOR  TO  BE  COMPATIBLE  WITH 
STABILIZATION  JITTER 

TV  VIDICON  BORESIGHT  MOUNT 

SENSOR  WEIGHT 

FPA  (TDI)  COMPATIBILITY 


. ,« 


3-5^m  and 

8-12y<M 

w* 

0.25  X  0.25 

1 

6°  X  1.7°  (ntfov) 

r  ^0“  X  3,4“  (WTFW) 

V 

iis-  '  ..  . 

'  ^  36?/sec 

3  a  Narrow  TFOV 

---•  !- 

0.3  0  WIDE  TFOV 

1 

<  SO/'HRAD 

•  j 

.  i 

Yes  (surface 

based) 

"  « 

<  60  LBS  (airborne 

- - 

configuration) 

YES 
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Table  I  (contii'iued) 


DETECTORS 


NO.  detectors/waveband 

120 

1/f  noise  shoulder 

MWIR  (HZ) 

<0.5 

lwir  (hz) 

<300  . 

NEAT  (pixel-to-pixel)^  -  ^ 

NO  near  zeros  in  jMTF) 

<0.15“C 

HEI  (W/cm2)  (72°/sec  'scan) 

MWIR 

lwir 

<  2  X  lO'l'* 

<  2  X  10'^ 

detector  &  electronics 

DYNAMIC  RANGE 

>90  DB 

PIXEL  REGISTRATION 

SINGLE  COLOR  ' 

BETTER  T.HANiiO.25,  .  ’ 
IFOV  (scan-to-scan) 

COLOR-TO-COLOR 

SYSTEMATIC  CONSTANT 

OFFSET  KNOWN  TO 

WITHIN  0.25  X  IFOV 

DETECTOR  ARRAY  GEOMETRY 

(suggested) 

LINEAR  ARRAYS^ 

STAGGERED^  CONTIGUOUS 

DC  RESTORATION; ON  OR  OFF  LINE 

YES 

•  »  .  *  * 

SELECTABLE  SPECIBAL  •;  ' 

FILTERS  IN  3-5/4M  BAND 

SWITCH ABLE 

AZIMUTHAL  SHAFT  ENCODER  303 

16  BITS  (0.1  MR 

resolution) 


Table  I.  (continued) 


RADIOMETRIC  OPERATION 


RADIOMETRIC  ACCURACY 
ABSOLUTE 
REPEATABILITY 


<10% 
<  3% 


TEMPERATURE  RANGE  OF  CLOUDS 

TO  BE  MEASURED  220° K  TO  320° K 

TEMPERATURE  RANGE  OF  TARGETS  ..  270° K  TO 

TO  BE  MEASURED  T.B.D.  (PLUMES) 


ELECTRONICS 

-PREAMPS  AND  POSTAMPS  ’  . 

^  DYNAMIC  RANGE  90  DB 

■  VIDEO  PASSBAND  "  '  ^  0,0  TO  7500  HZ 

■"  ■  (3-5 

0.2  HZ  TO  7500  HZ 
(8-12/tM) 


VIDEO  COUPLING 


DC  AND  AC, 


Table  I.  (comtimued)’ 


MODULAR  TELESCOPE  DESIGN 
TELESCOPE 
APERTURE 
COLD  SHIELDING 

photons/mil^-sec  X  10^ 

MWIR 
--  LWIR 


>95%  cold 


SHIELDING  EFFECTIVENESS 


ri 


<  50  X  10* 

T.B.D. 

T.B.D. 


^.5  FOCAL  PLANE  ARRAYS  .  ' 

3-5/tM  FOCAL  PLANE  ARRAY  (tDI); 

TYPICAL  CHARACTERISTICS 

AVERAGE  BACKGROUND  PHOTON  £  50  X  10^ 
(photons/mil^-sec) 

.FEUX  ALLOWED  ON  FOCAL  '  "  "  :V' 

PLANE  (3-5y^M) 


TYPICAL  FPA  (TDI)  WELL 
CAPACITY 


10*  TO  10^ 

ELECTRONS 


FPA  (TDI)  size 

FPA  (TDI)  DYNAMIC  RANGE 
OPERATIONAL  TEMP. 

DWELL  T I me/d ET ECTOR 

I  T;  ^ 

TYPICAL  FPA'  (TDI)dETECTOR 
ELEMENT  SIZE 

SCAN  LINEARITY 
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200  MILS  X  200  MILS 

14  BITS  (34  db) 
t<80°K 

200 /t  SEC 

1.5  MILS  (blur  circle 

SHOULD  BE  EQUAL  OR* 
SLIGHTLY  less) 

<  1/10  ifov/dwell  and 

<  1/10  IFOV  TOTAL  ACCUMU 
LATED  ERROR  OVER  FPA 
(tDI)  focal  PUNE  AREA 


^  U  .  I  I,  l  I.  1  i  I  V  .'  .  ■' 


PROGRESS 


o  Completion  of  dual~band  leased  sensor  by  Raytheon 
-  New  calibration  technique 

0  BHAP  Test^I^  and  Analysis^^)  Pl;\ns 

0  Successful  correction  of  resistive  coupling  defect^^^ 

o  Successful  accomplishment  of  first  surface-based  measurements 
at  Montauk  Point,  L-  I-  NY,  Aug  1983^5^ 

o  Effective  sample  error  correction  software  developed, 
validated  and  documented 

0  Quick  Look  Analysis  (QLA)  Software  development  validation  and 
documentation 

0  NSWC/BMAP  Sensor  Procurement  awarded  Jan  198^1  to  Raytheon 
Company,  $1-45  f1,  delivery  Jan  1986 

o  Production  of  research  tapes  and  dissemination  in  NATO  RSG 
IMAGE  transfer  CCT  FORMAT 


Tasks  to  be  accomplished  yet  next  year: 

0  Second  Surface  Based  Measurements,  late  Summer  1984 

o  Development  and  fabrication  of  HRL/BMAP  240-Channel  Data 
Acquisition  System  (DAS),  Fall  1984 

o  Modification  of  leased  Raytheon  scanner  for  wide  FOV  operation 
(20°  azimuth).  Fall  1984 

o  Absolute  calibration  of  Raytheon  sensor 

0  Development  of  forward-looking  IR  dome  window  for  P3 

o  Development  of  a  passive  line-ofts ight  stabilization  for 
airborne  operation.  Fall  1984 

o  Integration  of  the  dome  window,  stabilization  system  and 
Raytheon  sensor  into  P3  (Winter  1984) 

o  First  airborne  measurement  series,  late  Fall  1984 
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R-  Steinberg,  A-  Hirschman  and  A.  Blumenthal,  'Navy  Infrared  Backgrounds  Measurement 
AND  Analysis  Program— Test  Plan,'  NRL  memorandum  (in  publication).  Naval  Research 


< 


h- 

< 

z 

ca 

o 

1 

z 

-J 

1 

UJ 

oe: 

£ 

QC 

as 

< 

Z 

_l 

• 

u 

QC 

o 

z 

as 

z 

z 

o 

o 

O- 

» 

CA 

t- 

z 

< 

< 

z 

M 

u 

1- 

c/5 

w 

< 

> 

-1 

ca 

Q 

-1 

ca 

Z 

< 

a 

a 

Z 

z 

z 

UJ 

< 

c 

z 

z 

z 

< 

z 

Q 

z 

z 

z 

z 

< 

z 

z 

z 

H 

a 

< 

Z 

Q 

z 

UJ 

Z 

o 

z 

X. 

< 

z 

o 

LU 

z 

z 

z 

ot 

o 

a 

z 

CA 

z 

CA 

UJ 

o 

< 

z 

u 

LU 

< 

H 

_l 

z 

< 

cac 

M 

N 

CA 

as 

I-* 

M 

Q 

z 

Z 

Z 

• 

<c 

a 

% 

ca 

o 

z 

z 

< 

ec 

< 

o 

h* 

19 

CA 

CO 

a. 

z 

u 

z 

t- 

< 

CA 

c>o 

z 

oa 

u 

19 

CA 

z 

Q 

>- 

o 

CO 

UJ 

-1 

1 

o 

a 

< 

z 

z 

ca 

< 

Z 

o 

o 

z 

*1^ 

1 

u. 

»- 

z 

z 

z 

< 

< 

z 

z 

o 

1- 

z 

o 

H 

< 

-J 

iD 

>- 

ca 

z 

1- 

Z 

> 

1 

z 

< 

< 

1 

-1 

> 

u 

z 

z 

LU 

(/I 

« 

< 

« 

z 

CA 

z 

< 

z 

o 

CA 

OQ 

3C 

C9 

% 

< 

a 

C9 

O 

C9 

t- 

O- 

z 

Z 

Z 

z 

< 

« 

>- 

lU 

O- 

z 

o 

z 

ec 

CQ 

CO 

M 

o 

Z 

CA 

z 

z 

M 

OQ 

< 

UJ 

CA 

z 

a 

U 

H 

oe 

>■ 

l- 

z 

a 

Z 

o 

co 

-1 

C/9 

_j 

O 

CQ 

< 

z 

z 

< 

• 

Z 

• 

• 

z 

QC 

<c 

oe 

oa 

QC 

A.  Hirschman,  et  al,  NSWC/RA2,  'BMAP  Surface-Based  Background  Measurement 
Activities  at  Montauk  Pt*,  N*Y-,  Aug  1983' 


TITLE:  AN  EMPIRICAL  CLOUD  MODEL  FOR 
THE  SWIR  &  LWIR  PASSBANDS 


Authors:  TBE,  Mike  Scarborough,  Bob  Pilgrim 

Sponsor:  NV&EOL,  Ground  Systems,  Ted  Doepel 


ABSTRACT 


An  empirical  approach  has  been  used  to  synthesize  imagery  of 
clouds  as  recorded  in  two  infrared  spectral  bands.  Statistical 

parameters  characterizing  cloud  textures  in  two  IR  bands  were  derived 
from  measurements  acquired  by  the  U.S.  Army's  Ground  Based  Measurement 
sensor.  A  simple  two-dimensional  autoregressive  model  was  used  to 
approximate  the  texture  of  cloud  images,  and  a  quasi-fractal  model  was 
used  to  generate  synthetic  cloud  boundaries.  A  procedure  for  generating 
sequences  of  synthetic  stochastic  cloud  images  exhibiting  credible 
temporal  correlation  was  developed  and  implemented. 


CLOUD  SCENE  SYNTHESIS 


1.  IMAGE  MODEL 

The  generation  of  synthetic  imagery  to  simulate  measured  infra¬ 
red  cloud  imagery  requires  the  use  of  a  suitable  mathematical  model. 
Since  cloud  structure  is  highly  variable  and  unpredictable,  it  is  natu¬ 
ral  to  choose  for  this  purpose  some  type  of  stochastic  model.  Although 
clouds  are  inherently  three-dimensional  objects  and  three-dimensional 
stochastic  models  might  be  considered  for  their  representation,  it  seems 
more  efficient  to  attempt  to  represent  the  two-dimensional  cloud  images 
directly  rather  than  to  obtain  them  by  projection  of  three-dimensional 
cloud  structure  onto  a  plane.  This  approach  affords  a  choice  from  among 
a  large  number  of  two-dimensional  stochastic  models  which  have  been 
developed  for  purposes  of  texture  analysis  and  texture  synthesis  in  the 
fields  of  digital  image  processing  and  computer  image  generation. 

The  literature  of  texture  analysis  and  texture  synthesis  was 
surveyed  at  some  length  in  an  effort  to  identify  a  stochastic  model 
which  combines  the  attributes  of  simplicity  of  implementation  and  ade- 
quate  fidelity  of  modeling.  Since  these  requirements  are  obviously 
mutually  antagonistic,  a  compromise  is  to  be  expected  in  selection  of 
the  model. 

Study  of  the  literature  suggested  that  a  simple  two-dimensional 
autoregressive  model  should  satisfy  the  requirements,  particularly  since 
it  has  been  used  with  success  by  others  to  model  cloud  imagery  (Ref.  1). 
For  this  model  (as  for  most  others),  the  image  is  regarded  as  a  two- 
dimensional  rectangular  array  of  discrete  elements  (pixels),  each  of 
which  is  described  by  a  pair  of  indices  specifying  its  location  within 
the  array  and  a  value  corresponding  to  the  image  "brightness  or  gray 
level"  at  the  point  represented.  Using  the  notation  Ui,j  to  represent 
the  gray  level  of  the  pixel  in  the  ith  row  and  jth  column  of  the  array, 
the  general  first-order  autoregressive  model  may  be  written  as 

“i,j  “  ai«i-l,  j*a2Ui,j-l'"«3'*i-l,  j-l+  *  i,  j  (1) 


3 


whi-  r.' 


,ind  33  are  constants,  and  is  a  white  noise  process  such 
that  K  1 1  ,  ^  jt  J  =  B  ,0.  Here  g  is  a  constant,  and  *^01,0 
represents  the  Kronecker  delta  function. 

The  constants  ai,a2  and  33  determine  the  autocorrelation  charac¬ 
teristics  of  the  process  defined  by  ui^j.  In  this  study,  it  is  assumed 
that  the  structure  of  cloud  imagery  is  approximately  isotropic.  Al¬ 
though  the  simple  autoregressive  model  of  Equation  1  is  inherently  ani¬ 
sotropic,  it  can  be  made  to  approximate  the  ipotropic  ideal  by  choosing 
ai  =  32  =  p  with  33  =  -p2.  In  this  case,  3  =  1-  p  Thus,  a  quasi¬ 
isotropic  form  of  the  model  of  Equation  1  is 


“i.j  ■  <■“1-1, ““i.j-i 


(2) 


In  this  model,  and  exhibit  the  same  variance. 

It  can  be  shown  that  the  process  U£^j  generated  by  the  model  of 
Equation  2  has  an  autocovariance  function  of  the  form 

R(k,«.)=E  [(uij-Vi)(ui+k,jH-p)]  = 

where  E  [•]  denotes  the  expected  value  of  the  enclosed  random  variable, 
and  p  =  E[u^^jj,  independent  of  i  and  j.  The  corresponding  spectral 
density  function  is 


S(C.n)  = _ 

[  l-pC-prH-p^Cn]  [  l-p  -p  +pJ^T 

5  n  Cn 

2.  IMAGE  ANALYSIS 


(4) 


To  employ  the  model  of  Equation  2  to  synthesize  imagery  repre¬ 
sentative  of  a  given  cloud  type,  it  is  necessary  to  supply  numerical 
values  for  the  model  parameters  p  and  a  defined  in  the  previous  sec¬ 
tion.  These  values  were  determined  by  empirical  means  using  the 
available  data  recorded  by  the  GBH  IR  telescope.  This  instrument  is 
located  at  the  Army  Optical  Station  on  Roi  Namur  Island  in  the  Quajalein 
Atoll  where  it  was  used  to  observe  infrared  phenomena  associated  with 
the  re-entry  into  the  atmosphere  of  ballistic  missile  warheads. 


Incidental  to  this  mission  it  collected  a  large  amount  of  two-color 
infrared  data  on  clouds  which  happened  to  be  present  in  the  area  at  the 
times  of  re-entry.  Using  this  measured  data  it  is  necessary  only  to 
compute  values  of  the  autocovariance  function  R(k,0)  as  defined  in 
Equation  3  for  a  range  of  values  of  the  lag  parameter  k  and  then  to 
select  values  for  p  and  a  which  give  the  best  fit  of  Equation  3  to  the 
computed  autocovariance  function.  The  value  of  0  2  may  be  read  off 

directly  as  r(0,0),  and  p  may  then  be  identified  with  R(1,0)/  a2.  Al¬ 
ternatively,  the  value  of  may  be  taken  to  be  exp  (-1/k'),  where  k'  is 

the  value  of  the  lag  such  that  R(k',0)/  2  =  exp(-l). 

The  autocovariance  function  R(k,0)  was  computed  on  selected  GBM 
images,  suitably  filtered  to  resemble  data  as  it  would  appear  if  record¬ 
ed  by  the  IRST  sensor.  For  this,  purpose  a  standard  statistical 

subroutine  (FTAUTO)  from  the  IMSL  Library  of  mathematical  applications 
routines  was  employed.  The  autocovariance  of  the  signal  from  each  of 
several  selected  detectors  was  computed  for  each  detector  separately, 

and  then  an  average  over  detectors  was  performed  to  yield  a  composite 

value  for  R(k,0),  k  =  1,2,..., 80.  The  resulting  composite  values  were 

plotted  against  the  lag,  k.  The  plot  was  then  used  to  read  off  values 
of  p  and  a  as  outlined  above. 

3.  IMAGE  SYNTHESIS 

Inspection  of  the  autoregressive  model  of  Equation  2  reveals 
that  samples  of  the  process  may  be  generated  recursively  provided 

that  values  of  j  =  1,2, ...J  and  ui^l,  i  =  1,2,. ..,1  are  available. 

That  is,  given  values  of  uj^^j  for  all  pixels  in  the  first  row  and  the 
first  column  of  the  image  array,  the  values  of  all  remaining  pixels  may 
be  directly  determined  by  successive  application  of  Equation  2.  The 
remaining  problem  is  then  to  provide  values  for  pixels  in  the  first  row 
and  first  column  of  the  array. 

The  approach  taken  in  this  study  is  to  supply  the  required  pixel 
values  from  a  combination  of  measured  data  and  synthesized  values  as 
follows.  A  typical  sample  of  data  from  a  single  GBM  detector  (suitably 
filtered  to  simulate  IRST  data)  consisting  of  a  sequence  of  93  succes¬ 
sive  pixels  is  used  as  the  input  "time  series"  to  SUBROUTINE  FTCMP  of 
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the  IMSL  Library.  This  routine  ^lerfornis  an  analysis  of  a  tim*-  series 
based  c  i  a  stochastic  autoregressive  integrated  moving  average  (AKIMA) 
model.  It  develops  values  of  required  model  parameters  from  the  given 
time  series  and  then  uses  these  values  to  synthesize  a  continuation  of 
the  stochastic  process  for  future  time.  With  the  93-pixel  sample  as 

input,  FTCMP  was  used  to  extrapolate  the  sequence  to  give  a  total  of  512 
samples.  These  512  samples  were  then  used  as  the  values  of 
=  1,2,.. .,512  and  again  as  values  of  u^^l,  i  =  1,2,. ..,512.  A 
straightforward  application  of  Equation  2  using  values  of  P  and  O 
obtained  as  described  in  the  preceding  section  then  permitted  filling  of 
the  remainder  of  the  array. 

The  procedure  described  above  may  be  used  to  synthesize  homo¬ 
geneous  cloud  images,  but  of  greater  interest  in  the  present  application 
are  celestial  scenes  containing  a  mixture  of  cloud  and  clear  sky.  The 
approach  used  in  this  study  to  generate  such  images  is  sometimes 
referred  to  as  the  Background-Foreground  model.  In  this  model,  the 
image  array  is  suitably  partitioned  into  regions  of  two  types.  Although 
the  terms  "background"  and  "foreground"  are  conventionally  used  to  label 
these  regions,  "sky"  and  "cloud"  will  serve  better  here.  Homogeneous 
images  of  sky  and  cloud  are  generated  as  described  above,  and  the  par¬ 
titioned  array  is  used  as  a  guide  in  constructing  a  composite  image 
using  the  following  simple  algorithm.  Pixels  in  the  partitioned  image 
are  examined  sequentially,  and  their  type  ("sky"  or  "cloud")  is  noted. 
If  the  pixel  type  is  "sky,"  the  value  of  the  pixel  at  the  corresponding 
location  in  the  homogeneous  sky  image  is  inserted  into  the  composite 
image  at  the  corresponding  location.  Similarly,  if  the  type  "cloud"  is 
noted,  the  pixel  at  the  corresponding  location  in  the  composite  image  is 
assigned  the  value  of  the  pixel  at  the  corresponding  position  in  the 
homogeneous  cloud  image. 

Use  of  the  Background-Foreground  approach  introduces  two  new 
problems:  How  is  the  partitioning  to  be  accomplished  to  give  a  realis¬ 

tic  (irregular)  boundary  line  separating  "cloud"  and  "sky"  regions?  How 
can  the  abrupt  transition  which  will  occur  at  the  boundary  be  made  more 
realistic  without  additional  processing? 

The  partitioning  problem  is  solved  by  employing  a  stochastic 


F  ij  s  s  e  1 1 


interpolation  scheme  proposed  by  Fourni'-r,  Fussell,  and  Carpenter 
(Ref.  2)  to  approximate  fractal  curves.  It  has  been  shown  (Ref.  3)  that 
cloud  perimeters  are  fractal  curves  of  fractal  dimension  D  =  1.35.  The 
method  of  Reference  2  permits  the  generation  of  a  fractal  curve  of  the 
appropriate  fractal  dimension  between  any  pair  of  points  in  a  plane.  In 
the  present  application,  one  of  these  points  is  chosen  to  lie  on  the 
upper  boundary  of  the  image  plane,  while  the  second  point  is  chosen  to 
lie  on  the  lower  boundary.  The  fractal  curve  connecting  these  two 
points  thus  divides  the  image  roughly  into  left  and  right  regions  which 
may  be  identified  with  "sky"  and  "cloud",  respectively.  A  slight  modi¬ 
fication  of  the  method  of  Reference  2  permits  the  introduction  of  a  con¬ 
trolled  degree  of  convexity  into  the  generated  curve  to  simulate  the 
"fluffy"  appearance  of  clouds.  Without  this  feature,  any  curves  gene¬ 
rated  would  be  statistically  symmetrical,  and  the  assignment  of  "sky" 
and  "cloud"  values  to  left  and  right  regions  would  be  completely  arbi¬ 
trary.  The  convexity  parameter  introduced  is  simply  a  constant  bias 
value  added  to  the  random  variable  used  in  the  stochastic  interpolation 


algorithm  of  Reference  2.  Since  no  empirical  value  for  this  bias  has 
been  reported,  a  value  (0.02)  was  selected  which  appears  visually  to 
yield  credible  cloud  contours. 


Plots  of  infrared  cloud  data  for  scans  which  intersect  cloud 
boundaries  typically  exhibit  a  gradual  decrease  in  cloud  radiance  as  the 
cloud-sky  boundary  is  approached  from  the  "cloud"  side.  In  addition,  as 
the  cloud  becomes  more  tenuous  near  the  boundary,  there  may  be  "holes" 
in  the  cloud  through  which  the  sky  may  be  seen.  To  represent  these 
effects,  the  "cloud"  portion  of  the  image  is  modified  as  follows. 

Let  C(x,y)  and  S(x,y)  represent  the  cloud  and  sky  radiances  at 
point  (x,y)  in  the  two  homogeneous  regions,  respectively.  Consider  a 
single  horizontal  scan  line  (row  of  pixels)  from  the  image  corresponding 
to  a  constant  value  of  y,  and  suppose  that  the  cloud  boundary  intersects 
this  scan  line  at  the  point  (xi,y)  with  "cloud"  corresponding  to  x  >  xi , 
"sky"  to  X  <  x^.  We  define  a  modified  cloud  radiance  C'(x,y)  for  x  >  xj^ 
by  C'(x,y)  =  greater  of  (S(x,y)  and  C(x,y)  -  C(xi,y)  exp  ,-(x-xi)/d) . . 
Here,  d  is  a  characteristic  length,  determined  empirically,  which  pro¬ 
vides  a  measure  of  the  rate  of  decrease  of  cloud  radiance  as  the  cloud- 
sky  boundary  is  approached.  In  the  case  where  the  horizontal  scan  line. 


y.  i  i  t  s  tlu'  boundary  at  two  points,  xj  and  X2i  with  the  interval 

(x_),xt)  r  or  re  s  pond  1  ng  to  "cloud",  the  expression  for  C'(x,y)  with  xj  <  x 


C'(x,y)  -  greater  of  (S(x,y)  and  C( x, y ) -CCxj , y )exp( -(x-xi ) /d) - 
C(x2,y)exp(-(x2~x)/d)). 

figures  1  and  2  are  graphic  representations  of  cloud  scenes 
generated  using  the  method  described  above.  They  represent  clouds 
generated  based  on  the  data  for  LWIR  and  SWIR  passbands,  respectively. 

4.  TEMPORAL  PROPAGATION 

In  addition  to  the  ability  to  synthesize  textures  resembling  the 
spatial  distribution  of  infrared  radiation  observed  in  clouds,  it  is 
also  desirable  to  be  able  to  simulate  the  evolution  of  these  distribu¬ 
tions  in  time.  Given  a  synthetic  radiance  distribution  at  time  t  =  t^ 
represented  by  an  array  of  pixels  a£^j,  such  a  model  would  be  capable  of 
providing  a  credible  distribution  over  the  same  set  of  pixels  at  a 

later  time  t  =  tj,.  Any  model  proposed  for  generating  the  array  b£^j 
should  exhibit  certain  statistical  characteristics.  Thus,  in  general, 
the  array  bi,j»  while  differing  in  a  random  manner  from  ai,jt  should 
nevertheless  be  correlated  with  to  a  degree  dependent  on  the  value 

of  the  time  difference  tj,  -  tg.  In  particular,  when  t  =  tg,  the  model 

must  predict  b£^j  =  ai,j  for  all  i  and  j.  For  tj,  greater  than  tg,  the 
model  should  predict  a  monotonically  decreasing  value  of  the  quantity  E 
[ai,jbi,j]  as  the  difference  tj,  -  tg  increases,  and  in  the  limit  as  tj, 
approaches  infinity,  the  expected  value  E  [®i,jbi,j]  should  approach 
zero.  In  addition,  the  spatial  covariance  properties  of  the  b-array 
should  be  identical  to  the  corresponding  properties  of  the  a-array  for 
all  values  of  tj,  and  t^,  at  least  for  time  differences  tj,  -  t^  on  the 
order  of  a  second.  In  particular,  we  must  have  E  [bi^j]  =  E  [ai,j]  and 

^  [’’i,  j^’i+h,  j  +  k]  -  M  ^i+h,j-*-k]  integral  values  of  i,  j,  h,  and 

k. 

Consider  a  model  for  bi^j  of  the  form 
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where  E  fli,j  ~  0|  ***<1  ^  array  generated  using  the  same 

stochastic  process  as  that  used  in  generating  the  a-array  but  statis¬ 
tically  completely  independent  of  the  a-array.  That  is, 


^  (®i,  j®i+h,  jk]  “  ^  [ci+h,j+k] 

for  any  i,  j,  h,  and  k,  but 

^  f^i, jci+h, j+k}  =  0 

for  all  i,  j,  h,  and  k.  Then  we  have 

['’i, j^i+h, j+k]  ~  (®i, j®i+h, j+k] 

+  (1  p2)  E  jCi+h, j+k]  ~  ^  [^i, j^i+h, j+k] • 

The  parameter  p  is  a  function  of  the  time  t  at  which  the  distri¬ 
bution  is  given  by  bj^^j.  The  array  may  be  thought  of  as  the  radi¬ 

ance  distribution  existing  after  an  infinite  time  has  elapsed  following 
observation  of  the  distribution  aj^^j  at  time  tg.  A  reasonable  model  for 
the  dependence  of  p  on  the  time  difference  t^  -  tg  is 

P  (tb  “  tg)  =  exp  [-(tb  -  tg)/T]. 

where  T  is  the  characteristic  time  associated  with  the  evolution  of  the 
radiance  distribution.  The  value  of  T  could  be  determined  empirically. 
With  this  expression  for  p,  the  bodel  for  becomes 

•»i,j  =  exp  (-(tb  -  tg)/T)  aij 

+  /l  -  exp  (-2(tb  -  tg)/T  cij 

which  satisfies  the  requirements  discussed  above. 


All  in  L  c  r  |>(>  1  ,i  t  ini  was  >',i’iu‘rat  mi  as  sliown  in  Fi^'.uri  i,  .m  h  ti.al  lln 


correlation  cool  I  ic  imit  ])- .H  was  assummi  for  the  int  ermed  i.at  e  ima^'e.  i'oi 
this  proceedure  to  |)rovide  a  legitimate  time  sequence  at  images,  t  tie 
relationship  between  tlie  correlation  coeffient  and  elapsed  time  must  be 
determined.  Current  analysis  suggests  that  this  relationship  is  derivable 
1  rom  statistical  image  measurables. 
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Assessment  of  the  L0WTRAN6  Cirrus  Model 
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John  Hornstein 
Optical  Sciences  Division 
Naval  Research  Laboratory 
Washington,  DC  20375 


L0WTRAN6  provides  a  model  of  transmittance  through 
cirrus  clouds  (but  not  emission  or  reflection  from  these 
clouds),  for  wavelengths  spanning  the  near  ultraviolet 
through  the  thermal  infrared.  The  model  rests  on  three 
assertions,  whose  most  general  forms  are:  (1)  extinction  in 
cirrus  clouds  depends  only  on  condensed  water  content;  (2) 
extinction  varies  negligibly  over  the  cited  wavelength 
range;  (3)  condensed  water  content  depends  only  on  cloud 
thickness.  The  present  paper  asseses  the  validity  of  these 
assertions,  both  in  their  general  forms  and  as  they  are 
implemented  in  LOWTRAN.  The  first  assertion  is  supported  by 
both  theory  and  experiment,  but  only  over  a  restricted 
spectral  range  whose  location  and  width  depend  on  the  sizes 
of  the  condensed  particles.  Because  cirrus  particles  are 
large,  this  valid  spectral  range  starts  at  wavelengths 
somewhat  larger  than  the  upper  wavelength  limit  claimed  for 
the  model.  The  second  assertion  is  not  true  over  the  large 
spectral  range  of  the  model:  several  combined  lidar  and 
infrared  experiments  on  cirrus  found  cirrus  to  be  four  to 
six  times  more  transparent  in  the  visible  than  in  the  far 
infrared.  The  third  assertion  is  contradicted  by  most 
published  measurements  on  optical  variations  within  single 
cirrus  clouds.  Physical  arguments  suggest  that  it  is  correct 
in  a  rough  and  qualitative  way  for  comparisons  between 
average  or  peak  values  between  different  cirrus  systems,  but 
there  is  no  reason  to  hope  for  widely  applicable  numerical 
values  of  the  parameters.  Because  the  asserted  relation  is 
not  applicable  within  a  single  cirrus  cloud  or  system,  the 
current  L0WTRAN6  model  is  not  appropriate  for  modeling 
texture  within  a  cirrus  cloud.  The  LOWTRAN  cirrus 
transmittance  model  can  be  corrected  without  great 
difficulty,  but  this  will  require  revision  of  its  built-in 
numerical  parameters.  In  addition  to  the  reasons  already 
given,  revision  is  needed  because  the  current  values  are 
based  on  flux  transmittances,  whereas  LOWTRAN  needs  radiance 
transmi ttances .  A  simple  way  of  making  this  correction  will 
be  described.  Some  effects  of  the  nonspherlcal  shape  of 
cirrus  particles  can  also  be  estimated  simply. 


The  L0WTRAN6  Cirrus  Model 


F.F.  +VUJ.^  JR. 

M.J.  POST 

R.A.  RICHTER  NOAA 

6.M.  LERFALD 
V.E.  DERR 


-  O.IA  (az)^/^ 


TRANSMITTANCE 


T  =o.W(az)2/^ 


OPTICAL  DEPTH 


=o.w(Az)/ 


EXTINCTION  COEFFICIENT 
PER  UNIT  VOLUME:  KM"^ 


31  550  cm' 


rl 


0  cm"^  1000  cm'^ 


ilOOOO  cm' 


-1 


28.6  10^ 


I  0.25  pL 
0.317  ' 


EBAL  FOf^ 


Extinction  depends  only  on  condensed  water  content^ 

AND  IS  ONLY  WEAKLY  DEPENDENT  ON  WAVELENGTH  OVER  CITED  RANGE. 

Condensed  water  content  depends  only  on  aouD's  thickness. 


IN  LIMRAN6 


Extinction  is  directly  proportional  to 

CONDENSED  WATER  CONTENT. 


Condensed  water  content  is  directly  proportional  to 
aouD's  thickness. 
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Fio.  1.  Schanktic  behavior  of  the  nonsaUaed  eatioctioo  oou 
tectiea  (cfidcncy  for  eatinctieo)  Q,wt  a»  a  functioD  of  the  liae 
paraaetcr  a.  la  the  ic|ioa  x»\,  Q,tt  ia  uauall)'  appionaated 
by  ita  aayiaptotic  value  If  tte  droplet  atae  totributioo 

apaaa  fioin  s^O  to  s*«a«,  Qmfa)  caa  be  appiosaatcd  [for 
t^  purpote  of  cvaloatioa  of  aa  intcfral  ia  Eq.  (4)3  by  a  atrmifht 
hae  (?»•(*) 


Chylek,  Jl.  Atmos.  Sci.  296  (1978) 
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Fro.  2.  Xonnaliud  extinction  cnu  section  Qui  for  mter 
droplets  et  «nve]eo(th  X"O.SS  fim  (index  of  refimction  «ii“ U33) 
nod  its  appioximstion  by  a  stnifht  line  0.n(x)-ex  for 
The  spproziinAtioo  ovetestimstei  the  exact  value  of  Q,n  at  small 
siie  parameten  x,  and  undercstiinates  it  at  larger  x  (still  with 
x$xjr).  These  two  errors  tend  to  cancel  out  in  the  evaluation  of 
the  integral  in  Eq.  (4). 


Fic.  3.  As  in  Fig.  2  except  at  Xi*ll  nm.  Due  to  ika  stfaeg 
absorption  (mwl.lS3^.0968i)  at  this  sravelcngth  fb»— rver 

of  the  Qn>(z}  curve  differs  from  that  at  XwO^SSsan.  Still 
curve  can  be  approxiinated  by  a  straight  line  for 

If  the  considered  size  disttibation  ^ans  np  la  aw  the  cffeci 
underestimating  Qm  in  one  region  and  overcatiasatiag  U 
another  region  cancels  out  in  the  inugral  in  Eq.  (4).  li  only  seBall 
dtopleu  are  present,  canoellation  does  aot  take  place  and  the 
(?«tw  cx  approximation  b  less  aocnmtc. 


Chylek,  Jl.  Atmos.  Sci.  ^  296  (1978) 
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Fic.  3.  Particle  spectra  taken  in  the  dira-nshcar  rccinn  near 
the  hase  of  the  head  of  a  cirrus  uncinus  cluud  (A)  ami  near  the 
top  of  a  cirrostralus  deck  (0). 


A.  Heymsfield,  Jl.  Atmds.  Sci.  2. .  799  (1975) 
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SHOWS  THAT  a  CWC  WHERE  EXPECTED^  AND  NOT  WHERE  NOT  EXPEaED.  FOR 
SMALLER  SIZED  PARTICLES  IN  WATER  CLOUDS^  EXPECT  VALIDITY  FOR  10  u  <  X  <  IQO 
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LIQUID  WATER  CONTENT  (g  m'*) 

no.  I.  VariUM  ef  cxlinciioa  eocCBcicni  »iib  liquid  water  cooicat  is  alaiMphcric  fog  aad  haie 
for  Ml  litc  dinribiilioa  otcasuraacau  Bade  si  diffcrcol  feofraphic  locsin  and  aadrr  a  variety 
af  f  irnmliniral  coaditioM.  la  the  infraiad  tpeeual  fapoR  arouad  A  -  II  «iai  (dl  Ibere  caitu  a 
fmar,  aiar-diitribiiuawMcpcadcai  iclatioa  bctweco  the  vohune  caiioclioa  eocOcaeal 
Md  Ifec  liquid  water  coaltal  W  (g  at  **)  af  the  form  of  Eq.  (I).  Cooiaqueiilly.  Ike  fosulu 
of  aB  maaeuramcatt  aia  doic  «a  a  titaiilit  boe.  Tbc  predicted  iclaiioo  bciweaa  caiiaetioo 
w,  aad  liquid  waiar  aoatcal  W  aBBordiag  to  Eq.  <11  ie  abowa  by  the  itiaigbi  Uae.  Oe  We  ether 
haad  el  A  •  0  JS  «tai  (a),  the  Q,  •  ce  apprountaiioo  i>  eot  laiitScd  aad  ao  naaBbiciiaus  lals- 
tioe  banraca  the  asttaedon  and  ^uid  water  coaicai  cuu.  The  larie  eptcad  of  the  poiatf  ta  the 
graph  abowt  that  the  ailittctaoe  eoefBdeai  it  a  fimetioa  a(  the  sue  diiuibulioa  at  a>«B  at  of  the 
Kqtiid  water  eeeicnt.  As  the  wavekaglh  h  taereated  to  A  •  U  iuo  (b)  aad  A  •  4  am  (c)  tbc  Q,  -  rx 
approaiaiatioo  it  islitficd  for  latgcr  dropicu  and  the  rclatioa  (I)  shown  by  tbc  tuaigbi  liect  it 
lirrnaiini  a  aiore  realistic  appreeiauiioa  for  hazes  aad  fogs. 


Jl.  Atmos.  Sci.  IL  1577  (1379) 


NOTE  ALSO;  HAS  SMALL  SCAHER  ABOUT  A  FUNCTION  OF  CWC.,  EVEN  WHERE 

THEY  ARE  NOT  SIMPLY  PROPORTIONAL, 


WAVELENGTH  DEPENDENCE 

LowtranB  cirrus  model  assumes  wavelength  dependence  to  be  weak. 

Putt  and  Davis  find  cirrus  to  be  much  more  transparent  in  the  visible 

THAN  IN  THE  INFRARED.  (ThIS  MUST  BE  SO,  FOR  LIDAR  TO  WORK  WELL  AS  A 
PROBE  OF  CIRRUS.) 

Reasons  for  expecting  a  difference; 

-  SCATTERING  IS  MORE  IMPORTANT  THAN  ABSORPTION  IN  THE  VISIBLE^ 

VWILE  ABSORPTION  IS  MORE  IMPORTANT  THAN  SCATTERING  IN  THE  IR. 

-  THE  IMPORTANCE  OF  VARIOUS  PARTICLE-SIZE  POPULATIONS  VARIES  WITH 


WAVELENGTH. 


Hall,  Jr.,  et  al  (1983),  in  lowtran6  guide: 


o  o  o 


oooo  o  o 


VARLEY  ET  AL.  MEASURED 
PARTICLE  SHAPES  AND 
SIZES;  OBSERVATIONS 
WERE  NOT  SUITABLE  FOR 
DERIVING  CLOUD  THICKNESS. 


AF6L  LWC  vs 
CLX)UD  THICKNESS 


o 


o 
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CIRRUS  LIQUID  WATER  CONTENT  (g-  m») 


Figure  22.  Cirrus  Cloud  Thickness  vs  Liquid  Water  Content 
(LWC)  From  AFGL  Measurements 


Figure  23.  NOAA  Measurements  of  Cirrus  Transmittance. 
Solar  irradiance  was  measured  using  a  nine-channel  photom¬ 
eter/radiometer  and  cloud  thickness  with  a  pulsed  ruby  lidar 


miSSIVITY 

(10-12  MICRONS  ) 


Fic.  2.  Emiuixiiy  «  and  optical  ihicLneu  tj  of  cirrus  as  a  function 
of  cloud  vertical  ihickncat  k. 


FLUX  EMISSIVm^  BUT  LACK  OF  CORRELATION  APPLIES 
TO  RADIANCE  EMISSIVITY,  TOO. 


NEITHER  AVERAGE  THICKNESS  NOR  REPRESENTATIVE 
MAXIMUM  THICKNESS  YIELD  A  CORRELATION. 


FROM:  C.  M.  R.  PLATT.  Jl.  Atmos.  Sci.  2Q  .  1191  (1973) 


. 


▼  NK  4C  tTJ 

re  1.  Lidar  (ixnc-height  profiler  of  the  cimn  bacKscaticr  coefiic»€Di  (0-6W/im  wavelength)  ic 
icnjih  on  27  November  1970.  Simulunecnn  cUjod  emmiviiy  «  (10-I2ftm)  is  plotted  bckjw. 


Platt.  Quart.  Jl.  Roy.  Meteorol.  Soc.  IQL  119  (1975) 


PlaH/  Jl.  Atmos.  Sci.  1191  (1973) 
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J _ I  tittlll - 1 _ I  I  ■  ■  ««ll _ I _ I  I  I  I  I  III  I _ ■  I  t  I  ml 

10-^  IO-»  I0‘*  10"'  10* 

ICE  WATER  CONTENT  ( g  m"*) 


ICE  WATER  CONTENT  (g  m”*) 


Fic.  7.  Una  (•)  and  (b)  oysUl  lengths  plotted  egaiast  ice  water 

conlenL  Mazimuin  crystal  length  ooiresponds  to  particles  ia  eonoentiatiaas  of 
1  per  mm  liie  claa.  Sampling  temperature  is  indicated. 


Heymsfield.  Jl.  Atmos,  Sci.  5^  367  (1977) 


HOW  LOWTRAN  MIGHT  DO  IT 


User  specifies  base  height  amd  thickness  of  cloud. 


Lowtran  obtains  in-cloud  temperatures  from  atmospheric  profile. 


User  specifies  severity  of  vertical  convection^  probably  by 

SPECIFYING  CIRRUS  TYPE. 

LomRAN  INFERRS  IWC  (aND  OTHER  MICROPHYSICS)  FROM  TEMPERATURE^  VELXITY 
AND  POSITION  IN  CLOUD.  ThIS  DETERMINES  OPTICAL  PROPERTIES. 

Random-number  generator  may  be  used  for  texture. 


Lowtran  integrates  along  path. 


Initially,  warn  user  that  validity  is  assured  only  in  the  far  IR.  until  the 
ROLE  OF  2  -  20  micron  PARTICLES  CAN  BE  ASSESSED. 


Cirrus  vs.  Contrails 


c  w  c 

CCm]N  LENGTHS 
POPULATION  1 
POPULATION  2 

PARTICLE  NIMBERS 
POPULATION  1 
POPULATION  2 

PHASE 


CIRRUS  CONTRAILS 


^0.5  gm/m^  (max)  ^  0.1  gm/m^ 

0.005  TO  0.05 

IS  MORE  REPRESENTATIVE 


1  TO  loop 

W  TO  700^  150  TO  ^50  ^ 


10  cm"^ 

10*^  cm’^ 

ICE  ICE^  WATER  OR  MIXED 


THICKNESS 


^  TO  3  KM 


^  TO  1  KM 


02  Ofc  06  oa  10 

Abiortitcon  Emittonc*  Eqx 


0  02  0  6  0  6  0  8  10  1  2  1  6  16 
Absorplion  tmitlonce  EoaIHA'’’') 


Fic  9.  The  1 1  >im  flux  emiltance  as  a  function  of  beam  "ab-  Fio.  10.  As  in  Fig  9  except  for  a  midlaiitude  summer  atmosphere 

sorption  emiltance  '  for  both  upward  and  downward  dtrections 
and  for  the  three  case  studies  indicated.  The  calculations  were 
performed  for  a  high  cloud  positioned  in  a  model  tropica) 
atmosphere.  The  three  different  sets  of  points  represent  the 
three  different  ice  cylinder  radii,  as  given  in  Table  1.  Dashed 
line  curve  for  «„  “  e,,;  doited  line:  curve  for  the  simple  dif¬ 
fuse  relation  of  E()  (I6i. 


DOTTED  CURVES  FITS  =  2e-^  (Tj) 


Platt  and  Stephens,  Jl.  Atmos.  Sci.  3Z,  231A  (1980) 

.t 
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PROFILES  OF  OPTICAL  EXTINCTION  COEFFICIENTS  CALCULATED 
FROM  DROPLET  SPECTRA  OBSERVED  IN 
MARINE  STRATUS  CLOUD  LAYERS 

ABSTRACT 

Because  marine  stratus  clouds  are  persistent  in  the  eastern  portion  of 
subtropical  high-pressure  systems,  models  of  their  optical  properties  are 
pstticularly  needed  for  multiple  purposes  such  as  atmospheric  radiation 
balance  studies  and  EO  system  assessments. 

Airborne  measurements  of  the  water  droplet  spectra  n(r)  (r  is  radius; 
0.23^im  <  r  <  150pm)  were  made  at  about  14  levels  in  stratus  layers  during  May 
and  August  1981,  130  km  SW  of  San  Diego.  Five  parameters  were  calculated  from 
15  n(r)'s  observed  along  each  level  extending  6.44  km  and  from  the  average 
n(r)  over  the  entire  run.  These  parameters  are  the  number  of  droplets  N,  the 
mean  droplet  radius  r,  the  cross-sectional  area  A,  the  liquid  water  content  w 
and  the  optical  extinction  coefficient  k.  for  wave  lengths  K  of  0.53,  3.75  and 
10 .59|jm. 

The  vertical  profile  of  these  parameters  revealed  features  expected  for 
changes  in  n(r)  experiencing  adiabatic  cooling  in  vertically  rising  air. 
However,  the  vertical  profiles  for  the  May  and  August  data  differed  in  a  sense 
consistent  with  a  differnce  in  the  chemical  make-up  of  the  nuclei  generally 
associated  with  marine  and  continental  air  mass,  respectively.  The  division 
of  the  data  into  "marine"  (May)  and  "continental"  (August)  was  supported  by 
other  data.  The  capability  of  a  droplet  model  (Fitzgerald,  JAM,  1975,  v.  14, 


p.  1044)  to  duplicate  closely  w  and  A  below  the  cloud  base  on  a  selected 
"marine"  and  "continental"  day,  using  NaCl  and  (NH^)2  SO^  nuclei, 
respectively,  provided  strong  support  on  the  air  mass  classification  and  on 
the  adiabatic  nature  (well-mixed  layer)  of  the  stratus  layers. 

The  data  supported  the  conclusion  that  the  saturation  level,  the  cloud 
base,  was  near  the  elevation  where  w=0.02  gm 

Vertical  profiles  of  )c  (3\'s)  were  constructed  as  a  function  of  the 
distance  from  the  defined  cloud  base  z*  in  the  range  -300m  <  z*  <  200m.  An 
approximation  to  the  )c's,  a  power  function  of  w,  was  found  to  be  superior 
compared  with  approximations  involving  linear  functions  of  w  or  A.  However, 
the  power  function  varies  with  air  mass. 

The  horizontal  variability  of  the  )t's  was  greater  near  cloud  base  and 
top,  but  the  )t's  were  highly  correlated.  The  k's  were  approaching  the  same 
value  near  z*=200m.  The  scale  size  of  the  horizontal  variations  was  near  3  )an 
on  two  selected  days.  Both  vertical  and  horizontal  optical  paths  would  be 
influenced  by  variations  at  this  scale  size. 

Saimpling  problems  including  artificial  variability  created  by  small 
sampling  volumes  of  the  spectrometer  and  the  path  lengths  required  for 
reliable  estimates  of  n(r)  are  also  considered. 


PROFILES  OF  OPTICAL  EXTINCTION  COEFFICIENTS 
CALCULATED  FROM  DROPLET  SPECTRA 
OBSERVED  IN  MARINE  STRATUS  CLOUD  LAYERS 

V.  RAY  NOONKESTER  (NOSC) 

OUTLINE  OF  PRESENTATION 

•  OBSERVATIONS 

•  INTERPRETATION;  ABOVE  AND  BELOW  CLOUD  EASE 

•  EXTINCTION  COEFFICIENT  PROFILES  (3  ''D 

•  HORIZONTAL  VARIATIONS 

•  SUMMARY 
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Elevations  of  the  two  minute  horizontal  runs  (6.44  km)  (dots)  in  the  stratus 
layers  130  )an  SW  San  Diego  and  elevations  of  the  cloud  base  z 
(w  =  0.02  g  m“^)  and  top  z,.* 


Droplet  spectra  n(r)  at  selected  levels  and  distances  z*  from  cloud  base 
(w  =  0.02  g  m”^)  on  18  August  1981  when  a  'continental'  air  mass  was  concluded 
to  be  present. 


36 


'CLOUD  BASE',  z. 


<W>  -  0  02 


—  _ 

18  AUG  , 

t 

^  / 

^  29  MAY 

J 

REGIONS  FOR 

« 

t 

t 

/ 

29  MAY 

'  -r 

i 

t 

t 

1 

0,  •  model  (FIT2GERALD, 

T-  'T 

FLIGHT 

t 

JAM  14,  1975) 

MAY  AUG 

~  LEVELS 

<W>,  LIQUID  WATER  CONTENT,  9m  m'^ 


Vertical  profile  of  the  liquid  water  content  on  29  May  and  18  Aug  81. 
w(r  )  is  the  profile  (average  for  29  May  and  18  Aug)  expected  for  a  moist 
adiabatic  cooling  of  a  rising  air  parcel.  The  large  deviation  at  level  5  on 
29  May  is  an  atypical  unexplained  perturbation.  The  closed  and  open  dots  are 
predicted  by  Fitzgerald's  model  if  ^*^4)230^  and  MaCl,  respectively,  are  the 
nuclei  and  if  the  relative  humidity  is  100%  at  the  cloud  base.  The  capability 
of  the  model  to  reproduce  closely  the  profile  w(z*)  strongly  suggests:  (1) 

the  layers  are  well  mixed  (adiabatic);  (2)  the  saturation  level  is  near  the 
defined  cloud  base  and  (3)  the  nuclei  on  29  May  and  18  Aug  are  commensurate 
with  nuclei  expected  in  a  marine  and  a  continental  air  mass,  respectively. 


AEROSOL  MODEL  BELOW  CLOUD  BASE 


BASIC  EQUATIONS  (FITZGERALD,  JAM,  V.  lA,  P.  1055,  1975) 


•  'Wet'  (r)  and  'dry'  (r^  )  aerosol  radii  are  related 

6  fo  =  F  (NUCLEI  TYPE,  RH) 

R  =  aR  ^ 

[e  =  F  (RH) 

•  'Dry'  AEROSOLS  (nuclei)  have  Junge  distribution 

- ( V+1 ) 

N(Rj  )  =  (LOGe)  C  Rj  ;  O.li-M^Rj^  5uM 

•  'Wet'  aerosols  have  distribution 


N(R)  =  (Looe)  c  <»  R  ;  81%  ‘  R,<99.5% 

t 

APPLICATION  ON  18  AUGUST  and  29  MAY 

•  RH  ACCORDING  TO  WELL-MIXED  ADIABATIC  LAYER 

-  PH  =  100%  AT  CLOUD  BASE  (w=0.02  GM  '^ ) 
t  Nuclei  parameter  -<  fof- 

(NHi^)2  SOtj:  IS  August 
NaCl:  29  August 

•  C,  ^  FROM  LCW-LEVEL  N(R) 


Integrate  n(r)  from  r  =  Ium  to  20um  for  w 


AEROSOL  STRUCTURE  ABOVE  CLOUD  BASE 
•CLASSICAL*  IN-CLOUD  AEROSOL  HDDELS 

•  SUPPORTIVE  MODELS:  KEIBURGER  AWO  OCIB, 
6E0PHYS.  MON,  NO.  5.  A6U.  1960; 

FITZGERALD.  JAS.  V.  41.  P.  1044.  1975. 

t  GENERAL  FEATURES  OF  MODa£ 

-  ADIABATIC  COOLING  OF  RISING  AIR  PARCELS 
CARRYING  CCN  THROUGH  SATURATION  LEVEL 

-  ACTIVATED  PARTICLES  GROW  IN  SUPERSATURATED 
AIR  FORMING  A  MODE  IN  N(R)  AND  MODE  RADIUS 
INCREASES  WITH  ELEVATION  WHEN  NUCLEI  LIKE 
NACL  ARE  INVOLVED 

COMPARISON  WITH  DATA  ON  29  MAY  AND  18  AUGUST 

t  DATA  ON  29  MAY  FOLLOWS  CLASSICAL  MODEL 

•  DATA  ON  18  AUGUST  APPROXIMATES  CUSSICAL 
MODEL  EXCEPT  A  DISTINCT  MODE  IS  ABSENT. 

-  MANY  SMALL  AEROSOLS  ON  18  AUGUST. 


In- cloud  profiles  demonstrates  differences  expected  for  marine  (29  May)  and 
continental  (18  Aug)  air  masses* 


EXTINCTION  COEFFICIENTS 


'TRUE'  EXTINCTION  COEFFICIENT 

(x)  =  fo  H)  N(R)DR 

APPROXIMATE  EXTINCTION  COEFFICIENTS 


•  EMPIRICAL  REUTION 

K,  =  M 
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.  DISTANCE  FROM  CLOUD  BASE 


X  *  10.59  fitn 


Profiles  for  the  May  and  August  days  and  for  all  days  combined  at 
X  =  10 .59  mn. 


Wavelength.  X  (pm) 


3.75  10.59 


20  m  z  <  200  ( 

1  1 

^1 

8  0 

9.7 

6.1 

^2 

7.0 

7.0 

- 

^3 

- 

25.2 

5  9 

8.2 

9.9 

6.3 

1  1  1 

-40  m  .i  2  <  10  m 

1  1 

0  30 

032 

P.10 

^2 

0.19 

0.19 

- 

^3 

- 

062 

0.15 

0  32 

0  28 

on 

Optical  depths  for  various  extinction  coefficients  and  regions  of  the  average 

cloud  layer.  The  s 

lubscripts  are  associated  with 

the  following  extinction 

coefficients; 

=  71  / 

Qr^  n(r) 

dr  (  • 

true' ) 

^2 

=  2A 

( small 

ka 

37tC 
2p  X 

■  w  ( large  X) 

w 

1^4 

=  aw^ 

( empirical) 

These  -cs  indicate  that  is 

a  superior  approximation  to  k.|  for  3  Xs. 

:w-v 
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hOP123MAL  VARIABILITY  OF 

RELATIVE  VARIABILITY 

«  DEFINED  TO  EE 

FOR  i  '  0.53,  3.75,  10.59um 

N  =  15  (A2S.M  FOR  EACH  n(a)) 

ARTIFICIAL  (INSTRUHEK'TAL)  VAR  1  ABILITY 

•  cr  INCnE..;SES  when  the  EKOFLET  COULT  H  is 

1 

SKALL  Hi  A  REGION  OF  N(fi)  CONTROLLING  k, 

N  <  50  FOR  ASSP-IOC  0C.23-M  -J^t 

N  15  FOR  OAP-200  (1A.2-M  .?i  155 -H) 

•  ARTIFICIAL  variability  EXISTED  AT  HOST  L£VC:.S 

UP  TO  ABOUT  GOm  BELOW  CLOUD  EASE. 

TFERD.^SCALE  SIZE 

t  TRENDS  ARE  PRESENT  ALONG  A  HORIZONTAL  RUN  ..HE: 
THE  SCALE  SIZE  ■,  ■  TrE  RU'  _E: 'TN. 

-  Cr  IS  MEANINGLESS  ahES  ,  5. 

•  -vS  ESTi'.ATEZ  “C  EE  '■■■■  USI'.o  Tr.E  L.'  LEE 

OF  OEVIATIJNS  OF  FFlI  <  >  ALZ'-Z  A  , 

•  A  ,  =  Zk::  is  APPROf’RlATE  FOR  •IjuERATE 
l.;nvz:t::-:  ^fitzja'fal:.  .a^*-  v.  iz,  ::: 


JOi 


.TANCE  FROM  CLOUD  BASE 


Profiles  of  parameters  describing  the  horizontal  variability  and 
interrelationships  for  the  extinction  coefficients  on  18  Aug.  The  average 
variability  was  about  25%  in  the  cloud  for  29  May  and  18  Aug. 


A  =  0  53  /jw 


OPTICAL  DEPTH 


EQUIVALENT  ERROR  e  m  z^,  m 


OPTICAL  EXTINCTION,  km  ’  k  =  16+0.26  z^j 
RANGE  r+  .  =  r(1  ±0.25) 


Example  illustrating  the  range  in  x  at  a  X.  =  0.53  if  varied  by 
±25%  from  the  true  x  along  a  vertical  path  as  a  function  of  cloud 
depth  (k.|  for  average  of  May  and  Aug  data).  The  indicated  range 
would  be  expected  if  the  increases  and  decreases  in  the  cloud 
occurred  simultaneously  (in  phase)  along  the  vertical.  This  maximum 
range  would  be  repeated  every  3  )an,  the  horizontal  scale  size 
indicated  by  the  data. 


The  equivalent  error  in  estimating  the  cloud  depth  to  give  the  range 
of  xs  for  the  25%  variability. 


g] 


A  =  10.59  /ifn 


RANGE 


CLOUD  BASE  =  110  m 


HORIZONTAL  RANGE,  km 


Fyample  illustrating  the  error  in  transmittance  along  a  path  from  the  surface 
Lv,  -loud  base  as  a  function  of  range  for  \  =  10.59  (im.  ThuS/  if  the  true 
base  is  at  z  =  100  m  at  a  range  of  2  km,  the  transmittance  would  be 
calculated  (for  May  and  Aug  data  combined)  to  be  0.3,  whereas  if  the  cloud 
bas^B  -►•re  incorrectly  estimated  to  be  10  m  above  or  below  the  true  cloud  base, 
the  li  ".raittance  would  be  calculated  to  be  0.46  and  0.1  respectively 


1 

REPRESENTATIVE  MEASUREMENTS 


SAMPLE  LENGTH  L  IN  TURBULENT  LAYER 


WHEN  TURBULENCE  IS  THE  ONLY  SOURCE  OF 
VARIATIOIEALONG  A  RUN  (LUMLEY  AND 
PANOFSKY,  THE  STRUCTURE  OF  ATMOS.  TURE, 
196Ai. 

L  ~  2. 

.  SCALES  WITH  MIXED  LAYER  DEPTH,  Z,^ 

IS  ACCEPTABLE  RELATIVE  ERROR 


L  AVERAGED  36  KM  FOR  --lOTo  ON 
29  MAY  AND  IS  AUGUST 


L  c:  2  KM.  ABOUT  -HO  KM 
AM  EXCESSIVE  AMOUNT  OF  TIME  WOULD  EE 


REQUIRED  TO  USE  L  2:  30-H0  KM  WHEN 
TEMPORAL  CHANGES  WOULD  OCCUR  AND 
■■ESOSCALE  VARIATIONS  MIGHT  EE  INCLUDED. 


SUMMARY 


A  PROTOTYPE  DROPLET-EXTINCTION  MODEL  HAS 
BEEN  DEVELOPED  FOR  MARINE  STRATUS  LAYERS. 
AN  IMPORTANT  EO  ENVIRONMENT 

THE  MAJOR  INPUT  IS  THE  ELEVATION  OF 
SATURATION 

THE  "AIR  MASS"  CAN  BE  INCLUDED  IF 
KNOWN 

EFFECTS  OF  LOCAL.  LIGHT  WINDS  WERE 
NOT  EVIDENT 

HORIZONTAL  VARIABILITY  CAN  BE  APPRECIABLE 
IS  MINIMUM  IN  MID-CLOUD 
NECESSITATES  LONG  SAMPLING  PATH 


A  rOETBAH  prograa  Chat  utlaataa  t1)«  probability  of  elood-fraa  Intarvala  (CFI) 
haa  baaa  daralepad  aa  aa  aid  la  aaaaaalag  tha  vulnarabillty  of  aircraft  to 
spaca-baaad  laaara.  Tha  softvara  Intaraetivaly  loads  tba  analyst  through  tbs 
daflaltlon  of  tha  aircraft  path  (aay«fhara  la  tha  Mortharm  Haal^hara)  and  othar 
parsBStars  of  tha  problaa  and  caleulatas  and  displays  tha  results  ia  real  tljsa. 
This  pamits  tba  analyst  quickly  to  anslysa  tha  sansltlTlty  of  cloud  bloekaga 
affaets  to  such  rarlablas  as  look  angla,  saason  of  yaar»  geographic  region, 
path  length,  altitude  of  target,  and  length  of  cloud-free  interval  needed  for 
alssioa  success*  Incorporating  Hoate-Carlo  tachaiquaa,  the  software  slaulates 
flying  along  the  defined  path  nany  tiaas  while  looking,  on  each  flight,  for  the 
occurrence  of  cloud-free  Intervals  greater  than  or  aqu^  to  a  length  of  Interest 
The  nunber  of  flights  during  which  at  leaat  one  such  occurrence  is  detected 
divided  by  the  total  nunber  of  flights  beconas  tha  (71  probability  for  that 
particular  set  of  input  parameters .  Tha  program's  technical  basis  rests  upon  a 
probability  distribution  nodal  relating  cloud-free-llne-of-slght  to  cloud-free 
interval.  This  model  and  some  other  aspects  of  the  software,  such  as  the 
treatment  of  spatial  correlation  of  cloudiness,  were  derived  from  the  earlier 
work  of  Kaliek  and  Allen  at  the  Stanford  Research  Institute. 
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MONTE  CARLO  ITERATION  SEQUENCE 


Si  DISTANCE  ALONG  PATR,  0  EVENT  SV.'ITCHES 

ALTERNATE  CLOUDY  /  CLEAR  INTERVALS 

DRAW  CLOUD  COVER  FROM  FREQUENCY  DISTRIBUTION 

CORRElJ\TE  CLOUD  COVER  (SPATIALLY) 

CALCULATE  P  C  F  L  0  S 

CALCULATE  MEAN  RUN  LENGTH  FROM  P  C  F  L  0  S 

DRAW  RUN  LENGTH  FROM  EXPONENTIAL  DISTRIBUTION 

ADD  RUN  LENGTH  TO  DISTANCE  COVERED 

SET  EVENT  SWITCH  TO  1  IF  RUN  LENGTH  MEETS  EVENT 
SUCCESS  CRITERIA  ▼ 


-NO - 


END  OF  PATH' 


ADD  SWITCH  SETTINGS  TO  EVENT  COUNTERS 


NO — 


m  OF  ITERATtoNS? 


COMPUTE  EVENT  PROBABILITY 

(^ITERATIONS  IN  WHICH  EVENT  OCCURRED>'(TOTAL  ^ITERATIONS) 
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SOP'iE  RESULTS 


Qsiyvd 
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GROUND  -  TO  -  SPACE 
(conversion  of  time  ro  distance  intervals) 


^c"^  (1.99  X  10^)  H. 


-7 


)(-  =  HORIZONTAL  DISTANCE  INTERVAL  AT  TOP 
OF  CLOUDS  (km) 

T  =  TIME  INTERVAL  NEEDED  FOR  LOOK  OR  SHOOT  (SEC 

A  =  RADIUS  OF  EARTH  (km) 

ALTITUDE  OF  SATELLITE  (km) 

ALTITUDE  OF  TOP  OF  CLOUDS  (km) 


SPACE  -  TO-  AIRCRAFT 
(conversion  of  to  distance  u^tervals) 


V  =  velxity  of  aircraft  (km/sec) 


3D  -  B8W  -  MODEL 


>  13  mm  0  05  1.6 

>  25m  m  0  014  2.2 

>  38mm  0-0035  2.7 


3.0 


IN  3D-  modeling: 


COS^a  +  COS^^  +  COS^y  *  I 

CHOOSE  AT  RANDOM;  COSy,  X, 

WAVE  HEIGHT  =  FRACh-*-  D) 
D=U-COSa  +  v  COS/9  +  w  COSy 
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2S  FOR  S  <1 


A  STOCHASTIC  SNAPSHOT 
(Po*0.5) 

WAVELENGTH  :  256  km 
(No  FRACTALS) 

AREA  :  (50km)^,  SKY- DOME  SIZE 


AREA  :  APPROXIMATELY  SKY- DOME  SIZE 


A  STOCHASTIC  SNAPSHOT 
WITH  FRACTALS 

BASIC  wavelength:  340km 
SMALLER  WAVELENGTH:  5  km 
CORRELATION  :  0-95 

AREA  ;  (20  km)^ 
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VERTICAL  CLOUD  SIMULATION 


r- 

K 

k 


u 


i 

UJ 


<p  K  ^  K)  O 
2)  !-  rj  -  w 

-  —  ‘  II 


CVJ 

ro 

r- 

00 

O 

rj 

r- 

\n 

Q 

llJ 

UJ 


o 

Q 

z 

< 

cc 


Id  ^01 


-  1H9I3H 


411 


<D  N-  ^  fO  lO  O 
If)  —  r»  CM  —  N 

—  *  *  I*  r 


Id  gOI  -1H0I3H 


.  J 


VERTICAL  CLOUD  SIMULATION 
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CORRELATION  FUNCTIONS  FOR  CLOUD  COVER 
•  AND  CEILING 

Ralph  Shapiro 

Systema  and  Applied  Sciences  Corporation 

ABSTRACT 

Short-tarra  (0-12  hour)  prediction  of  the  weather  at  a  point  depends 
priaarily  on  the  innadiate  past  meteorological  history  at  the  point.  This 
dependence  is  sufficiently  pronounced  to  permit  the  development  of  a  highly 
auccesaful  ahort-term,  point  forecast  methodology  that  maJcea  exclusive  use 
of  the  current  values  of  the  meteorological  parameters  at  the  point  in  ques¬ 
tion.  This  procedure,  knotm  by  the  acronym  GEM  (Generalized  Exponential 
Markov) ,  was  devised  for  operational  use  by  Bob  Hiller  at  the  Techniques 
Oevelopamnt  Laboratory  of  the  National  Weather  Service.  Unfortunately,  for 
many  military  applications,  forecasts  or  specifications  are  required  for 
locations  without  observing  stations,  or  where  the  observations  are  otherwise 
unavailable.  GEM  is  not  useful  for  such  applications,  and  it  is  necessary 
to  infer  the  weather  at  the  point  in  question  by  other  means. 

Objective  analysis  procedures,  especially  in  data-sparse  regions, 

Buike  use  of  the  fact  that  weather  is  coherent  in  space  as  well  as  time. 
Knowledge  of  the  weather  at  one  point  implies  some  knowledge  of  the  weather 
at  nearby  locations.  The  rate  at  which  the  information  at  one  point  decays 
with  distance  is  a  function  of  the  typical  scale  size  of  the  parameter  in¬ 
volved  and  is  efficiently  summarized  by  means  of  spatial  correlation  func¬ 
tions.  Such  functions  have  been  obtained  for  a  variety  of  meteorological 
parameters  for  several  widely  separated  geographical  regions— Germany,  Korea, 
the  Middle  Bast,  and  Central  America.  Correlation  functions  for  ceiling  and 
cloud  cover  will  be  illustrated  for  these  regions  and  will  be  compared  with 
correlation  functions  for  other  parameters  such  as  pressure  and  teag>erature , 
possessing  greater  coherence. 

Correlation  functions  are  obtained  with  hourly  data,  except  where 
Indicated  as  three-hourly  data. 
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1.  Tliis  series  of  slides  deals  with  the  joint  probabilities  and  areal  persis¬ 
tence  in  the  mesoscale  in  Central  Europe  (Letter  Report  RR-84-16). 

2.  After  an  introduction  of  the  station  locations  and  the  cloud/visibiity 
classes  used  in  this  study  (slides  1  and  2),  slides  3-6  shov/  a  caiparison  of  the 
period  of  record  of  the  10-year  period  1966-76  with  earlier  periods.  As  expected, 
some  climatic  variations  are  displayed.  The  significance  of  the  finding  is  the 
fact  that  1966-76  is  a  time  frame  in  which  adverse  weather  (visibility-clouds) 
exhibits  a  smaller  occurrence  than  during  other  periods  of  record.  Thus,  con¬ 
clusions  about  the  frequency  of  occurrence  of  adverse  weather  may  underestimate 
the  probability.  In  other  climatological  periods  the  chances  may  be  higher. 

For  more  details,  see  reference  10. a. 

3.  The  next  slides  (7-12)  provide  the  maximum  and  minimum  of  the  single  station 
probabilities  of  adverse  v^eather  for  the  period  1966-76  and  the  diurnal 
variation  in  fall,  winter,  and  November  for  two  mesoscale  areas.  The  fall  and 
winter  seasons  and  the  month  of  November  have  been  selected  because  they  portray 
the  seasons  (months)  with  high  frequency  of  adverse  weather.  More  detaiLs  can 
be  found  in  references  10. a.  and  10. b. 

4.  Slides  13  and  14  disclose  the  diurnal  variation  of  adverse  v/eather  at  single 
stations  for  fall.  Slide  13  displays  the  occurrence  of  fog  (visibility  ;^1  Imi) 
and  slide  14  provides  a  comparison  of  the  probabilities  for  ceiling  <8000  ft 
alone  (second  cohjmns  under  station)  and  the  contribution  by  the  additional 
requirenent  that  the  visibility  >  5  miles  (first  columns  under  stations).  These 
two  slides  show  the  basic  conditions  for  evaluating  the  reductions  of  joint 
probabilities  and  the  area  persistence  for  the  mesoscale  study  (for  more  details 
see  references  10. a.  and  lO.b.). 

5.  Slide  15  exhibits  the  joint  probabilities  of  two  stations  for  fall.  Tne 
stations  are  arranged  in  order  of  distance;  SA-HA  are  only  40  miles  apart  and 
SA-FR  have  the  largest  distance  with  90  miles.  Again,  the  two  coluims  under  one 
station  pair  represent  the  joint  probability  for  the  most  stringent  requirements 
ceiling  <8000  feet  (second  col’jmn)  and  ceiling  <  8000  feet  and/or  visibility 

<  5  miles  (first  column).  These  numbers  reflecting  the  percentage  of  joint 
occurrence  of  adverse  weather  are  important  in  some  applications.  They  do  not 
provide,  however,  information  about  the  areal  persistence  viiich  is  presented  in 
the  next  series  of  slides. 

6.  The  areal  persistence  could  be  studied  by  calculation  of  the  linear 
correlation  coefficient  of  adverse  vioather  between  two  stations.  Besides  the 
difficulty  in  interpretation  of  the  linear  correlation  coefficient  of  cloud 
amounts  because  of  their  "U-distribution, "  other  problems  exist.  Ordinary 
calculations  of  the  ILnear  correlation  coefficient  are  restricted  to  the  total 
sot  of  data  of  two  stations  without  permitting  evaluation  of  the  vertical 
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structure.  Another  problem  is  the  expansion  to  include  more  stations  by  use  of 
the  multiple  correlation  coefficient,  v^ch  is  determined  with  reference  to  one 
key  station.  Again,  it  is  unsuitable  to  show  the  vertical  structure.  Thus  this 
author  decided  to  utilize  the  2X2  contingency  table  as  the  basis  and  calculated 
the  tetrachoric  correlation.  The  technical  details  are  given  in  slide  16. 

Since  our  data  sample  comprises  N  -  1000,  the  5%  significance  threshold  of  the 
correlation  is  r  =  0.055.  Any  value  higher  than  the  threshold  must  be  considered 
different  from  zero,  although  the  magnitude  of  the  correlation  may  not  be  of  great 
practical  value  unless  it  is  at  least  r  >  0.3.  It  should  be  noted  that  r  =  0 
if  in  the  contingency  table  the  field  a  =  R^S  j,  where  Rj  and  Sj  (marginal 
distribution)  are  the  probabilities  of  the  occurrence  of  adverse  weather  (or 
ceiling  alone)  at  station  one  or  two,  respectively. 

7.  Slides  17  and  18  display  areal  correlations  of  two  station  combinations  in 
fall.  This  season  was  selected  for  its  highest  probabilities  of  fog  (see 
references  10. c.,  d.  and  e.),  but  the  other  seasons  disclose  little  change  of 
the  principal  features  of  areal  persistence,  although  the  actual  correlations 
;.iay  differ  slightly.  'The  two  station  correlation  of  areal  persistence  decreases 
with  distance  as  expected,  but  scatter  aroimd  the  general  trend  of  decrease. 

For  fog  (visibility  <  1  km,  slide  17)  the  lowest  areal  persistence  appears  in 
the  afternoon  hours. ~The  maximum  occurs  in  the  morning  hours  only  for  the  larger 
distances,  while  closer  stations  show  a  maximum  areal  persistence  in  the  late 
evening  hours. 

The  vertical  stiructure  is  evident  from  slide  18.  It  can  be  deduced  that  low 
ceiling  and  fog  are  localized  phenomena  with  limited  areal  persistence,  and  the 
joint  probabilities  correspond  to  random  coincidence.  However,  clouds  with 
ceiling  above  8000  feet  can  be  accounted  for  in  a  widespread  area  in  the  iteso- 
scale.  This  implies  also  that  probably  cirrus  clouds  (ceiling  >15,000  feet) 
have  the  higtest  area  persistence.  These  findings  agree  with  general  expectations 
from  )<nowledge  about  characteristics  of  synoptic  phenomena  in  the  mesoscale. 

8.  The  three  station  areal  persistence  for  fall  is  found  in  slides  19  and  20. 

The  joint  probability  for  two  close  stations  (slide  19)  and  largest  distance 
(slide  20)  was  chosen.  Correlation  coefficients  of  all  hoiars  were  averaged. 

In  slide  20  the  maximum  correlation  was  added,  selected  irrespective  of  the  hour 
of  the  day.  As  indicated,  the  lower  the  ceiling  or  visibility  the  smeller  the 
area  persistence.  This  implies  that  the  addition  of  a  third  station  can  reduce 
the  joint  probability  of  occurrence  of  adverse  weather  to  virtually  the  minimum 
frequency  as  expected  for  independence.  Only  a  small  reduction  of  the  adverse 
weather  occurrence  is  gained  if  the  conditions  require  ceiling  >2000  feet  and 
visibility  >  4  miles  unless  the  two  stations  are  a  greater  distance  ^art 
(slide  20). 

9.  Finally,  slides  21-23  provide  information  on  areal  persistence  in  fall  with 
four  station  comibinat ions .  The  correlations  in  slide  21  are  based  on  empirical 
probabilities  for  three  stations  and  adding  a  fourth  station.  This  fourth  station 
is  listed  in  the  heading,  i.e.,  the  stations  not  listed  in  the  column  heading 

are  the  throe  stations  of  basic  data.  Again,  the  correlation  coefficient  has 
been  averaged  for  all  hours.  Although  the  correlations  vary  for  the  individual 
combinations  of  stations,  the  general  nature  of  the  three  station  concept 
reappears.  Low  ceiling  and  visibility  show  only  smell  area  consistency,  v^le 
clouds  with  ceiling  >  4500  feet  provide  the  best  area  cover. 
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Slide  22  exhibits  an  alignment  of  two  stations  which  can  be  considered  as 
an  evaluation  of  the  relationship  between  East -West  or  North-South  of  the  meso- 
scale  area.  In  this  mesoscale  area  the  north-south  persistence  appears  signifi¬ 
cantly  higher  than  the  west-east  combination,  but  only  for  the  conditions  <4500 
feet.  Apparently  the  movement  of  clouds  from  west  to  east  produces  a  more 
homogeneous  field  between  the  north  and  south  part  of  the  area  in  this  particular 
case. 

Finally,  the  diurnal  variation  of  the  areal  persistence  for  the  four  stations 
is  given  in  slide  23  for  conditions  of  low  ceiling  and  visibility  and  high  ceiling 
and  visibility  of  five  miles  requirement.  The  minimum  persistence  appears  in 
the  middle  of  the  day  at  10  to  13,  except  for  the  last  column.  Again,  the  highest 
persistence  is  found  in  the  late  evening  hours  or  midnight. 
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COriPARISON  OF  ADVERSE  WEATHER  PROBABILITIES 

(Annual  Summary) 
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SINGLE  STATION  PROBABILITY  OF  ADVERSE  WEATHER 
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A  DISCUSSION  OF  THE  SYNTHETIC  IR  CLOUD 
SCENES  OF  PHOTON  RESEARCH  ASSOCIATES 


Edward  J.  Stone 
Naval  Research  Laboratory 


A  Discussion  of  the  Synthetic  IR  Cloud  Scenes 
of  Photon  Research  Associates 


Edward  J.  Stone 

Code  6520,  Naval  Research  Laboratory,  Washington  DC  20375 


ABSTRACT 

(U)  The  Photon  Research  Associates  cloud  reference  data  scenes 
consist  of  12  synthetic  images  of  high  and  low  cloud,  in  four 
different  infrared  wavelength  bands  with  different  sun  conditions,  as 
seen  by  a  space  observer.  Several  of  these  are  displayed  and 
examined  in  the  light  of  present  understanding  of  cloud  radiance. 
Specific  questions  considered  include: 

•  Absolute  values  of  apparent  cloud  radiance,  and  distribution  of 
radiance  values 

•  Edge  radiance  profiles,  partially  cloud-filled  pixels,  and  the 
apparent  distribution  of  optical  thickness 

•  Scale  and  sharpness  of  internal  cloud  structure 

Examples  of  measured  infrared  cloud  images  from  the  Daedalus  and 
Landsat  Thematic  Mapper  instruments  are  displayed  and  compared  to  the 
synthetic  scenes. 


INTRODUCTION 

One  of  the  most  obvious  and  pressing  needs  in  development  of 
infrared  air  vehicle  detection  techniques  is  the  acquisition  of  a 
sufficient  data  set  for  performing  the  necessary  simulations.  For 
clouds  in  the  downlooking  case,  a  f ew  sets  of  data  have  been 
acquired.  The  most  important  sensors  are  Landsat-D,  Daedalus,  and 
Hi-Camp.  They  have  provided  a  good  deal  to  work  with,  but  for  the 
purposes  of  development  of  advanced  sensors,  they  often  have 
inadequacies,  such  as  small  dynamic  range,  low  sensitivity,  and  poor 
spatial  resolution. 

In  addition  to  these  problems,  which  can  be  overcome,  there  is 
an  inherent  difficulty:  images  acquired  from  high  altitude  sensors 
fail  to  give  the  three  dimensional  structure  of  the  clouds  viewed. 

To  obtain  a  three-dimensional  structure  one  must  make  assumptions 
about  the  optical  thickness  and  microphysics  of  the  cloud  in  each 
pixel  and  the  general  structure  of  the  cloud.  Without  a  three- 
dimensional  reconstruction  of  the  scene,  it  is  impossible  to  do 
simulations  with  the  scene  at  different  aspect  or  illumination 
geometries  or  to  include  the  effect  of  parallax  due  to  cloud  motion 
and  sensor  motion. 
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This  problem  can  be  overcome  by  beginning  with  a  three- 
dimensional  structure  for  a  scene  and  then  calculating  the  radiance 
values  from  first  principles.  The  resulting  scene  is  synthetic,  but 
it  is  more  amenable  to  simulations.  The  disadvantages,  of  course, 
are  that  the  scene  is  only  "real"  to  the  extent  that  the  synthesized 
forms  are  similar  to  those  in  nature,  that  one's  knowledge  of  optics 
is  complete  and  that  one  is  able  to  perform  the  detailed  radiative 
transfer  calculation  thoroughly. 

Some  of  the  requirements  which  a  synthetic  scene  must  meet  are: 

•  Spatial  resolution  at  least  as  great  as  that  of  the  sensor  whose 
performance  is  to  be  modeled.  A  spatial  (angular)  resolution  the 
same  as  that  of  the  sensor  permits  direct  simulations  against 
background  for  that  static  scene  only,  such  as  a  flythrough  viewed  by 
a  staring  earth-synchronous  satellite.  To  adequately  model  the 
effects  of  spatial  changes  of  various  sorts,  such  as  jitter,  a  moving 
sensor  platform,  or  motion  of  the  background  features,  however,  one 
should  have  a  modeled  scene  at  substantially  higher  resolution  than 
that  of  the  sensor  being  modeled. 

•  Adequate  sensitivity,  that  is,  having  a  high  resolution  in 
radiance.  Precision  itself  is  rarely  a  problem,  since  quantities  are 
usually  calculated  with  much  more  precision  than  they  can  be 
measured.  The  requirement  is  for  the  physical  quantities  in  a  scene 
to  vary  from  pixel  to  pixel  in  a  realistic  way,  to  include  both 
gradual  and  precipitous  variation.  For  example,  an  entire  cloud  or 
terrain  region  cannot  be  characterized  by  a  single  temperature.  (A 
related  problem  in  modeling  is  the  inclusion  of  the  right  amount  and 
kind  of  noise  on  the  data.) 

•  Three-dimensional  structure.  The  simulator  needs  to  be  able  to 
model  the  effects  of  in-scene  motion  and  sensor  motion. 

•  'Realistic'  forms.  The  true  requirement  here  is  that  the  modelled 
scene  with  implanted  target  respond  to  various  signal  processing 
schemes  with  the  same  number  of  false  alarms,  and  true  and  missed 
detections,  as  the  real  scene  would.  In  principle  there  is  no 
requirement  that  the  scene  'look'  realistic.  However,  there  is  no 
practical  basis  for  comparison  other  than  the  natural  forms. 


P.R.A.  SYNTHETIC  SCENES  OF  IR  CLOUDS 

Photon  Research  Associates,  Inc.,  of  La  Jolla,  California  have 
been  commissioned  by  DARPA  to  generate  infrared  images  of  the  Earth 
and  atmosphere  from  space,  which  may  be  used  for  simulations  of 
military  satellite-based  surveillance.  A  large  set  of  scenes,  known 
collectively  as  the  Reference  Scene  Data  Base,  has  been  generated, 
which  includes  scenes  of  mountains,  coastlines,  urban  areas,  tundra, 
and  clouds. 
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These  scenes  are  based  on  a  detailed  pixelized  specification  of 
an  area,  including  orography,  spectral  reflectance  and  emittance, 
thermal  conductivity,  and  so  forth  for  each  pixel  in  the  image.  A 
computer  program  then  generates  an  optical  image  based  on  some  choice 
of  wavelength  band,  sun  position,  time  of  day,  model  atmosphere,  and 
viewing  geometry. 

The  clouds  scenes  shown  here  were  generated  by  starting  with  a 
photographic  image  and  generating  from  it  a  three-dimensional 
structure  for  the  cloud.  This  was  done  by  fitting  geometrical 
ellipsoidal  forms  to  the  photographic  image.  "Texture",  meaning  a 
random  variation  in  altitude  of  the  viewed  surface  on  a  scale  smaller 
than  the  scale  of  the  ellipsoid,  was  then  applied  to  the  cloud  form 
using  a  correlated  Gaussian  method.  The  cloud  forms  were  then  cut 
off  below  a  certain  altitude,  to  simulate  the  condensation  level, 
leaving  thin  places  and  "islands"  in  the  cloud.  A  radiance 
calculation  was  performed  which  included  light  reflected  by  the  cloud 
from  the  sun,  light  transmitted  through  the  cloud,  light  emitted 
thermally  by  the  cloud,  and  the  path  radiance  and  transmittance  of 
the  atmosphere. 

Two  physical  cloud  scenes  were  generated,  one  to  represent  low- 
altitude  cumulus  ("Low  cloud"),  and  one  to  represent  a  cumulonimbus 
( thunderhead)  with  a  large  range  of  cloud  altitudes  in  view  ("High 
Cloud").  The  calculation  was  performed  in  several  different 
wavebands  and  for  different  sun  positions  to  generate  12  scenes  in 
all.  The  images  are  512  x  512  pixels,  and  are  available  from  PRA  on 
computer-compatible  tapes. 


FIGURES 

The  figures  presented  in  this  'hardcopy'  are  substitutes  for 
color  image  transparencies  used  in  the  oral  presentation. 


Figure  1  ( A  and  B) .  High  Cloud  over  Ocean,  Middle  Infrared 

Photon  Research  Associates  scene  of  high  cloud  over  ocean  in  the 
middle  infrared,  shown  as  a  coarsely-sampled  3-D  plot  (A)  and  as  a 
contour  plot  (B).  In  the  3“D  plot,  'up'  is  colder,  'down'  is  warmer, 
so  that  the  radiance  plot  is  similar  to  the  physical  geometry  of  the 
cloud.  This  is  a  6^1x64  representation  of  the  original  512x512  image. 

The  histogram  shows  the  number  of  pixels  with  given  radiant 
intensity,  intensity  increasing  left  to  right  (cold  to  warm).  The 
spike  in  the  histogram  at  the  extreme  right  is  the  radiance  of  the 
ocean  background.  The  gap  in  the  histogram  between  the  ocean  and 
cloud  radiance  indicates  that  no  pixels  in  the  scene  were  considered 
to  be  partially  cloud-filled,  nor  did  any  parts  of  the  cloud  have 
very  low  optical  thickness. 
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The  image  is  intended  to  be  an  example  of  cumulonimbus 
( thunderhead) .  It  displays  a  wide  range  of  radiances,  and  thus 
altitudes,  and  might  represent  the  early  congestus  stage  of  such  a 
storm.  No  capillatus  or  cirrus  anvil  appears  to  be  present,  nor 
associated  cumulus  or  stratus  formations. 


Figure  2.  (A  and  B) .  Low  Clouds  over  Ocean,  Middle  Infrared 

PRA  scene  of  low  cloud  over  ocean,  shown  in  the  middle  IR  with  a 
low  sun,  shown  as  both  3"D  plot  (A)  and  contour  plot  (B).  Again, 

'up'  is  cold,  and  'down'  is  warm.  Both  hard  and  soft  edges  are 
apparent,  and  the  histogram  of  intensity  shows  some  intensities 
approaching  those  of  the  warm  ocean.  These  two  plots  are  64x64 
arrays,  coarsely  sampled  from  the  original  512x512  image. 

F i gur e  3 •  Detail:  High  Cloud  over  Ocean,  Long-wave  Infrared 

Enlargement  of  portion  of  Photon  Research  Associates  scene  of 
high  cloud  over  ocean.  All  of  the  edges  appear  hard,  indicating  no 
partially  filled  pixels  and  no  thicknesses  approaching  zero.  This  is 
a  64x64  section  of  the  original  512x512  digital  image. 


Figure  4.  Detail:  Low  Cloud  over  Ocean,  Middle  Infrared 

Enlargement  of  part  of  the  Photon  Research  Associates  synthetic 
scene  of  low  cloud  over  ocean  in  the  middle  IR.  Most  of  the  edges 
show  gradual  variation  to  the  radiance  level  of  the  ocean.  There  are 
few  structures  smaller  than  about  100  meters.  At  this  altitude, 
smaller  features  are  expected. 


Figure  5.  Detail:  High  Cloud,  Long-wave  Infrared 

Enlargement  of  portion  of  PRA  Scene  of  high  cloud  in  the  middle 
IR.  Note  the  filamentary  structures,  with  dimensions  as  small  as  1 
pixel  across,  which  are  synthetic.  While  these  features  do  not 
appear  to  be  "realistic",  they  may  nonetheless  be  suited  for 
simulations  as  generators  of  false  alarms. 

Figure  6.  Detail:  Daedalus  Image  of  Cloud  over  Sierra  Nevada 

This  portion  of  a  Daedalus  long-wave  (10-12  microns)  image  of 
patchy  cloud  over  the  Sierra  mountains  shows  small  features  with  a 
characteristic  dimension  of  2  or  3  pixels.  Some  of  these  features 
appear  to  be  filament-like  in  the  displayed  image.  This  should  be 
compared  to  the  synthetic  texture  in  the  previous  figure.  The 
projected  pixel  sizes  are  approximately  the  same. 
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Table:  Rough  estimates  of  scale  of  texture  as  a  function  of 

alti tude . 

These  estimates  are  based  on  the  assumption  that  the  condensed 
water  in  a  volume  of  cloud  is  equal  to  the  water  vapor  content  (it  is 
usually  less).  It  assumes  an  average  cloud  droplet  radius  of  10 
microns.  The  actual  droplet  size  distribution  in  a  cloud  depends  on 
many  things,  including  cloud  age,  altitude,  and  concentration  of 
aerosol  condensation  nuclei.  In  general,  the  smaller  the  droplets, 
the  more  opaque  the  cloud. 

It  is  possible  at  any  altitude  to  have  diffuse,  unstructured 
formations  with  extents  of  may  kilometers:  it  is  also  possible  to 
have  short  distances  corresponding  to  significant  optical  depths,  for 
example,  by  viewing  a  thin  feature  edge-on.  These  numbers  therefore 
serve  only  as  a  guide.  But  a  scene  of  cloud  at  a  given  altitude 
would  be  be  likely  to  have  features  in  some  distribution  about  the 
dimensions  indicated. 


■  REFERENCE 

1.  Photon  Research  Associates,  La  Jolla,  CA,  Two-Dimensional  Data 
Assessment  and  Applications  Study.  Final  Report  (U),  Volume  I: 
Reference  Scene  Data  Base  (August  1983),  Report  R-058-83,  SECRE 
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Figure  2A.  PRA  synthetic  scene  of  low  cloud  In  the  middle  IR,  and  histogram  of 
the  radiance  distribution. 


Figure  2B.  Contour  plot  of  the  low  cloud,  mid-IR  scene  of  Figure  2A. 
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Figure  Detail  of  scene  of  low  cloud  in  the  middle  IR:  Cloud  edge. 


Figure  5.  Detail  (64  x  64  pixels)  of  PRA  scene  of  high  cloud  in  the 
middle  IR.  Top  of  the  cloud,  showing  texture.  Dark  is  cold.  Note  the 
apparent  filamentary  structures. 


Figure  6.  Section  (64  x  64  pixels)  of  a  Daedalus  measured  image 
of  cloud  over  mountain  in  the  thermal  IR.  Some  possible  filamentary 
features  appear  (it  helps  to  stand  away  from  the  image). 


CHARACTERISTIC  DIMENSION  OF  OPTICAL  DEPTH  IN  CLOUD 


Extinction  K  =  2  x  Condensed  Water  Content  /  Average  Droplet  Radius 


(Assumed 

average 

droplet  radius 

=10  microns) 

Altitude 

Temp 

Pressure 

H20  SVP 

CWC 

Ext.  k 

Length 

(km) 

(K) 

(mb) 

(mb) 

(g/m3) 

(cm-1) 

(m) 

1 

282 

899 

11.4 

9.16 

18. (-3) 

0.5 

*> 

275 

796 

7.1 

5.7 

11. (-3) 

0.9 

3 

269  ' 

701 

4.5 

3.6 

7.2(-3) 

1.4 

262 

617 

2.65* 

2.1 

4.2(-3) 

2.3 

5 

256 

541 

1.37 

1.1 

2.2(-3) 

4.5 

6 

249 

472 

0.69 

0.55 

l.l(-3) 

9.0 

7 

242 

411 

0.35 

0.28 

5.6(-4) 

18. 

8 

236 

356 

0.18 

.145 

2.9(-4) 

35. 

9 

230 

308 

.091 

.073 

1.5(-4) 

68. 

10 

223 

265 

.039 

.031 

6.2(-5) 

161 

Caveat  '  these  numbers  are  based  on  certain  assumptions.  Very  wide 
deviations  from  these  numbers  occur,  in  both  directions. 

The  condensed  water  content  is  obtained  by  assuming  that  the  amount  of 
water  or  ice  in  the  cloud  volume  is  equal  to  the  amount  of  vapor.  It 
nay  be  much  more  or  less.  While  we  indicate  CWC's  of  less  than  0.1  gram 
per  cubic  meter  at  high  altitude,  some  measurements  have  reported  cirrus 
contents  as  high  as  0.25  g/ra3,  which  in  Cum  would  indicate  shorter 
characteristic  optical  lengths,  and  thus  smaller,  sharper  features  in 
cirrus. 

The  atmosphere  used  is  the  US  standard  atmosphere.  The  temperature 
profile  may  vary  considerably  from  this,  of  course.  And  at  low 
temperatures,  the  saturation  humidity  is  a  very  strong  function  of 
temperature. 

The  droplet  size  is  common  for  moderately  aged  continental  clouds, 
l.liere  clouds  have  formed  very  recently  (consider  the  example  of  the 
clouds  formed  by  dropping  a  piece  of  dry  ice  into  a  glass  of  water) 
droplets  can  be  much  smaller  and  characteristic  optical  lengths  much 
shorter. 
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CLOUD  DATA  BASES  ROM  THE  INTERNATIONAL  SATELLITE  CLOUD 
CLIMATOLOGY  PROJECT  (ISCCP) 


The  establishment  of  an  adequate  global  cloud  climatology  is  the  first 
program  sanctioned  by  the  World  Climate  Research  Program.  The  understanding 
of  the  role  of  clouds  in  climate  and  the  proper  representation  of  cloud- 
radiation  feedback  in  models  is  of  primary  importance  to  climate  prediction. 

The  scientific  objectives  of  the  ISCCP  are  to  1)  produce  a  reduced  re¬ 
solution  visible  and  infrared  calibrated  and  normalized  global  radiance  data 
set  from  which  cloud  parameters  can  be  derived,  2)  foster  basic  research  on 
techniques  for  inferring  the  physical  properties  of  clouds  from  the  data  set 
and  3)  to  promote  research  using  the  ISCCP  data  in  an  improved  understanding 
of  the  earth  radiation  budget  and  inferring  other  climate  parameters. 

Data  specification  for  30  day  averages  calls  for  obtaining  the  fraction 
of  low,  middle,  cirrus  and  deep  convection  clouds  to  better  than  5%  and  for 
the  total  cloud  cover  to  better  than  3%.  Also,  the  heights  are  to  be  deter¬ 
mined  to  better  than  1  km  for  cloud  layers  other  than  low  level  clouds  where 
0.5  km  is  the  precision  level.  Cloud  top  temperatures  are  to  be  determined 
to  better  than  l^K.  Spatial  averages  are  to  be  over  approximately  250  km  by 
250  km  boxes  with  sampling  8  times  a  day  for  a  5  year  period.  At  present 
GOES-Uest,  GOES-East,  METEOSAT,  and  GMS  are  the  participating  meteorological 
satellites  and  the  NOAA  7  polar  orb  iter. 

It  has  been  concluded  that  surface  and  atmospheric  correlative  data  are 
essential  to  the  determination  of  the  cloud  climatology  which  will  be  archived 
along  with  the  reduced  radiance  data  set  and  cloud  parameters  set. 

The  First  ISCCP  Regional  Experiment  (FIRE)  is  the  national  program  in 
response  to  some  of  the  research  goals  of  the  ISCCP.  Its  emphasis  is  on  the 
way  clouds  are  parameterized  in  climate  models  due  to  an  inadequate  knowledge 
of  cloud  formation,  maintenance  and  dissipation.  An  interagency  program  has 
been  set  up  to  deal  with  1)  validating  the  cloudiness  and  radiation  para- 
meterizations  in  climate  models,  2)  characterizing  the  marine  boundary  layer 
clouds,  3)  determining  the  characteristics  that  maintain  cirrus  clouds,  4)  de¬ 
riving  the  statistics  of  the  cloud  physical  radiation  properties  and  5)  assem¬ 
bling  the  data  set  in  a  central  facility  to  permit  independent  study. 

While  initial  activity  has  already  begun  on  the  FIRE,  NOAA  is  hopeful 
in  playing  a  major  role  in  the  time  period  FY86  -  FY90  if  our  initiative  is 
passed  by  Commerce,  0MB  and  the  Congress. 
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IMPORTANT  CLIMATE  PARAMETERS  SUCH  AS  PRECIPITATION. 


DATA  SPECIFICATION  FOR  THE  INTERNATIONAL  SATELLITE 

CLOUD  CLIMATOLOGY 

PARAMEJERS -SPATIAL  AND  TEMPORAL  AVERAGES  AND  VARIANCES  (OR  ANOTHER 
STATISTICAL  MEASURE  OF  THE  SHAPE  OF  THE  TEMPORAL  DISTRIBUTION!  ARE 


TIME  AVERAGING-THE  GLOBAL  CLOUD  CLIMATOLOGY  SHOULD  CONSIST  OF  30-DAY 
AVERAGES  FOR  EACH  OF  THE  8  OBSERVING  TIMES  PER  DAY. 


ISCCP  DATA  FLOW  CONCEPT 


ISCCP  CENTRAL  ARCHIVE 

•  REPOSITORY  FOR  PRIMARY 
DATA  PRODUCTS 


INTERNATIONAL  SATELLITE  CLOUD  CLIMATOLOGY  PROJECT 


o  —  T 

<i>'s 

o H  < 
_l  ^  z  < 
O  ^  LU  O 
O  3  H  Z 

<<  -  < 
to  0;  Q  c/3 
D  T  ?  D 


m  to  CO  ^ 

I-  UJ  LU  00 

UJ  o  O  S  ^ 

^  (J  (J  (J  ± 


to  <  < 
00  ^ 
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UNAI  APPl 


ARCHIVED  DATA 


ll>rrERNATK>NAL  SATELLITE  CLOUD  CUMATOLOQY  PROJECT  (ISCCP) 

DATA  CONDENSATION  (VIS) 


INTERNATIONAL  SATELLITE  CLOUD  CLIMATOLOGY  PROJECT  (ISCCP) 

DATA  CONDENSATION  (IR) 


SURFACE  PROPERTIES 
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INTERNATIONAL  SATELLITE  CLOUD  CLIMATOLOGY 

PROJECT  (ISCCP) 

FIRST  RESULTS  (1  JULY  1983)— VISIBLE  (B2) 


NASA  HQ  EE84.736  (1) 


Why  is  FIRE  needed? 


Elements  of  FIRE 


Program  for 


FIRE  MARINE  STAl  lbUCAL  SURVEY  AND 

STRATUS/STRATOCUMULUS  PROCESS  STUDY 


tuilies  FlKt  requires  obsei  valiuns  front  bolb  Cul  S 


Program  for  Cirrus  Clouds 


of  U. 


INTERNATIONAL  SATELLITE  CLOUD  CLIMATOLOGY  PROJECT 


EL  NINO/SO.  OSCILLATION  REGION 


Data  Collection 


Notes  on  the  Viewgraphs 


2  -  The  establishment  of  a  global  cloud  climatology  was  the  first 

sanctioned  program  under  the  World  Cl Imate  Research  Program. 

3  -  The  reduced  resolution  radiances  are  being  saved  so  that  the  cloud 

climatology  can  be  recomputed  if  a  signi ficantly  better  algorithm 
is  found. 

5  -  INSAT  data  will  not  be  available  from  July  1,  1983  through 

December  31,  1984.  NOAA  7  &  8  polar  orbiter  data  will  attempt  to 
fill  this  gap.  Beginning  January  1,  1985,  limited  lower  resolution 
INSAT  data  will  become  available. 

6  -  NOAA/NESDIS  provides  the  correlative  data  on  surface  and  atmospheric 

properties  required  by  the  GPC  for  the  production  of  the  cloud  data 
sets. 

7  -  Canada  actually  became  the  GOES-E  SPC  in  June,  1984. 

9  &  10  -  Upper  left  panel  is  a  view  of  a  smaller  portion  of  the  region  than 

shown  in  the  lower  panels.  The  region  of  interest  is  off  the  northwest 
coast  of  South  America  plus  some  of  South  America. 

21  -  San  Nicholas  Island  has  been  proposed  as  a  site  for  the  MBL  field 
phase.  The  Department  of  the  Navy  has  been  cooperating  with  us 
on  this,  particularly  NRL. 

25  -  Categories  are  for  example  cirrus  alone,  stratus  alone,  double 
layers,  broken  clouds,  etc. 
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2 


27  -  Field  experiments  are  for  approximately  a  3  week  period  in  fiscal 
years  separated  by  a  1  year  gap  for  the  purpose  of  giving  time  for 
analysis  and  making  changes  for  the  second  field  phase. 

NOTE:  FIRE  is  being  supported  by  the  NSF,  NASA,  NOAA,  DOE,  and  to  a 

limited  extent  by  the  DOD.  More  DOD  support  is  encouraged 
(participation  by  individuals  and  funding). 
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CLOUD  RETRIEVAL  FROM 
NIMBUS-7  SATELUTE  DATA  SETS 


PRIMARY  DATA  BASE 


•  SURFACE  TEMPERATURE 

GLOBAL.  EVERY  THREE  HOURS 

•  II5  MICRON  RADIANCE 

GLOBAL,  TWICE  A  DAY 

f  ULTRAVIOLET  REFLECnVTTY  (.36 -  .38/x) 
GLOBAL,  ONCE  A  DAY 

•  SNOW/ICE  MAPS 

GLOBAL,  DAILY 
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SUPPORTING  DATA  BASE 


•  CUMATOLOGICAL  TEMPERATURE 

•  LOWER  LEVEL  ( 0  TO  2  KM) LAPSE  RATE 

•  UPPER  LEVEL  (ABOVE  2  KM)  LAPSE  RATE 

•  TERRAIN  HEIGHT 

t  LAND/WATER  PERCENTAGES 

COVERAGE:  GLOBAL 
RESCjLUTlON:  160  KM  X  160  KM 
FREQUENCY:  MONTHLY  AVERAGE 
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SURFACE  TEMPERATURE  DATA  BASE 


SOURCt:  AIR  FORCE  3D-NEPriANALySIS 

RESOLUTION:  2S  KM  AT  EQUATOR 

50  KM  AT  POLE 

FREQUENCY:  EVERY  THREE  HOURS 

DATA  PERIOD:  STARTING  APRL  m 

INFORMATION:  SHELTER  TEMPERATURE  OVER  LAND 
SKIN  TEMPERATURE  OVER  HATER 

PROBLEMS:  INCORRECT  TEMPERATURES  OVER  ANTARCHCA 


2^ 

115  MICRON  RADIANCE  DATA  SET 

SOURCE:  TEMPERATURE  HUMIDITY  INFRARED  RADIOMETER  aHIR) 

CHANNEL:  115  MICRON 

67  MICRON  (NOT  USED) 

COVERAGE:  GLOBAL  APPROX.  LOCAL  NOON  AND  MIDNIGHT 

RESOLUTION:  67  KM  x67  KM  AT  NADIR 

12  KM  X  2*4  KM  OFF-NADIR 

DATA  PERIOD:  N0V78-  CONHNUING 

(USED  STARHNG  APRL  197® 
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CLOUD  RETREVAL  FROM  INFRARED  ONLY 


ATTENUAnON  CORRECTION 

POLYNOMINAL  IN  SURFACE  TEMP  «  SATELUTE  ZENITH  ANGLE 
BIAS  ADJUSTMENT 


ESnMATlON  OF  UNCERTAINTY  IN  CORRECTED  TEMPERATURE 
20  K  ERROR  IN  ATTENUATION  CORRN 
TEMPERATURE  INHOMOGENOTY  ESHMATED  FROM  AIR  FORCE 
DATA 


DETERMINATION  OF  CLOUD/NO  CLOUD  THRESHOLD 


ESHMATION  OF  MID  &  HIGH  CLOUD  THRESHOLDS 


^6 


THRESHOLDS  FOR  IMPROVED  HInBUS-7THIK  CLOUD  ANALYSIS  (BLCE) 


6  5a  5b 


Q)  SURFACE  TEMPERATURE  FROM  AIR  FORCE  ANALYSIS 
Q)  ADJUSTMENT  FOR  ATMOSPHERiC  ATTENUATION  (2-8C) 

(5)  BIAS  IN  ATTENUATION  ADJUSTMENT  (20 

Q{)  TNICE  STANDARD  DEVIATION  OF  VARIANCE  DUE  TO  HUMIDITY  AND 
HORIZONTAL 

(5)  ADJUSTnENT  FOR  PARTIALLY  HaED  HELDS  OF  VIEN  (20 
(5a)  cloud/no  cloud  THRESHOLD 

(5b)  cloud/no  cloud  threshold  mhen  cumatological  inversion 

PRESENT 

(B)  LON/MID  CLOUD  THRESHOLD 
(7)  MID/HIGH  CLOUD  THRESHOLD 
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ULTRAVIOLET  REFLECTIVITY  DATA  SET 


SOURCE: 

CHANNELS: 

COVERAGE: 

RESOLUTION: 


TOTAL  OZONE  HAPPING  SPECTROHETER  aOHS) 


4  FOR  TOTAL  OZONE  (31  -  34/1)  -  NOT  USED 
2  FOR  REFLECTIVITY  (36  4. 38/1) 


SUNLIT  PORTION  OF  GLOBE  AT  APPROX  LOCAL  NOON 


50  KH  AT  NADIR 
150  KH  OFF-NADIR 


DATA  PERIOD:  SAHEASTHIR 


r-” 


rT*'j •  .  » «  j  1  *'■»  «  <■  •_ 


I.- 


r: 


2-8 


COnPUTATlON  OF  ULTRAVIOLET  REFLECTIVITY 


GEOMETRICAL  ALBEDO  =  F(eo.  e,4>»Pj,0.R) 
WHERE: 

^0 
e 
4> 

Pt 

/5 

R 


=  SOLAR  ZENITH  ANGLE 


=  SATELLITE  ZENITH  ANGLE 


=  AZIMUTH  ANGLE 


»  TERRAIN  PRESSURE 


=  RAYLQ6H  SCATTERING  COEFF. 


=  REFLECTIVITY 


APPROXIMATION  MADE: 


0  SURFACE  ASSUMED  LAMBERnAN  FOR  ATMOSPHERIC  SCATTERING 
COMPUTATION 
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PROPERTIES  OF  TOMS  DERIVED  REFLECTIVITY 

SIMLARinES  TO  VISIBLE 

CLOUDS/WATER  CONTRAST  AS  HIGH  AS  IN  VISIBLE 

CLOUD/SNOW/ICE  CONTRAST  AS  POOR  AS  IN  VISIBLE 

DISSIMLARITIES  WITH  VISIBLE 

NO  SIGNHCANT  LAND/WATER  DIFFERENCE  IN  REFLECnVITY 
EVEN  THE  DESERTS 

SOLAR  GUNT  EFFECTS  ARE  WEAKER 

NO  SIGNFICANT  ATMOSPHERIC  ABSORPTION 


MUCH  BETTER  ABSOLUTE  CAUBRAnON 


F»gur«*  !! 


2-lS 


COMPARISON  OF  THIR  &  TOMS 
CLOUD  RETRIEVAL  ERRORS 


THIR 

LESS  NOISY 


CLOUDS  WARMER  THAN  THRESHOLD 
UNDETECTABLE 


MAY  OVERESHMATE  COLD  SU&-PIXEL 
CLOUDS(CIRRUS).  UNDERESTIMATE 
WARM  SUB-PIXEL  CLOUDS(CUMULUS) 


MAY  UNDERESHMATE  THIN  CLOUDS 
OVER  LAND  WITH  HIGH  SKIN 
TEMPERATURE  (DESERTS) 


DEPENDENT  ON  THE  ACCURACY  OF  AIR 
FORCE  TEMPERATURE 


STRONG  ATMOSPHERIC  ATTENUAHON 
-  HARD  TO  ESTIMATE  ACCURATELY 


TOMS 

MORE  NOISY 


ALL  CLOUDS  DETECTABLE.  EXCEPT 
THOSE  THAT  ARE  VERY  THIN 


MAY  UNDERESTIMATE  THIN  CLOUDS 
(CIRRUS).  OVERESHMATE  THICK 
CLOUDS(CUMULONIMBUS) 


WLL  DETECT,  BUY  MAY 
UNDERESTIMATE 


NOT  DEPENDENT  ON  EXTERNAL  INFO. 
EXCEPT  FOR  SNOW/ICE 


STRONG  SCATTERING  EFFECTS  -  BUT 
WELL  ESTIMATED 
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COMBINING  THIR  WITH  TOMS 


RECOMMENDED  CLOUD  AMOUNT  =  THIR  * a-W)  + TOMS*  W 

WHERE: 


W  =  Wi  W2 


Wl  =  _!L.  N  =  NO.  OFTOMSSAHRI  lNA9JB-PIXa.AREA 
1  +  N 

W2  =|THIR-TOMSl 
THIR  +  TOMS 


3  GIVES  LESS  W06HT  TO  TOMS  WHEN  THIR  AND  TOMS  ARE  IN  GOOD 
AGREEMENT.  AND  CLOUDS  ARE  PRESENT 


OP  STA 


2-17 


CLOUD  CLASSinCATION 


nVE  TYPES:  GBASED  STRICTLY  ON  RADIATIVE  TEMP) 


•  LOW 

TOP  BELOW  2  KM 


•  MID 

BETWEEN  LOW  &  HIGH 


•  HIGH 

ABOVE  7  KM  AT  EQUATOR 
ABOVE  4  KM  AT  POLES 


•  WARM 

DURING  TEMP  INVERSIONS  ONLY 


•  CIRRUS 

MID  +  HIGH  AMOUNT  WHEN  THIR  SEES  CLOUDS  BUT  TOMS 
HAS  LOW  REFLECTIVITY 

WILL  NOT  BE  IDENTIFIED  IN  PRESENCE  OF  OTHER  MID  OR 
HIGH  CLOUDS 
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ATP  FCPCZ  A-.^TypIf- 

SOURCZ  OF  DAT* 


Cr’A'F'TIO:;'-.!  SHELTER  TH?-'''ZRAT'JFE  RZPORTf  FRPF  ':KC  A*.'D  AIR  WZ-'7:-:Z- 
crpviCE  REOFTVFD  EVERY  'ITEREE  KfRlE.". 

SE~  ER.fRFiCE  TE?-'EERi'Y'JRE  RE^RTF  DERIVED  FROF  KICRCV'VS  FATE!  LITE  A" 
0;,;tp  RE^CRTF  received  E\'EIRY  fix  HC'URF  FRCV  T;AVY  fleet  LAD-'REIr  I  C  a  I 

weather  get? at. 


AIAIYFIF  PROCEDiE^E' 

UFE  FIVE  LAYER  FODEL  A*'ALYF:F  CI  WHOLE  fOESH  GRID  (320  F}-'\ 

I’TTERA ^DATE  to  HLALF  ^CEFE  GRID  (ifC  !<?■''  ATT  EXTRAHYALATE  FROL'  7}“  3?.' DIE” 

leve:  (bc  ME  .AB':""  terba:’:'  tc  ufitg  climatctcgica.i  lapfe  ?a-te: 

A-?LY  DIiJR^'A.l  TEL^E“4TL?.”  CCRRECT'IC’.'T  DERIV'ED  EK-IRI CALLY , 

ITCA'RXR.ATE  MOFT  C'J^-Z'r  FHELTE^  TSy?E?.AIl»REr  ILTC  TKS  AN.AIYSIF  UFTTG 
MODTEIEI  e.ar’:ef  tecrt.'tolt: . 

T’'CCRPr;RATE  YAVY'F  FE‘-  FIHIFACE  TEMPERATURE  ANALYFIS  OT:  K.ALF  MSFH  GRID. 

TA'OORROR.ATE  MOFT  CITLRETT  FEA  SURFACE  TSMPERATUREF  AND  INTERPOLATE  THT 
G'XRAT  FIELD  TO  THE  EIGHTH  MRFH  GRID  ( UO  KM). 

ONCE  AGAIN  APFIY  ALL  CUFRSNT  LAND  FHEITSR  TSMPERATUREF  BUT  TC  THE  EIGHT 
MEFH  FIELD. 


DC  THIF  EVEJ^Y  THREE  HOlTiF  AND  ARCHIVE  AT  NATION  A]  CLIMATE  CENTS? 
ESGINI.TNG  APRIL  %  1979. 


.  AFH\T 


CORRELATION  COEFFICIENT 

0.60 

0.93 

0.91 

0.95 

RMS  DIFFERENCE 

41 

14 

17 

12 

MEAN  OF  ALGORITHM 

(%  CLOUD) 

65 

41 

49 

46 

MEAN  OF  ANALYST 

(%  CLOUD) 

40 

40 

40 

40 

SAMPLE  SIZE 

96 

96 

90 

96 

WATER 

12-3-79 

17Z 

ALGORITHM 

CLE 

BCLE(4) 

TOMS 

NCLE(  6) 

CORRELATION  COEFFICIENT 

0.83 

0.94 

0.91 

0.94 

RMS  DIFFERENCE 

10 

13 

9 

MEAN  OF  ALGORITHM 
(%  CLOUD) 

81 

78 

81 

78 

MEAN  OF  ANALYST 
(%  CLOUD) 

79 

79 

79 

79 

SAMPLE  SIZE 

134 

134 

133 

134 

WATER 

12-4-79 

05Z 

ALGORITHM 

CLE 

BCLE(4 ) 

TOMS 

NCLE( 6  ) 

CORRELATION  COEFFICIENT 

0.86 

0.95 

RMS  DIFFERENCE 

19 

11 

MEAN  OF  ALGORITHM 
(%  CLOUD) 

77 

75 

MEAN  OF  ANALYST 
(%  CLOUD) 

71 

71 

SAMPLE  SIZE 


143 


143 


LAND  6-11-79  17Z 


ALGORITHM 

CLE 

BCLE( 4) 

TOMS 

NCLE( 6) 

CORRELATION  COEFFICIENT  0.93 

0.96 

0.95 

0.96 

RMS  DIFFERENCE 

18 

12 

13 

12 

MEAN  OF  ALGORITHM 
(%  CLOUD) 

51 

37 

42 

39 

MEAN  OF  ANALYST 
(%  CLOUD) 

42 

42 

42 

42 

SAMPLE  SIZE 

68 

68 

68 

68 

LAND  12-3-79 

17Z 

ALGORITHM 

CLE 

BCLE( 4 ) 

TOMS 

NCLE( 6) 

CORRELATION  COEFFICIENT  0.62 

0.67 

0.85 

0.66 

RMS  DIFFERENCE 

38 

39 

22 

39 

MEAN  OF  ALGORITHM 
(%  CLOUD) 

34 

25 

44 

28 

MEAN  OF  ANALYST 
(%  CLOUD) 

51 

51 

53 

51 

SAMPLE  SIZE 

59 

59 

56 

59 

LAND  12-4-79 

05Z 

ALGORITHM 

CLE 

BCLE( 4 ) 

TOMS 

NCLE( 6  ) 

CORRELATION  COEFFICIENT  0.04 

0.97 

RMS  DIFFERENCE 

63 

11 

MEAN  OF  ALGORITHM 
(%  CLOUD) 

89 

50 

MEAN  OF  ANALYST 
(%  CLOUD) 

50 

50 

SAMPLE  SIZE 

38 

38 

CORRELATION  COEFFICIENT 


0.65 


0.87 


0.94 


0.95 


RMS  DIFFERENCE 

34 

22 

14 

12 

MEAN  OF  ALGORITHM 

(%  CLOUD) 

61 

39 

51 

46 

MEAN  OF  ANALYST 

(%  CLOUD) 

48 

48 

48 

48 

SAMPLE  SIZE 

79 

79 

79 

79 

ALGORITHM 

COASTAL 

12-3-79 

CLE 

17Z 

BCLE(4) 

TOMS 

NCLE(6) 

CORRELATION  COEFFICIENT 

0.59 

0.  62 

0.91 

0.75 

RMS  DIFFERENCE 

34 

35 

■H 

26 

MEAN  OF  ALGORITHM 

(%  CLOUD) 

51 

43 

56 

48 

MEAN  OF  ANALYST 

(%  CLOUD) 

58 

58 

58 

58 

SAMPLE  SIZE 

104 

104 

103 

104 

COASTAL  12-4-79  05Z 


ALGORITHM 

CLE 

BCLE(4) 

TOMS 

NCLE(6) 

CORRELATION  COEFFICIENT 

0.71 

0.94 

RMS  DIFFERENCE 

39 

14 

MEAN  OF  ALGORITHM 
(%  CLOUD) 

70 

49 

MEAN  OF  ANALYST 
(%  CLOUD) 

45 

45 

SAMPLE  SIZE 

93 

93 

Zonal  Average  Cloud  Cover 


Zonal  Average  Cloud  Cover 
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Abstract 


A  NEW  DATA  BASE 
OF  CLOUD  VARIABLES 
AT  SUBFREEZING  TEMPERATURES* 


Richard  K.  Jeck 
Naval  Research  Laboratory 
Washington,  DC  20375 


About  9000  miles  of  airborne  measurements  at  various 
altitudes  above  the  freezing  level  in  a  wide  variety  of  cloud 
types  and  weather  conditions  have  been  computerized  to  form  a  new 
data  base  of  cloud  variables  at  subfreezing  temperatures.  About 
half  of  the  data  is  from  the  aircraft-icing  research  flights 
conducted  by  the  National  Advisory  Committee  for  Aeronautics 
(NACA)  in  1946-50.  The  other  half  is  from  recent  flights  by  a 
variety  of  government,  university,  and  corporate  research  aircraft 
equipped  with  modern  cloud  physics  instrumentation.  Most  of  the 
data  is  from  research  flights  over  the  conterminous  United  States 
and  nearby  offshore  areas.  The  data  base  Includes  liquid  water 
content  (LWC),  cloud  droplet  median  volume  diameter,  true  outside 
air  temperature  at  flight  level,  altitude  and  Information  on 
precipitation,  ice  particle  concentrations,  cloud  type  and 
distribution,  weather  conditions,  date,  geographic  location,  and 
other  data. 

Currently  the  data  base  emphasizes  wintertime  clouds  below 
10,000  ft  AGL,  but  data  are  now  being  added  for  high  altitude 
measurements  in  winter  storms  and  In  summertime  convective  clouds. 
A  small  amount  of  foreign  data  (from  northern  Spain)  is  already 
included  and  other  foreign  measurements  will  be  added  In  the 
future.  Data  on  size  distributions,  shapes,  and  other 
characteristics  of  snow  and  ice  particles  are  planned  for  later 


*Research  Sponsored  by  the  Federal  Aviation  Administration 


inclusion  too.  Measurements  from  warm  clouds  can  be  added  easily 
if  there  is  a  requirement. 

The  data  base  is  particularly  useful  for  applications  such 
a  s : 


a)  Verifying  the  USAF  3r)-NEPH  procedure  and  the  Smi th-Feddes 
LWC  and  droplet  size  add-on  procedure,  by  providing  air  truth  data 
for  a  number  of  individual  cases. 

In  addition,  the  basic  assumptions  underlying  the 
Smith-Feddes  model,  such  as  the  dependence  of  LWC  and  dropsize  on 
temperature  and  altitude  above  cloud  base,  can  be  evaluated 
statistically  for  accuracy  and  reliability. 

b)  Characterizing  the  optical  thipkness,  LWC,  droplet  size, 
cloud  composition,  or  other  variables  as  a  function  of  cloud  type, 
airmass,  weather  condition,  altitude,  temperature,  and  geographic 
location. 

c)  Performing  frequency  of  occurrence,  probability  of 
exceedance,  and  multivariate  analyses  on  cloud  variables, 

d)  Determining  extreme  values  of  cloud  variables  and 
identifying  the  weather  conditions,  temperatures,  altitudes, 
geographic  locations,  etc.,  where  they  may  be  encountered. 

Actual  examples  of  the  above  types  of  applications  will  be 
given. 


Security  class.:  Unclassified 


TABLE  A-1.  DATA  CODING  SCliEME 


First  Card  of  each  Data  Suite 


Ca  rd 
Columns 


Mission  Identifier  (Flight  No.  or  Cloud  No.) 

Date  of  Measurement 
Geographic  Location 
Source  of  Data: 

Agency 

Reference  (Publication  ID  or  Report  No.) 
Altitude  Reporting  Convention: 

Scale  used  for  data  in  this  suite 
Elevation  of  local  surface  (hundreds  of  feet, 

ASL) 

Cloud  Information: 

Cloud  type 
Uniformity 

Cloud  base  altitude  (hundreds  of  feet) 

Cloud  base  temperature  (^C) 

Cloud  top  altitude  (hundreds  of  feet) 

Cloud  top  temperature  (^C) 

Weather  Factors: 

Airraass  type 
Weather  description 
Card  1  indicator 


XXXXXXX 


Time  of  Day  for  Event 
Icing  Event  Information: 
Event  No 

Defining  criterion 
Duration  (min.) 

Distance  in  cloud  (nmi) 
Aircraft  State  during  Event: 
Sampling  maneuver 
Airspeed  (kt) 

Altitude  (ft) 

Outside  Air  temperature  (OC) 
LWC  Meter  Data: 

LWC  (g/m^) 

Probe  ID 


MMDDYY 

8-13 

GG  or  GGG 

14-16 

LLL  or  LLLL 

17-20 

ZZ - ZZ 

21-36 

A 

37 

FF  or  M 

38-39 

CC  or  CCCc 

40-43 

G 

44 

BBS 

45-47 

+WW.W 

48-52 

HHH 

53-55 

-YY.Y 

56-60 

aAa  or  Ma 

61-63 

CC - CC 

64-79 

1 

80 

h  Data  Suite 

Card 

Code 

Columns 

hhmnrFh 

1-6 

ee 

7-8 

f 

9 

mm  or  .m 

10-11 

dd  or  .d 

12-13 

s 

14 

wv 

15-17 

aaaaa 

18-22 

+t  t.t 

23-27 

w*w  or  .ww 

28-30 

pp 

31-32 

TABLE  A-1.  DATA  CODING  SCHEME 

(Continued) 

Ca  rd 

Item 

Code 

Columns 

Droplet  Probe  Data: 

LWC  (g/m3) 

w.w  or  .ww 

33-35 

MVD  (ym) 

uu 

36-37 

Max  droplet  diameter  (/iin) 

zz 

38-39 

N  (no/cc) 

nnn 

40-42 

Probe  ID 

PP 

43-44 

Precipitation  or  Other  Large  Particle  Probe  Data: 

PMS  2D-C  Probe  concentration 

n.n  or  .nn 

45-47 

(Particles  per  liter) 

n.n  or  .nn 

4  8-50 

PMS  2D-P  Probe  . 

nn  or  .n 

51-52 

PMS  ID  (OAP)  Probe  . 

nn  or  .n 

53-54 

Other  Ice  Particle  Counter 

P 

55 

Probe  ID 

Ice  Meter  Data: 

LWC  (g/m3) 

Icing  rate  (g  cm“2  hour“l  or 

w.w  or  .ww 

56-58 

cm/hour) 

i.ir  or  .iir 

59-62 

Probe  ID 

PP 

63-64 

Ice  Meter  Data  (2nd  probe,  if  used) 

- 

LWC  (g/m3) 

Icing  rate  (g  cm“2  hour" 3-  or 

w.w  or  .ww 

65-67 

cm/ hour) 

i.ir  or  .iir 

68-71 

Probe  ID 

PP 

72-73 

Weather  Factors: 

Precip.  during  measurements 

(type  &  intensity) 

cfP  or  ,qq 

74-75 

State  of  cloud  particles 

ssss 

76-79 

'  1  rd  2  i  nd  icator 

2 

80 

TAhiL  A-A.  rut  NEW  iUHEK'CUULLU  CLOULi  L/ATA  bAsi  (Lontinuec1> 

Da  ta  E  1  1  e  No  . 

.  - - Lard  Column  No. - 

11111111 1  l<ld22i:.12^223333doJJ3oAA44‘(44444l>d5Sl)51>'6'o5bbfe6&^66667  7  7777  77  7'?0 
i:  ,Mbt>  7  Ago  1 2345b  7bS'Jl  2345b  78901 234ib7bS012345b7by0 1234^6  7890 123436  78901 23456  7890 

^  c  . 

No.  Card  CRecordI  Contents 

1  Cloudo  C13081DDC  UWY  Priv  Comm  55781S24StScC  25  -5.6  45  -8.4cP  Hp,Srl45  1 

?  lAlb4213P  N  NP158  4535  -8. 4 . 18JW. 21 12  m21b  13.4  .253  HN  N  »9  cRO  X  H  CS  U  U  2 

3  19171814P  H  NP163  4030  -8 . 1 . 31 JW. 2812  m283  ill  0  99  NN  N  »9  cRO  X  N  CS  U  82 

4  19:17  15E  3  7L151  3740  -  7 . 4 . 2 0 JW. 24 12  m258  F9.  1  0  m  NN  N  8.2cR0  X  H  CS  U  U  2 

5  19204817P  N  HP157  3190  -b. 8. 22JW. 1910  m405  IS.b  .186  NH  N  7.5cR0  X  N  CS  8  82 

6  19:21  18E.8  2L147  2620  -5.7.05JW.07  7  m322  Fll  . 2  n*  NN  N  «3  cRO  X  N  CS  8  82 

7  19225419P  N  HP138  3190  -6. 8. 1 7 JW. 2010  m343  F14  .198  NN  N  »7  cRO  X  N  CS  8  82 

8  19231820P  H  NP142  3620  -  7 . 3. 21 JW. 2612  m321  F12  0  99  NN  N  •?  cRO  X  N  CS  8  8  2 

9  19235421P  H  NP150  4115  -7 . 5 . 12JW. 2215  ml38  F12  0  99  NN  N  ■lOcRO  X  H  CS  8  82 

10  19242422P  N  NP154  4365  - 7 . 8 . 18JW. 1812  ml7b  F7.4  0  68  NN  N  8  cRQ  X  N  CS  8  8  2 

11  Cloud4  013081wKS  8Wy  Priv  Comm  5/81  S24StScB>26«-5 . 7»47»-8. 5cP  Np,Sri47  t 

12  19.30  E1B.5  1L169  4690  -7 . 4 . 1 7 JW. 1512  ml42  F  . 2  . 2  0  NN  N  •&  cRO  X  H  CS  8  82 

13  19.32  E28  1  3L171  4650  -8. 3. 20 JW. 1612  ml67  F4.0  0  19  HM  N  ml  cRO  X  N  CS  8  82 

14  19.34  E30  2  4L167  4645  -  7 . 5 . 12 JW. 1 1 1 1  ml47  F  .7  0  n*  NN  N  4.5cRQ  X  N  CS  8  8  2 

15  19.36  E4J  411L160  4640  -8. 4. 20 JW. 1812  ml76  Fll  0  m  NN  N  8.5cR0  X  H  CS  8  8  2 

16  19.41  E5J  513L15S  4650  -8. 5 . 22 JW. 1111  mib9  Fll  0  99  NN  N  S.OcRO  X  N  CS  8  82 

17  19.48  E6A1123L1S9  4645  -8.8. 19JW. 0910  ml78  Fll  0  n*  NN  N  5.ScRO  X  N  CS  8  82 

18  Clouds  013081I4CX  8WV  Prlv  Coimi  5/81  524  St  C  29  -8.7  53-10. 5cP  Hp.Sri53  1 

19  202048F1P  N  NP185  5315  -9. 1 . 23JW. 1915  ml03  F3.3  .228  NN  N  »7  cRO  X  N  CS  8  82 

20  202124F2P  N  NP178  4800- 11 . 2. 1 7 JW. 1214  m  88  F14  .199  NN  N  *6  cRO  X  N  CS  8  82 

21  202154F3P  N  NP175  4300  -9. 2 . 16JW. 0912  ml09  F15  .199  NN  N  »6  cRO  X  N  CS  8  U  2 

22  20.22  F4E  2  4L156  3850  -9.7.11JW.IU  9  m250  Fll  .199  NH  N  5.5cR0  X  N  CS  8  8  Z 

23  202506f5P  K  NP163  3325  -9.2.Q5JW.03  7  ml37  F9.2  .179  NN  N  <1  cRO  X  N  CS  0  82 

24  202S42F6P  N  NP161  3005  -8.7.03JW.02  6  m230  F3.5  .222  NN  H  <1  cRO  X  H  CS  8  82 

25  203006G1P  N  NP151  3110  -9.1.04JW.02  7  ml55  Fll  .593  HN  N  «6  cRO  X  H  CS  8  82 

26  203024G2P  N  NP147  3610-10 . 1 . 10 JW. 10  9  m320  H3  .299  NH  H  »6  cRO  X  H  CS  8  82 

27  203100G3P  N  NP155  4125  -9. 4 . 05JW. 0912  *107  F9.1  .168  HN  N  »6  cRO  X  H  CS  8  82 

28  20.32  G4A1429L157  4810  -9. 9. 16JW. 1915  ml08  Fl2  0  99  NN  H  8.5cR0  X  H  CS  8  82 

29  Cloudl  013181L8F  8WY  Prjv  Com*  5/81  S28  Sc  8  51  -8.3  64-11. 2cP  8,C1L0  1 

30  18:13  A6E  513L166  6200-10 . 9.  14  jW.  Ibl2  *165  t-2.0  . 3  *  NN  N  7.2cRO  m  N  CS  8ffiW*12 

31  181848A7P  H  NP179  5625  -9. 5 . 1 1 JW. 091 1  *121  E5.3  .738  NN  N  ■/  cRO.lO  N  CSS-mW*l2 

32  182006A6P  N  HP174  5110  -8.3.09JW.04  9  m  95  Fb.2  .853  NN  N  8  cRO.05  N  CS5-mW»l2 

33  182824B1P  N  NP144  5205  -8.6  0  JW,  02  6  ml  7b  E2.3  . 6  4  NN  N  0  cRO  0  N  CS  8mW«-!2 

34  18284882P  N  NP151  5765  -9.9.04JW.06  7  *314  Fl.O  . 2  9  NN  N  ■ElcRO.05  H  CS  8mW  2 

35  182906B3P  N  HP148  6235- 10 . 8. 22 JW. 211 1  *305  FI. 3  .  3  6  NN  N  8  cRO.22  N  CS  8mW  2 

36  Cloud2  013181LBF  UWYPnv  Com*  5781  S28  51  V  75-11.0  90  -9. 2cP  8,ClL0Wkl75  1 

37  183000C1P  H  NP152  7500-11 . 0 . 01 JW. 03  7  m  76  .747  HN  N  <lcR0.02  N  CSS-  W»12 

38  183036C2P  N  NP159  8125- 10 . 2 . 02  JW.  02  9  m  56  F2.1  .819  NN  N  <.lcR0.03  H  CSS-  W*12 

39  183054C3P  N  NP151  8540  -9.6.04JW.05  9  mll9  P.i.61.236  HN  N  <lcR0.04  N  CSS-  W*I2 

40  183118C4P  H  NP156  8900  -9.4.04JW.05  9  *165  FI. 8  .717  NN  N  »<.lcR0.05  H  CSS-mM*l2 

41  Cloud3  013181LBF  8WV  Ppiv  Com*  5/81  528  St  VllO  -9. 7129-13. 3cP  8.ThCJL0  1 

42  183312D1P  N  HP1651105O  -9.7.02JW.03  5  *167  FI. 8  .8  9  NH  N  <lcRO.02  N  CSS-m«  2 

43  18333602P  H  NP1661 1540- 10 . 8. 10 JW. 10  9  *278  F2.4  .7  8  NH  N  8  cRO. 11  N  CSS-*M«-t2 

44  183406D3P  N  NPl 701 2035-11 . 6. 10 JW. 1 1  9  *i:6u  F2.7  .817  NN  N  *6  cR0.12  N  CSS-mH^IZ 

45  183430D4P  N  NP16312560-12. 6. 13JM. 1410  *254  Fif.O  .926  NH  N  8  cR0.16  N  CSS-mW^ia 

46  183454D5P  N  NPl 70 12940- 13. 3. 04JM. 05  9  *  91  FI. 4  .313  NH  N  ><lcR0.03  H  CS  8  M» 1 2 

47  Cloud4  013181LBF  UWY  Prlv  Com*  5/81  S28  St  1 144-16. 7150-1 7. 8cP  8,noi>100  1 

48  183624E1P  N  NPl 7314430-1 7 . 3. 01 JW. 01  6  *152  F4.01.415  NN  N  «0  cRO  0  N  CSS-W-JISZ 

49  183654E2P  N  NPl 7214930- 17 . 6. 02 JW. 02  6  *247  F2.5  0  4  NH  N  «0  cRO  0  N  CS  0W'>*I2 

50  CloudS  01318tCDR  8WV  Priv  Co**  5781  S33  Sc  C  8  8  8  0  cP  U,*ClLO  1 

51  19:43  GIG. 7  2L177  6 1 05- 1 1 . 4 . 02 JW. 01  7  *  59  F8.61.224  NN  N  *<lcR0.02  N  CSS-*1-S2 

52  13:44  &2D  411L174  6295- 12 . 0 . 05 JW. 05  9  *145  F4.1  .5  m  NN  N  3.0cRO  m  H  CS  UmW*l2 

53  19:47  G50  2  bLl78  6355- 1 2 . 0  .  13 JW.  1 1 1 0  *219  »■  2  .  7  .3»7  NN  N  5.0cRa  m  N  CS  l)mW»l2 

54  :  i  49  G4A  3  7L172  63 1 0  - 1 1 . 6 .  03  JW .  OJ  7  ml/j  r-.;.3  ,  >»}  cRO  m  N  CS  UmW*I2 

S'l  .oud6  013131SNY  UWY  Priv  Comm  5/81  S42  s;  1  U  U  J  8  cP  U  I 

56  17:48  A1J.9  3L169  61 70 -  1 0 . 5 . 11 JW. 1 0  9  *225  FI. 8  .  4  m  NN  N  8  RO  m  N  CS  UmW»l2 

57  17:51  A2J  2  4L172  61 75 - 1 0 . 3 . 08JW . 1 0  9  *234  FI. 8  .3  m  NN  N  »2  cRO  m  N  CS  UmW*l2 

111111111  12222222222j3333333334444444444555S555555b66bbbbbb67  777  77777  78 
1 2  3456  7890 1 23456  7 390  1  2  345  6  7  8  90  1  2  3  4  5  6  7890  1  2.54 -,b  7890  I  23456  7 890  123456  7890 1234567890 


OUTSIDE  AIR  TEMPERATURE 


Sample  Application  of  the  NRL/FAA  Data  Base. 

FIGURE  C-3.  DATA  MlliS  JIOTED  IN  THE  NEW  DATA  BASE  FOR  INDICATED  INTERVALS  IN 
THE  VICINITY  OF  EACH  *FEST  POINT  IN  THE  ARMY  ICING  MATRIX.  Results  apply  to 
the  altitude  range  0-10,000  ft  AGL.  The  percentages  indicate  the  fraction  of 
all  data  miles  below  10,000  ft  AGL  that  were  found  In  the  Data  Base  for  each 
of  the  +2.5®C  and  +0.125  gm/m^  Intervals  centered  on  the  test  "points” 
shown  In  Fig.  C-1.  Of  6685  data  miles,  4625  (69%)  fell  somewhere  within  the 
4x4  boxed  area.  Of  the  2060  (31%)  data  miles  which  fell  outside,  1840 
(27.3%)  were  at  lesser  LWCs,  50  (0.9%)  were  at  greater  LWCs,  185  (2.7%)  were 
at  higher  temperatures,  and  24  (0.4%)  were  at  lower  temperatures. 


Sample  Application  of  the  NRL/FAA  Data  Base 


o 


TABLE  C-4.  Modern  Uata  for  OAT  <.  -17. ;>  degC  (.Row  ^  of  Army  Icing  Test  Matrix) 


For  the  Altitude  Interval  of  O-SOOU  Ft  A(>L 


Cloud  Geog  Air  Month  OAT 
File/Rec  Type  Loc.  Masa  of  Yr  degC 


Data 
Ml  lea 


Data  Meather 
Source  Situation 


.  12 

1 . 0 

.  08 

3.  0 

.  10 

6.  0 

.  Ob 

3.0 

.  12 

1.0 

.  17 

1.0 

.  1  / 

1 . 0 

.  16 

4.  0 

.  14 

17.0 

.  11 

20.0 

.  to 

S.  0 

.  13 

6.0 

.23 

2.0 

.  13 

2. 0 

.  13 

18.  0 

.07 

2.0 

.  U8 

12.0 

.  10 

20.0 

.  12 

2.0 

.  08 

2.0 

.17 

3.  0 

.  09 

2.0 

.  08 

2.0 

.  10 

3.  0 

.  08 

2.0 

.  06 

1.0 

.  07 

1.0 

.  0/ 

6.  0 

.  12 

7.0 

.  17 

6.  0 

.23 

1.0 

.  10 

3.  0 

.08  4.0 

Hp»toFFmCf , 156 
U.Mk  155 
IJ.Wk  155 
U ,Mk 155 
U,Wkl55 
0,Mk  155 
U,Mkl55 
U,Mk 155 
U,Mk  155 
U.Mkl55 
U.Mk 155 
U,Mk  155 
U,Mk  155 
U.Wk 155 
U,Wk 155 
U 

U.pClO 

U,pC10 

0 

U 

U 

0 

U 

0 

u 

U  ,  164 

U.  164 

U,  164 

U.  164 

•J  ,  1 64 

U,pClLO 

U,MkI58.pClL0 

U,MkI58,pClL0 

U,Mk 158,pCiL0 

U,Mk 158,pClL0 

U .Wk 158.pClL0 

U ,  pClLO 

U.pClLO 

U,pClLO 

U.pClLO 

U.pClLCl 


Sample  Application  of  the  NRL/FAA  Data  Base 


Modern  Data  for  LWC  >.875  g/m^  (Column  4  of  the  Army  Icing  Test 
Matrix)  in  the  Altitude  Interval  of  0-10,000  ft  AGL. 
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FIGURE  47.  APPROXIMATE  EXTREME  VALUES  OF  LWC  AND  MVD  (XJMBI NATIONS  OBSERVED 
IN  SUPERCOOLED  CLOUDS  AT  ALTITUDES  UP  TO  10,000  FT  AGL.  The  curved  lines 
here  represent  the  approximate  extreme  values  of  LWC  and  MVD  observed  in  any 
supercooled  cloud  icing  event  up  to  10,000  ft  ACL  and  up  to  Che  temperatures 
ind Icated . 
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Sample  Application  of  the  NRL/FAA  Data  Base. 
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FIGURE  38.  SCATTERPLOT  OF  MVD  VS.  OAT  FOR  MODERN  DATA  FROM  SUPERCOOLED  LAYER 
CLOUDS  (Sc,  Sc,  Na,  As,  Ac)  UP  TO  10,000  FT  ACL.  The  various  plocclng 
symbols  represent  different  data  sources  as  indicated  in  the  key.  The  size 
of  each  symbol  Is  proportional  to  its  statistical  weight  (l.e.,  Che  observed 
horizontal  extent  of  Che  associated  icing  event)  as  shown  by  the  scale  above 
the  graph.  The  center  of  each  symbol  corresponds  to  the  average  (and 
approximately  constant)  value  of  MVD  and  OAT  observed  during  Che  icing 
event.  The  solid  line  bounding  the  data  points  represents  the  apparent  upper 
limit  to  MVD  as  a  function  of  temperature  for  supercooled  layer  clouds  below 
10,000  ft  ACL,  The  position  of  the  line  at  temperatures  above  -IS’C  is  based 
on  the  maximum  MVDs  in  the  modem  data  only,  but  below  -IS'C  the  line  is 
based  on  maximum  MVDs  from  both  the  NACA  and  modem  data  sets.  The  data 
points  plotted  at  1  iim  MVD  are  those  for  which  the  MVD  values  are  actually 
unknown.  A  total  of  2660  data  miles  is  represented  in  this  graph. 
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Sample  Application  of  the  NRL/FAA  Data  Base 
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FIGURE  40.  SCATTERPLOT  OF  MVD  VS.  OAT  FOR  MODERN  DATA  FROM  SUPERCOOLED 
CONVECTIVE  CLOUDS  (Cu,  Cb)  UP  TO  10,000  FT  ACL.  The  various  plaCCXng  symbols 
represent  different  data  sources  as  Indicated  In  the  key.  The  »Lz«  of  each 
symbol  Is  proportional  to  Its  statistical  weight  (t.e.,  the  observed 
horizontal  extent  of  the  associated  Icing  event)  as  shown  by  Che  scale  above 
the  graph.  The  center  of  each  synbol  corresponds  to  the  average  (and 
approximately  constant)  value  of  MVD  and  OAT  observed  during  the  icing 
event.  The  solid  line  bounding  the  data  points  represents  Che  apparent  upper 
and  lower  limit  to  MVD  as  a  function  of  temperature  for  supercooled 
convective  clouds  below  10,000  ft  ACL.  The  position  of  the  line  Is  based  on 
extreme  MVD  values  In  both  the  NACA  and  modern  data  sets.  The  data  points 
plotted  at  I  4m  MVD  are  those  for  which  the  MVD  values  ace  actually  unknown. 

A  total  of  980  data  miles  Is  represented  In  this  graph. 
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SECOND  TRI-SERVICE  CLOUD  MODELING  WORKSHOP 
26-27  JUNE  1984 


DEPARTMENT  OF  DEFENSE 

Office  of  the  Under  Secretary  of  Defense, 
Research  and  Engineering 
Office  of  Environmental  and  Lite  Sciences 
Room  3D129 
The  Pentagon 
Washington,  D.C.  20301 

Attn:  CAPT.  Edward  Harrison 

Defense  Advanced  Research  Projects  Agency/DEU 
1400  Wilson  Boulevard 
Arlington,  VA  22209 

Attn:  Lt.  Col.  R.P.  Benedict 

Dr.  Louis  C.  Marquet,  Room  911 

Defense  Advanced  Research  Projects  Ayency/STO 
1400  Wilson  Boulevard 
Arlington,  VA  22209 

Attn:  Col.  Wesley  Kurowski,  Room  941 

Lt.  Col.  Howard  Stears 
Dr.  Anthony  J.  Tether 

Defense  Advanced  Research  Projects  Agency/TTO 
1400  Wilson  Boulevard 
Arlington,  VA  22209 

Attn:  Mr.  John  Hanson 

Joint  Cruise  Missiles/ASG 
Washington,  DC  20363 

Attn:  Lt.  Col.  M.I.  Bloom,  Room  12E45 

NORAD/J3W 

Peterson  Field,  CO  80914 

Attn:  Maj .  George  Davenport 
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Defense  Nuclear  Agency 
Washington,  DC  20305 

Attn;  Maj .  Martin  R.  Allerding,  RAAE 
Mr.  Patrick  A.  Crowley,  RAAE 
Dr,  H.C.  Eitz,  Jr.,  RAAE 
Mr.  Peter  W.  Lunn,  RAAE 
Lt.  Col.  William  J.  McKechney,  RAAE 
Dr.  Randy  Rohr,  Aerospace  systems 
Mr.  Kenneth  Schwartz 
Dr.  Jim  Sommers,  Aerospace  Systems 
Lt.  Col.  Gilbert  W.  Ullrich,  RAAE 
Dr.  Leon  Wittwer,  RAAE 


U.S.  ARMY 

Headquarters,  Department  of  the  Army 
Chemical  Engineering 
Pulaski  Building 
20  Massachusetts  Avenue,  N.W. 

Washington,  DC  20314 

Attn;  Dr.  Richard  Gomez,  DAEN-RDM,  Room  6203 

U.S.  Army  Engineer  Topographic  Laboratory 
ETL-GS-A 

Fort  Belvoir,  VA  22060 

Attn;  Dr.  Paul  C.  Dalrymple 

Commander/Director 

USA  Atmospheric  Sciences  Laboratory 
White  Sands  Missile  Range,  NM  88002 
Attn;  Dr.  Elton  P.  Avara 
Dr.  Andreas  Goroch 
Dr.  Melvin  C.  Heaps 


BMD  Syscom 
106  Wynn  Drive 
Huntsville,  AL  35807 

Attn;  John  W.  Bowman 

DARCOM 

Harry  Diamond  Laboratory 
2800  Powder  Mill  Road 
Adelphi,  MD  20783 

Attn;  Dr.  Zoltan  G.  Sztankay,  DELHD-RT-GB 

U.S.  Army  Missile  Command 
Redstone  Arsenal,  AL  35898 

Attn;  Mr.  George  P.  Drake,  DRSMI-RAS 

Dr.  Oskar  M.  Essenwanger,  DRSMI-RRA 


646 


U.S.  Army  Engineer  iopographic  Lao 
ETL-GS-A 

Fort  Belvoir,  VA  22060 

Attn;  Dr.  Donald  Dery 

U.S.  Army  Missile  Command 
P.O.  Box  1500 
Huntsville,  AL  35809 

Attn:  Mr.  Al  Carmichael,  BMD  ATC-T 

Dr.  W.  Davies 

Mr.  Jerry  Hagood,  Chief,  DARPA/STO  Project  Office 

Cold  Regional  Research  and  Evaluation  Laboratory 

Geophysical  Sciences  Branch 
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Hanover,  NH  03755 

Attn:  Dr.  Roger  Berger 

Atmospheric  Sciences  Laboratory 
Directed  Energy  Branch 
White  Sands  Missile  Range,  NM  80022 
Attn:  Mr.  Tom  Hall 

U.S.  Array  Combined  Arms  Center 
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Fort  Leavenworth,  KS  66027 

Attn:  Mr.  Elbert  S.  Kennedy 

S.  NAVY 

U.S.  Naval  Air  Development  Center 
Warminster,  PA  18974 

Attn:  Dr.  Alan  H.  Blumenthal,  Code  3011 

Mr.  Michael  R.  Hess,  Code  3011 

Naval  Air  Systems  C'^mmand 
Washington,  DC  20361 

Attn:  Mr.  Barry  Dillon,  Code  340 

Dr.  Paul  Twitchell,  AIR  330-G 
Dr.  Paul  Wyman,  AIR  5493 

Naval  Environmental  Prediction  Research  Facility 
Naval  Postgraduate  school  Annex 
Monterey,  CA  93943 

Attn:  CDR.  Donald  Hinsman 

Mr.  Roland  E.  Nagle 


Naval  Observatory 
OPNAV  952D 

34th  and  Massachusetts  Avenue,  N.W. 
Washington,  DC  20390 

Attn:  CAPT.  Carl  W.  Hoffman 

Naval  Ocean  Systems  Center 
San  Diego,  CA  92152 

Attn:  Virgil  R.  Noonkester 

Dr.  Juergen  Richter 
Dr.  Larry  B.  Stotts,  Code  8114 

Office  of  Chief  of  Naval  Operations 
34th  and  Massachusetts  Avenue,  N.W. 
Washington,  DC  20390 

Attn:  CDR.  Rene  E.  Gonzalez,  Jr.,  OP-952 

Dr.  Richard  W.  James,  OP-952D1 

Office  of  Chief  of  Naval  Operations 
Washington,  DC  20350 

Attn:  CDR.  Tom  Piwowar,  OP-986 

Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  VA  22217 

Attn:  Mr.  Robert  Abbey,  Code  422MM 

Dr.  William  F.  Cross,  Code  422PO 
Dr.  Douglas  J.  DePriest 
Mr.  James  H.  Hughes,  422AT 
Dr.  Edward  Wegman,  Code  411 
Dr.  Alan  Weinstein 

Naval  Research  Laboratory 
Washington,  DC  20395 

Attn:  Dr.  Timothy  Coffey 

Dr.  John  Hornstein,  Code  6520 
Dr.  Richard  K.  Jeck,  Code  4113 
Dr.  Vince  Noble,  Code  7910 
Dr.  Lothar  Ruhnke 
R.  Steinberg,  Code  6509 
Dr.  Edward  J.  Stone,  Code  6520 

Naval  Sea  Systems  Command,  PMS-405 
Washington,  DC  20362 

Attn:  LCDR  Stanley  Grigsby,  Code  24 

Mr.  Ralph  Rudkin 
Mr.  A.L.  Stoessel, 
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Attn:  CAPT.  M.B.  Hollinyer 

Naval  Surtace  Weapons  Center 
Dafiijren,  VA  22448 
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Robert  Cawley 
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Mr.  Abraham  Hirschman 

Ur.  Barry  S.  Katz 

Or.  Bernard  V.  Kessler,  R-42 

Carl  W.  Larson 

Mr.  Bonald  E.  Matlack 

Martha  Melnik,  R-42 

Mr.  Spyros  K.  Petropoulos 
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Cnina  Lake,  CA  93555 

Attn:  Mr.  Edward  J.  Bevan,  Code  3912 
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Headquarters,  Air  Weather  Service 
Scott  Air  Force  Base,  IL  62225 

Attn;  Maj .  Marvin  G.  Coleman,  DNXA 
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Attn;  Capt.  Georye  Fisher 

Air  Force  CMD/SA 

Kirtland  Air  Force  Base,  NM  87117 
Attn:  Lt.  Brent  J.  Cote 

Dr.  Bruno  Manz 
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Maj.  Jeffrey  S.  Schleher,  DOWA 


U.S.  Air  Force,  Continued 


Air  Force  Studies  and  Analysis 
Washington,  DC  20330 

Attn:  Lt.  Col.  Royer  Christensen,  SAZ 

Lt.  Col  Thomas  E.  Moriarty,  SAZ 

Headquarters,  Air  Force  Systems  Command 
Andrews  Air  Force  Base,  MD  20334 
Attn;  J.  Anderson 

Lt.  Col  Wallace  H.  Chaplin 
Capt.  G.  Anderson  White,  WER 

Air  Force  Weapons  Laboratory 
Kirtland  Air  Force  Base,  NM  87117 

Attn:  Lt.  Col.  Vernon  L.  Bliss,  WE 

Col  Gilbert,  AR3 
Capt.  Allen  Ronn ,  WE 
Alfred  L.  Sharp,  NTAT 
Darrell  Spreen,  ARE 

Wright  Aeronautical  Laboratory 
Wr ight-Patterson  Air  Force  Base,  OH  45433 
Attn:  Capt.  Michael  D.  Abel,  WEA 

Mr.  vern  Best,  AART-1 
Vince  Diehl,  AART-2 
Maj .  Wayhe  Harper,  WE 
Ralph  Nieter,  AART-1 
David  Rardin,  AARAT-1 
Drexel  G.  Waggoner,  AARI-3 


OTHER  GOVERNMENT 

NASA  Headquarters 
Washington,  DC  20546 

Attn:  Dr.  Robert  Schifter,  Code  EE 

Dr.  Shelby  Tiltord 

NASA  Ames  Research  Center  -  Unclassified  Version  Onl 
Space  Science  Division  245-3 
Mountain  View,  CA  94035 

Attn:  Dr.  Edwin  Danielsen 

NASA  Goddard  Space  Flight  Center 
Greenbelt,  MD  20771 

Attn:  Dr.  Albert  Arkiny 

Dr.  Robert  Fraser 

Dr.  Paul  Hwang  -  Unclassified  Version  Onl 
Dr.  Joanne  Simpson,  Code  914 


other  Government,  Continued 

NOAA  Aeronomy  Lab/ERL 
325  Broadway 
Boulder,  CO  80303 

Attn:  Dr.  Freeman  Hall,  R4X2 

NOAA/NESDIS 

Federal  Office  Building  4 
Suitland  and  Silver  Hill  Roads 
Washington,  DC  20233 

Attn:  Dr.  Herbert  Jacobowitz 

Dr.  L.L.  Stowe 
Dr.  Harold  Yates 

NOAA  Air  Resources  Laboratory 
Eighth  Floor 
6010  Executive  Building 
Rockville,  MD  20852 

Attn:  Dr.  Lester  Machta 

NOAA  National  Climate  Program 
Room  108 

11400  Rockville  Pike 
Rockville,  MD  20852 

Attn:  Mr.  Martin  C.  Yerg 


INDUSTRY 

Aerodyne  Research 
45  Manning  Road 
Billerica,  MA  01821 

Attn:  Dr.  Morton  Camac 

Dr.  H.J.  Caulfield 
Dr.  Robert  B.  Lyons 

Aerojet  Electrosystems 
1100  West  Hollyvale  Street 
P.O.  BOX  296 
Azusa,  CA  91702 

Attn:  Mr.  Alain  L.  Fymat 

Mr.  Bernard  Lichtenste 
Ms.  Rosia  Y.  Pan,  Dept 
Mr.  Amiel  Shulsinger 
Ms.  Edna  M.  Sugihara 
Ms.  Ellen  S.  Thomas 

Aerojet  General  Corporation 
P.O.  Box  13222 
Sacramento,  CA  95813 

Attn:  Mr.  Thomas  G.  Lee 


,  ASB 


i  n 

.  4211,  Bldg.  160 


Industry,  Continued 

AeroMet  Inc. 

P.O.  Box  701767 
Tulsa,  Oklahoma  74170 
Attn:  Mr.  Ed  Gibeau 

Analytics,  Scientific  and  operations  Analysis  Group 
7680  Old  Sprinyhouse  Road 
McLean,  VA  22102 

Attn:  Dr.  Kenneth  D.  Shere 

Atlantic  Analysis  Corporation 
5  Koger  Executive  Center 
Suite  219 

Norfolk,  VA  23502 

Attn:  Mr.  Roger  S.  Staff 

Atmospheric  and  Environmental  Research  Inc. 

840  Memorial  Drive 
Cambridge,  MA  02139 

Attn:  Mr.  Ronald  G.  Isaacs 

Boeing  Aerospace  Company 
P.O.  BOX  3707 
Seattle,  WA  98124 

Attn:  Mr.  Richard  Geer 

Dr.  Martin  Savol 

Boeing  Aerospace  Company 
P.O.  Box  3999 

Seattle,  Washington  98124 

Attn:  Mr.  Tom  Sangston,  Mail  Stop  8H37 

Calspan  Corporation 
P.O.  BOX  400 
Buffalo,  NY  14225 

Attn:  Dr.  Paul  Marrone 

Control  Data  Corporation 
Research  Division 
P.O.  BOX  1249 
Minneapolis,  MN  55440 

Attn:  Mr.  Greg  Nastrom 


Environmental  Research  Institute  of  Michigan 

P.O.  Box  618 

Ann  Arbor,  MI  48107 

Attn:  Dr.  Anthony  J.  Larocca 


Industry,  Continued 


General  Dynamics 
Mail  Zone  2814 
P.O.  Box  748 
Fort  Worth,  TX  76101 

Attn;  Mr.  Tommy  R.R.  Augustsson 
Mr.  John  A.  Guinn,  jr. 

General  Dynamics 
P.O.  Box  2507 
Pomona,  CA  91769 

Attn:  Mr.  Richard  M.  Pietrasz/MZ  600-7 

General  Dynamics 
Corvair  Division 
P.O.  Box  80847  ■ 

San  Deigo,  CA  921J8 

Attn:  Mr.  J.S.  Toor,  MZ  41-6850 

General  Electric 

Aerospace  Electronic  Systems  Department 
French  Road 
Utica,  NY  13503 

Attn:  Mr.  Phillip  W.  Kane,  Mail  Drop  39 

General  Electric  -  AEBG 
1  Neumann  Way 
Cincinnati,  OH  45215 

Attn;  Mr.  S.D.  Bertke,  Mail  Drop  J-185 

General  Electric  -  AESD 
901  Broad  street 
Utica,  NY  13503 

Attn;  Mr.  Lome  H.  Ford,  Mail  Drop  900 
Mr.  Philip  W.  Kane,  Mail  Drop  900 

General  Electric  Company 
P.O.  Box  4840 
Syracuse,  NY  13221 

Attn;  Mr.  Edwin  L.  Post,  CSP-5,  Room  G7 

General  Motors  Research  LaD 
Physics  Department 
Warren,  MI  48090 

Attn:  Dr.  Ruth  Reck 

General  Research  Corporation 
Westgate  Research  Park 
7655  Old  Springhouse  Road 
McLean,  VA  22102 

Attn;  Mr.  John  L.  Gilbert 
Mr.  Robert  Poppe 
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General  Research  Corporation 

P.O.  Box  6770 

Santa  Barbara,  CA  93111 

Attn:  Mr.  Richard  Rein 

Mr.  Tim  Sullivan 

Grumman  Aerospace  Corporation 
R&D  Center,  Plant  26 
Bethpaye,  NY  11714 

Attn:  Mr.  Gil  Kelley,  MS  B06-04 

Dr.  Jerry  Krassner 
Dr.  John  E.A.  Selby 

Honeywell  EOD 
2  Forbes  Road 
Lexinyton,  MA  02173 

Attn:  Mr.  Raymond  Gelinas 

Honeywell ,  Inc. 

Minneapolis,  MN  55440 

Attn:  Dr.  Joseph  E.  Carroll,  MN  17-2328 

IRT  Corporation 
6800  poplar  Place 
McLean,  VA  22101 

Attn:  Dr.  David  M.  Harris 

Mr.  David  p.  Woodall 

ITT  Avionics  Division 
390  Washinyton  Avenue 
Nutley,  NJ  07110 

Attn;  Mr.  Robert  J.  Kraushaar 
Mr.  Kevin  B.  Ward 

Kaman  Sciences 

1500  Garden-of-the-Gods  Road 
Colorado  Springs,  CO  80933 

Attn:  Mr.  Jerry  L.  Harper 

Mr.  Jacob  j.  Hess 
Mr.  Wayne  j.  young 

Lockheed  Missiles  and  Space  Company 
0/92-20,  B/205 
3251  Hanover  street 
Palo  Alto,  CA  94304 

Attn:  Mr.  scott  Claflin 

Dr.  Georye  Gal 
Mr.  Norman  Kulgein 
Dr.  Alexander  Maksymowicz 
Mr.  Wayne  Rudolf 
Mr.  T.  Winarske 
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Lockheed  Missiles  and  Space  Company 
P.O.  Box  56 
Sunnyvale,  CA  94086 

Attn:  Mr.  Lee  G.  Dickinson,  Dept.  85-72 

Dr.  Rich  Hake,  Dept.  62-23 

Lockheed  California  Compnay 
P.O.  Box  551 
Burbank,  CA  91520 

Attn:  Mr.  A.C.  Brown 

Loton  Research 

3377  North  Torrey  Pines  Court 
La  Jolla,  CA  92037 

Attn:  Mr.  David  M.  Zimmerman 

McDonnell  Douglas 
A3-232-i3-3 

McDonnell  Douglas  Astronautics  Company 
5301  Boisa  Avenue 
Huntington  Beach,  CA  92647 

Attn:  Mr.  Robert  A.  Willett 

Ontar  Corporation 
129  University  Road 
Brookline,  MA  02146 

Attn;  Mr.  John  W.  Schroeder 

The  Optical  sciences  Company 
P.O.  Box  1329 
Placentia,  CA  92670 

Attn:  Dr.  David  Fried 

Dr.  David  L.  Hench 

The  Optical  Sciences  Company 
8301  Greensboro  Drive 
McLean,  VA  22102 

Attn:  Dr.  Thomas  F.  Wiener 


OptiMetrics  Inc. 

P.O.  Box  7766 

Ann  Arbor,  MI  48107 

Attn:  Mr.  William  O.  Gallery 

Dr.  Robert  E.  Meredith 
Mr.  Merle  J.  Persky 
Mr.  Frederick  G.  Smith 


Industry,  Continued 


Pacific-Sierra  Research  Corporation 
Suite  liOO 

1401  Wilson  Boulevard 
Arlington,  VA  22209 

Attn:  Dr.  Ronald  N.  Dewitt 

Pacific-Sierra  Research  Corporation 
12340  Santa  Monica  Boulevard 
Los  Angeles,  CA  90025 

Attn:  Mr.  Alan  Shapiro 

Perkin-Elmer  Company 
100  Wooster  Heights  Road 
Danbury,  CT  06810 

Attn:  Dr.  Robert  Arguello,  MS  845 

PhotoMetrics ,  Inc. 

4  Arrow  Drive 
Woburn,  I«1A  01801 

Attn:  Dr.  Irving  L.  Kofsky 

Photon  Research  Associates 
2223  Avenida  De  La  Playa 
Suite  301 

La  Jolla,  CA  92037 

Attn:  Dr.  David  Anding 

Dr.  William  Cornette 

Physical  sciences  Inc. 

P.O.  BOX  3100 
Research  Park 
Andover,  MA  01810 

Attn:  Dr.  Guy  Weyl 

Radiation  Research 
3550  Hulen  Street 
Fort  Worth,  TX  76107 

Attn:  Dr.  M.B.  Wells 

Rockwell  international 
11642  Wallingsford  Road 
Los  Alamitos,  CA  90720 

Attn:  Mr.  Michael  F.  Sentovich 

SRI  International 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94061 

Attn:  Dr.  J.H.  Allen 

Mr.  John  D.  Malick 

Mr.  Richard  H.  Monahan,  Room  AH-157 
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Science  Applications  Inc. 

P.O.  Box  1303 
McLean,  VA  22101 

Attn:  Dr.  John  Connolly 

Science  Applications  Inc. 

505  Marquette,  N.W. 

AlDuquerque,  NM  87102 

Attn:  Dr.  James  L.  Griggs,  Jr, 

Mr.  Ronald  J.  Nelson 

Science  Applications  Inc. 

P.O.  Box  2351 
La  Jolla,  CA  92083 

Attn;  Mr.  Dan  Hamlin 

Science  Applications  Inc. 

1010  Woodman  Avenue 
Dayton,  OH  45432 

Attn;  Mr.  Bob  Turner 


Scientific  Simulation  Inc. 

P.O.  Box  9331 
Albuquerque,  NM  87119 

Attn:  Dr.  L.R.  Ebbesen 

Mr.  Don  C.  Fronterhouse 

Spectral  Sciences  Inc. 

Ill  South  Bedford  Street 
Burlington,  MA  01803 

Attn:  Dr.  Lawrence  Bernstein 

Dr.  David  C.  Robertson 

Systems  and  Applied  Sciences  Corporation 
5809  Annapolis  Road 
Hyattsville,  MD  20784 

Attn;  Dr.  P.K.  Bhartia  -  Unclassified  Version  Onl 

Systems  and  Applied  Sciences  Corporation 
109  Massachusetts  Avenue 
Lexington,  MA  02173 

Attn;  Dr.  Ralph  Shapiro 
Mr.  C.N.  Touart 

TASC 

1  Jacob  Way 
Reading,  MA  01867 

Attn:  Dr.  Robert  F.  B rammer 

Dr.  Paul  Janota 
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TASC 

8301  Greensboro  Drive 
McLean,  VA  22102 

Attn:  Dr.  Thomas  Wiener 


TRW 

One  Space  Park 
Building  82/1710 
Redondo  Beach,  CA  90278 
Attn:  Mr.  bill  Momary 

Teledyne  Brown  Engineering 
300  S.  Parkman  Drive 
MS-19 

Huntsville,  AL  358-07 

Attn:  Mr.  Robert  A.  Pilgrim 

Mr.  Mike  Scarboro 


Titan  Systems 
P.O.  BOX  12139 
La  Jolla,  CA  92037 

Attn:  Jeffrey  j.  Puschell 


Verac,  Inc. 

2901  Juan  Tabo,  N.E. 

Suite  235 

Albuquerque,  NM  87112 

Attn:  Mr.  Joel  S.  Davis 

Mr.  Jack  Davis 
Mr.  Bryan  Thompson 


Verac,  Inc. 

10601  Lomas  Boulevard,  N.E. 
Suite  115 

Albuquerque,  NM  87112 

Attn:  Dr.  T.M.  Bomber 


Vis idyne 

5  Corporate  place 
South  Bedford  Street 
Burlington,  MA  01803 

Attn:  Dr.  Charles  H.  Humphrey 

Dr.  Henry  Smith 
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Aerospace  Corporation 
2350  El  Seyundo  Boulevard 
El  seyundo,  CA  90009 


B 

Attn;  Mr. 

Roy  M.  Chiulli 

Mr. 

David  Escoe 

Mr. 

Eric  Jensen,  Ml/041 

Mr. 

Paul  Kisliuk 

r  ' 

Dr. 

Huyh  Ruyye 

L_ 

Dr. 

Fred  Simmons 

9 

Colorado  State 

University 

Fort  Collins,  CO  80523 

Attn;  Dr.  Tom  Brubaker,  Electrical  Enyineeriny  Dept. 
Unclassified  Version  Only 

Prof.  Thomas  Vondernaar,  Atmospheric  Sciences  Dept 


f  , 


Jet  Propulsion  Laboratory 
4800  Oak  Grove  Drive 
Pasadena,  CA  91109 

Attn;  Dr.  Alex  Goetz 
Dr.  E.  Hinkley 

Johns  Hopkins  (APL) 

Johns  Hopkins  Road 
Laurel,  MD  20777 

Attn;  William  J.  Tropf 
Randolph  W.  Bruns 

Johns  Hopkins  (APL) 

9713  Clocktower  #202 
Columbia,  MD  21046 

Attn;  Andrew  N.  Vavreck 

Lawrence  Livermore  National  Laboratory 
P.O.  Box  808 
Livermore,  CA  94550 

Attn;  Dr.  Joseph  B.  Knox 
Dr.  M.  McCracken 

Los  Alamos  National  Laboratory 
International  Technical  Office 
P.O.  Box  503 
Los  Alamos,  NM  87545 

Attn;  Dr.  Sumner  Barr 

Dr.  Henry  G.  Horak 
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MIT  Lincoln  Laboratory 
P.O.  box  73 
Lexington,  MA  02173 

Attn:  Dr.  Herb  Kleiman 

Mr.  John  T.  Prohaska,  KB  237 

RAND  Corporation 
1700  Main  Street 
Santa  Monica,  CA  90406 

Attn:  Mr.  Ralph  Huschke 

Scripps  Institute  tor  Oceanography 
Visibility  Laboratory 
University  of  California,  San  Diego 
La  Jolla,  CA  92093 

Attn:  Mr.  Wayne  Hering,  P-003 

Mr.  Richard  W.  Johnson 

University  of  Maryland 

Institute  of  Physical  Science  and  Technology 
College  Park,  MD  20742 

Attn:  Dr.  Thomas  Wilkerson 

University  of  Washington 
Seattle,  WA  98195 

Attn:  Prof.  Peter  Hobbs 


IDA 

Dr.  Ernest  Bauer,  STD 
Mr.  Lucien  Biberman,  STD 
Dr.  Erwin  Kaufman,  SED 
Ms.  Elizabeth  Li 
Dr.  David  L.  Randall,  JED 
Dr.  Robert  C.  Oliver,  STD 
Dr.  Hans  Wolf hard,  STD 
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Commander /Director 

USA  Atmospheric  Sciences  Laboratory 
WSMR,  NM  88002 

Attn:  Dr.  Elton  P.  Avara 

BMD  Syscom 
106  Wynn  Drive 
Huntsville,  AL  35807 

Attn:  John  W.  Bowman 

DARCOM 

Harry  Diamond  Labs 
2800  Powder  Mill  Road 
Adelphi,  MD  20783 

Attn:  Dr.  Zoltan  G.  Sztankay,  DELHD-RT-LB 

U.S.  Army  Missile  Command 
Redstone  Arsenal,  AL  35898 

Attn:  Mr.  George  P.  Drake,  DRSMI-RAS 

Dr.  Oskar  M.  Essenwanger,  DRSMI-RRA 

S.  NAVY 

U.S.  Naval  Air  Development  Center 
Warminster,  PA  18974 

Attn:  Dr.  Alan  H.  Blumenthal,  Code  3011 

Mr.  Michael  R.  Hess,  Code  3011 

Naval  Air  Systems  Command 
AIR  330G 

Washington,  DC  20361 

Attn:  Dr.  Paul  Twitchell,  AIR  330-G 

Dr.  Paul  Wyman,  AIR  5493 

Naval  Environmental  Prediction  Research  Facility 
24  Greenwood  Vale 
Monterey,  CA  93943 

Attn:  Dr.  Roland  E.  Nagle 

Dr.  Paul  M.  Tag 


U.S.  Navy,  Continued 


Naval  Ocean  Systems  Center 
San  Diego,  CA  92152 

Attn;  Virgil  R.  Noonkester 

Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  VA  22217 

Attn:  Douglas  J.  DePriest 

Mr.  James  H.  Hughes,  422AT 

Naval  Research  Laboratory 
Washington,  DC  20395 

Attn:  Dr.  John  Hornstein,  Code  6520 

Dr.  Richard  K.  Jeck,  Code  4113 
Dr.  Edward  J.  Stone,  Code  6520 

Naval  Sea  Systems  Command 
Washington,  DC  20362 

Attn;  LCDR  Stanley  Grigsby,  PMS-405 

Naval  Surface  Weapons  Center 
Dahlgran,  VA  22448 

Attn;  Ms.  Kathleen  J.  Fairfax 
Ralph  J.  Fallin 
Susan  L.  Masters,  K-44 

Naval  Surface  Weapons  Center 
White  Oak 

Silver  Spring,  MD  20918 

Attn:  Barton  D.  Billard 

W.E.  Caswell,  R-42,  90-001D 
Robert  Cawley 

Ms.  Zarak  Hanks,  K-44,  1200 
Dr.  Bernard  V.  Kessler,  R-42 
Carl  W.  Larson 
Martha  Melnik,  R-42 

Commander  Naval  Weapons  Center 
China  Lake,  CA  93555 

Attn:  Mr.  Edward  J.  Bevan,  Code  3912 

Donald  A.  Kappelman,  Code  39403 


AIR  FORCE 


Headquarters,  Air  Weather  Service 
Det .  1 

The  Pentagon 
Washington,  DC  20330 

Attn;  Lt.  Col.  Randolph  W.  Ashby 
Capt.  Lauraleen  O'Connor 
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U.S.  Air  Force  Environmental  Technical  Applications  Center 
Scott  Air  Force  Base,  IL  62225 

Attn:  Dr.  Patrick  J.  Breitling,  CCN 

Lt.  Col.  Pershing  Hicks,  Jr.,  DN 

Air  Force  Geophysical  Laboratory 
Hanscom  Air  Force  Base,  MA  01731 

Attn:  Mr.  James  T.  Bunting,  LYS 

Dr.  Robert  Fenn,  OPA 
Mr.  Donald  Grantham,  LYT 
Mr.  Irving  I.  Gringorten,  LYT,H 
Mr.  John  H.  Schummers,  LSA 

U.S.  Air  Force  Space  Division 
P.O.  Box  92960 
Worldway  Postal  Center 

Los  Angeles  Air  Force  Station,  CA  90009 
Attn:  Edmund  D.  Daszewski,  WE 

Mark  W.  Lindsey 
Todd  M.  Niepke 

Lt.  Col.  Edward  M.  Tomlinson,  DAAX 

Headquarters,  Air  Force  Systems  Command 
Andrews  Air  Force  Base,  MD  20334 

Attn:  Capt.  G.  Anderson  White,  WER 

Air  Force  Weapons  Laboratory 
Kirtland  Air  Force  Base,  NM  87117 

Attn:  Lt.  Col.  Vernon  L.  Bliss,  WE 

Alfred  L.  Sharp,  NTAT 

Air  Force  Global  Weather  Central/TSIT 
Offutt  Air  Force  Base,  NE  68113 
Attn:  Lt.  Col.  William  M.  Cox 

Wright  Aeronautical  Laboratory 
Wright-Patterson  Air  Force  Base,  OH  45433 
Attn:  Capt  Michael  D.  Abel,  WEA 

Drexel  G.  Waggoner,  AARI-3 

OTHER  GOVERNMENT 

National  Climate  Program 

NOAA 

Room  108 

11400  Rockville  Pike 
Rockville,  MD  20852 

Attn:  Martin  C.  Yerg 
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Aerodyne  Research 
45  Manning  Road 
Billerica,  MA 

Attn:  Dr.  Robert  B.  Lyons 

Aerojet  Electrosystems 
1100  West  Hollyvale  Street 
P.O.  Box  296 
Azusa,  CA  91702 

Attn:  Alain  L.  Fymat 

Edna  M.  Sugihara 

Aerojet  General  Corporation 
P.O.  Box  13222 
Sacramento,  CA  95813 

Attn:  Thomas  G.  Lee 

Aerospace  Corporation 
2350  El  Segundo  Boulevard 
El  Segundo,  CA  90009 

Attn;  Roy  M.  Chiulli 

Aerospace  Corporation 
Mail  Station  Ml/041 
P.O.  Box  92957 
Los  Angeles,  CA  90009 

Attn;  Dr.  Eric  B.  Jensen 

Atmospheric  Sciences  BV/NOAA 
NESDIS 

5001  Silver  Hill  Road 
Washington,  DC  20233 

Attn:  Dr.  Herbert  Jacobowitz 

Atlantic  Analysis  Corporation 
5  Roger  Executive  Center 
Suite  219 

Norfolk,  VA  23502 

Attn;  Roger  S.  Staff 

General  Dynamics 
Mail  Zone  2814 
P.O.  Box  748 
Fort  Worth,  TX  76101 

Attn:  Tommy  R.R.  Augustsson 

General  Dynamics 
P.O.  Box  2507 
Pomona,  CA  91769 

Attn:  Richard  M.  Pietrasz/MZ  600-7 
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General  Electric 

Aerospace  Electronic  Systems  Department 
French  Road 
Mail  Drop  39 
Utica,  NY  13503 

Attn:  Phillip  W.  Kane 

General  Research 
Westgate  Research  Park 
McLean,  VA  22102 

Attn:  John  L.  Gilbert 

Grumman  Aerospace 
Plant  26 
R&D  Center 
Bethpage,  NY  11714 

Attn:  Dr.  Jerry  Krassner 

ITT  Corporation 

390  Washington  Avenue 

Nutley,  NJ  _ 

Attn:  Mr.  Robert  J.  Kraushaar 

Kaman  Sciences 

1500  Garden-of-the-Gods  Road 
Colorado  Springs,  CO  80933 

Attn:  Mr.  Jerry  L.  Harper 

Lockheed 
0/92-20,  B/205 
3251  Hanover  Street 
Palo  Alto,  CA  94304 

Attn:  Dr.  George  Gal 

Alexander  Maksymowicz 

Lockheed  Missiles  and  Space  Company 

P.O.  Box  56 

Department  85-72 

Building  154 

Sunnyvale,  CA  94806 

Attn:  Mr.  Lee  G.  Dickinson 

Loton  Research 

3377  North  Torrey  Pines  Court 
La  Jolla,  CA  92037 

Attn:  David  M.  Zimmerman 
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McDonnell  Douglas 
A3-232-13-3 

McDonnell  Douglas  Astronautics  Company 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 
Attn:  Robert  A.  Willett 

Ontar  Corporation 
129  University  Road 
Brookline,  MA  02146 

Attn:  Mr.  John  W.  Schroeder 

The  Optical  Sciences  Company 
P.O.  Box  1329 
Placentia,  CA  92670 

Attn:  Dr.  David  L.  Hench 

The  Optical  Sciences  Company 
8301  Greensboro  Drive 
McLean,  VA  22102 

Attn:  Dr.  Thomas  F.  Wiener 

OptiMetrics 

P.O.  Box  7766 

Ann  Arbor,  MI  48107 

Attn:  Mr.  Frederick  G.  Smith 

Pacific-Sierra  Research 
Suite  1100 

1401  Wilson  Boulevard 
Arlington,  VA  22209 

Attn:  Dr.  Ronald  N.  DeWitt 

PhotoMetrics ,  Inc. 

4  Arriw  Drive 
Woburn,  MA  01801 

Attn;  Dr.  Irving  L.  Kofsky 

Radiation  Research 
3550  Hulen  Street 
Fort  Worth,  TX  76107 

Attn:  Dr.  M.B.  Wells 

Rockwell  International 
11642  Wallingsford  Road 
Los  Alamitos,  CA  90720 

Attn;  Michael  F.  Sentovich 


Industry,  Continued 

SRI  International 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94061 

Attn;  John  D.  Malick 

Richard  H.  Monahan,  Room 

Science  Applications 
505  Marquette,  N.W. 

Albuquerque,  NM  87102 

Attn:  Dr.  James  L.  Griggs,  Jr. 

Mr.  Ronald  J.  Nelson 

Spectral  Sciences 

111  South  Bedford  Street 

Burlington,  MA  01803 

Attn:  Dr.  David  C.  Robertson 

Systems  Applied  Sciences  Corporation 
109  Massachusetts  Avenue 
Lexington,  MA  02173 

Attn;  Dr.  Ralph  Shapiro 

Teledyne  Brown 

300  S.  Parkman  Drive 

MS-19 

Huntsville,  AL  35807 

Attn;  Mr.  Robert  A.  Pilgrim 

Titan  Systems 
P.O.  Box  12139 
La  Jolla,  CA  92037 

Attn:  Jeffrey  J.  Puschell 

Verac,  Inc. 

2901  Juan  Tabo,  N.E. 

Albuquerque,  NM  87112 

Attn:  Mr.  Joel  S.  Davis 

Visidyne 

5  Corporate  Place 
South  Bedford  Street 
Burlington,  MA  01803 

Attn;  Dr.  Charles  H.  Humphrey 
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Mr.  Paul  Janota 
One  Jacob  Way 
Reading,  MA  01867 

University  of  California,  San  Diego 

Vislab  P-003 

San  Diego,  CA  92023 

Attn;  Mr.  Wayne  S.  Hering 
Attn:  Mr.  Richard  W.  Johnson 

Johns  Hopkins  (APL) 

Johns  Hopkins  Road 
Laurel,  MD  20777 

Attn:  William  J.  Tropf 

Johns  Hopkins  (APL) 

9713  Clocktower  #202 
Columbia,  MD  21046 

Attn:  Andrew  N.  Vavreck 
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TO:  AFWAL/AART-2 


1.  Activity  and  Place  Visited:  Naval  Surface  Weapons  Center,  White  Oak, 
Maryland. 


2.  Departure  and  Return  Date:  25-28  June  1984. 

3.  Purpose  of  Visit:  Second  Tri-Service  Cloud  Modeling  Workshop. 


4.  Persons  Contacted:  A  complete  list  of  attendees  will  be  published  with 
the  Workshop  Proceedings.  In  addition  to  the  presenters  and  session  chair¬ 
men  listed  in  the  Agenda  (Atch  1),  I  have  made  up  a  partial  listing  of 
attendees  (Atch  2). 


5.  Factual  Data:  As  stated  above,  the  meeting  agenda  is  given  in  Atch  1. 
There  was  a  wide  diversity  of  papers  presented.  All  three  Services  were 
well  represented.  The  work  being  done  today  seems  to  be  sponsored  by 
three  main  development  efforts:  space  or  ground  based  lasers,  infrared 
search  and  track,  and  ICBM  re-entry  (Navy).  Of  course,  better  understanding 
of  clouds  and  cloud  modeling  can  benefit  a  much  broader  group  of  systems/ 
development  efforts  including  CMAG,  MICOS,  IRST,  Night-in-Weather,  and 
HAVE  LACE,  just  to  mention  a  few  Avionics  Laboratory  programs.  I  have  some 
notes  on  all  but  the  classified  papers.  Atch  3  is  a  very  good  list  of 
Cloud  Free  Line-of-Sight  references  given  out  at  the  Workshop. 


6.  Conclusions :  The  Cloud  Modeling  Workshop  is  relatively  new  and  there 
are  perhaps  three  general  areas  where  problems  exist  (as  I  see  it):  The 
leadership  issue  is  probably  the  least  troubling.  A  number  of  Navy,  Army, 
and  Air  Force  laboratories  have  supported  this  effort  from  its  beginning 
in  late  1981.  Currently,  most  leadership  is  coming  from  the  Air  Force 
Geophysics  Laboratory.  Next  year's  workshop  is  being  scheduled  to  follow 
the  Tri-Service  Transmission  Conference  at  AFGL;  however,  when  compared 
with  the  atmospheric  transmission  area,  "clouds"  do  not  have  the  same 
support  and  direction  from  the  Pentagon. 

The  Workshop  was  most  concerned  about  the  money  issue. 

A  number  of  DoD  systems  are  affected  (but  no  one  knows  to  what  extent)  by 
clouds.  Unfortunately,  no  one  wants  to  fund  generic  research  in  this  area 
(it  had  been  hoped  that  DARPA  money  would  be  available).  There  is  some 
money  to  find  answers  to  specific  questions,  but  this  limited  effort  may 
not  really  be  enough  to  develop  good  cloud  models  for  realistic  evaluation 
of  weapon  system  performance.  The  laboratories  in  the  past  may  have  taken 
the  lead  and  funded  such  work  out  of  their  own  monies,  but  today  they  are  in 
a  'contractor- type'  position;  paid  by  specific  development  projects  to  do 
specific  tasks.  Somehow,  the  Workshop  will  have  to  design  its  efforts  to 
work  efficiently  within  this  less  than  desirable  financial/management 
framework.  Finally,  there  is  the  issue  of  what  should  the  Workshop  be 
doing  (besides  meeting  once  a  year)?  Below  are  listed  some  tasks  brought 
up  at  the  meeting  (which  I  endorse): 


a.  Search  and  document  DoD  user  requirements.  First,  locate  all 
possible  customers  and  then  rank  their  possible  cloud  questions/problems 
by  importance  using  agreed  upon  criteria. 

b.  Define  as  best  as  possible  the  kind  of  cloud  model(s)  that  is 
needed  and  the  kind  of  data  necessary  to  support  its  development  and 
validation. 

c.  Survey  and  document  existing  models,  data  sets,  and  data  measure 
ment  programs.  Identify  shortfalls  based  on  items  (a)  and  (b). 

d.  Set  down  standards  for  data  measurements  and  archiving. 

e.  Promote  the  exchange  of  new  data  processing  technologies. 

7.  Recommendations :  Staffmet  involvement  with  this  Workshop  can  benefit 
many  AFWAL  programs  including  CMAG.  We  need  to  try  to  define  the  type  of 
cloud  model  and  data  that  would  be  most  useful  for  our  analysis  needs. 

8.  Actions  to  be  Taken:  Coordinate  with  Mr  Ron  Kaehr,  CMAG  Program 
Manager,  and  any  other  interested  program  manager. 

MICHAEL  D.  ABEL,  Captain,  USAF  3  Atch 

Staff  Meteorologist  1.  Agenda 

Avionics  Laboratory  2.  Some  Attendees 

3.  CFLOS  References 

cc:  AFWAL/AS 
AFWAL/WE 
AFWAL/AARF-1 
AFWAL/AAX 


Joel  S.  Davis 
Verac,  Inc. 

Albuquerque,  NM 

ATTACHMENT  I 

COMMENTS  FOR  THE  MINUTES:  2nd  TRI-SERVICE  CLOUD  MODELING  WORKSHOP 

Much  of  the  conversation  following  the  formal  presentations  made  the 
point  that  there  wasn't  a  lot  of  funding  available  for  data  collection  and 
analysis  related  to  cloud  problems,  specifically  tor  the  development  and 
validation  of  cloud  free  1 ine-of-sight,  cloud  free  field-of-view,  and  cloud 
free  interval  models.  Barring  a  major  change  in  policies  seen  to  date,  it 
is  unlikely  that  major  funding  for  such  investigations  will  ever  be 
available.  Cloud  analysis  and,  in  general,  weather  analysis  is  rarely 
funded  in  isolation  because,  as  a  problem,  it  is  not  important  in  isolation. 
Rather,  weather  is  generally  one  of  many  problems  facing  an  end  user.  Thus, 
resources  to  address  weather  issues  will  typically  have  to  be  taken  from  the 
resources  of  some  larger  investigation  or  project.  In  the  systems 

evaluations  and  concept  analyses  with  which  I  am  familiar,  this  typically 
translates  into  a  few  man  months,  at  most,  available  for  analyses 
specifically  related  to  the  weather  aspects  of  system  performance. 

Given  only  a  few  man-months  for  specifically  addressing  weather  issues, 
the  scientists  and  analysts  involved  are  very  limited,  both  in  time  and 

resources,  in  the  data  and  models  they  can  run  down.  Unfortunately,  the 

data  bases  and  models  needed  for  such  analyses  are  in  different  places, 
typically  in  different  media  and  formats,  and  are  not  well  publicized.  In 
addition,  even  where  they  are  publicized,  there  are  often  no  clear  or 
convenient  channels  for  obtaining  them. 

Consider  the  problem  of  developing  a  dynamic  cloud  free  1  ine-of-sight 
model  based  on  available  data.  Ideally,  such  a  model  would,  for  a  given 
type  of  engagement,  predict  dynamic  cloud  free  1  ine-of-sight  probabilities 
for  cloud  conditions  as  specified  in  a  normal  weather  report.  One  way  this 
might  be  done  would  be  to  take  satellite  based  cloud  pictures  with 
associated  altitude  information,  digitize  these  pictures,  determine  cloud 

free  intervals  and,  for  a  great  many  such  pictures,  in  a  great  many  such 
locations,  for  a  great  many  types  of  intervals,  correlate  them  with  the 
weather  reports  at  the  locations  measured. 
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This  approach  requires  the  weather  reports  including  clouds  from 
locations  of  interest:  surface  weather  reports  or  the  3-D  nephanalysis  or  RT 
nephanalysi s .  These  data  bases  are  probably  among  the  best  known  of  those 
available  within  the  community,  but  learning  to  read  them  and  use  them  is  no 
small  task.  As  archived,  they  are  stored  in  a  peculiar  kind  of  IBM  binary, 
no  small  task  for  most  computer  systems  to  read.  The  mere  format  of  the 
data  itself  may  take  weeks  to  learn  and  truly  understand.  One  must  know  the 
data  is  there,  understand  its  content,  discover  how  to  procure  the  data, 
obtain  the  data  itself,  understand  its  format,  learn  to  read  the  tapes, 
develop  software  to  access  data  from  the  tapes,  and  interface  that  data  with 
whatever  software  one  has  available  for  statistical  analysis  before  one  can 
successfully  utilize  this  data  base  at  all. 

One's  task  is  not  yet  done.  There  are  a  variety  of  sources  of 
satellite  pictures,  each  with  its  own  organization,  channels  for  access, 
storage  locations,  etc.  The  data  are  archived  at  various  universities  and 
at  various  government  sites  often  controlled  by  NASA.  To  begin  with,  one 
must  be  aware  that  each  of  these  data  bases  simply  exists;  second,  one  must 
know  of  their  content;  third,  one  must  be  aware  of  how  to  obtain  the  data; 
fourth,  one  must  be  able  to  either  obtain  the  data  in  digitized  form  or  be 
able  to  get  this  data  digitized.  There  are,  of  course,  organizations  such 
as  METSAT,  Inc.,  which  accomplish  digitization,  but  it  is  merely  another 
step  in  the  process  of  trying  to  utilize  the  data.  Only  after  both  of  these 
data  bases  are  identified,  obtained,  read,  integrated  and  analyzed  can  one 
begin  the  analyses  which  integrates  them  to  begin  devising  a  dynamic  cloud 
free  1 ine-of-sight  model.  This  is  asking  a  lot  for  a  few  man-months. 

Barring  the  discovery  of  the  so  called  "sugar  daddy",  this  working 
group  is  unlikely  to  be  able  to  fund  major  data  gathering  or  analyses  in  and 
of  itself.  However,  by  doing  a  certain  amount  of  front  end  work  it  can  open 
the  door  for  many  analyses  to  be  performed  by  analysts  on  different  projects 
for  different  users  across  the  country,  by  making  it  convenient  and  most 
important,  not  terribly  resource  intensive  for  these  analysts  to  obtain  the 
data  and  models  they  need.  At  a  minimum,  I  would  envision  this  working 
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group  endeavoring  to  publicize  as  widely  as  possible  the  data  and  models 
available.  It  could  undertake  to  publish,  as  a  standard  package,  manuals 
for  the  different  data  bases,  including  their  content,  format,  how  they  can 
be  obtained  and  so  forth.  Similar  information  would  be  published  for  the 
models.  An  even  better  solution  would  be  for  this  organization  to,  either 
directly  or  working  with  ETAC  and/or  the  Global  Weather  Center,  establish  a 
central  clearing  house.  In  this  case,  actual  copies  of  the  models  and  of 
the  data  on  some  standard  media  (for  example,  9-track  tape)  could  all  be 
obtained  through  the  clearing  house  in  a  standard  procedure. 


Consider  the  advantages  of  such  a  scheme.  By  one  inquiry  to  one 
well-known  location,  an  analyst  could  discover  the  available  data  and  models 
for  weather-related  problems.  If  this  did  not  include  all  the  resources 
he/she  needed  to  address  the  weather  aspects  of  the  problem,  at  least  all  of 
the  available  data  and  models  would  be  available  with  a  minimum  of  effort. 

It  seems  clear  to  me,  at  least,  that  much  more  could  be  done,  given  the 
data  and  models  available,  than  is  now  being  done.  I  attribute  this 
difference  to  the  large  amount  of  leg  work  required  to  simply  run  down 
what's  there  already.  By  minimizing  this  leg  work,  this  committee  could 
encourage  a  lot  more  analysis,  a  lot  more  model  building,  and  a  lot  more 
validation,  without  having  to  identify  major  funding.  The  funding  will 
already  be  there,  in  the  individual  projects  and  programs  which  have  at 
least  a  few  months  to  spare  for  weather  related  analysis.  I  strongly 
encourage  the  members  of  this  working  group  to  consider  some  way  in  which 
such  a  clearing  house  might  be  established. 

Joel  S.  Davis 
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John  Hornstein 

Naval  Research  Laboratory 

Washington,  DC 
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9  July  198^4 


Dr .  Ernest  Bauer 
Institute  for  Defense  Analysis 
1801  N.  Beauregard  Street 
Alexandria,  VA  2231 1 


Dear  Ernie, 

Here  are  some  notes  and  comments  on  the  discussion 
session  at  the  Second  Tri-Service  Cloud  Modeling  Workshop. 

The  CLOUDS  Program 

The  CLOUDS  utility  package  should  be  designed  so  that 
a  particular  type  of  output  file  (an  analog  would  be  a 
LOWTRAN  transmittance  file)  can  be  produced  by  any  of 
several  modules  providing  the  same  types  of  information  but 
with  different  accuracies  and  degrees  of  detail.  For  this, 
the  format  of  a  given  type  of  output  file  must  be 
essentially  independent  of  the  program  that  produced  it. 
This  can  be  accomplished  by  having  the  file  contain  flags  to 
tell  the  next  program  how  to  read  it.  The  file  itself  would 
consist  of  a  sequence  of  "tables",  such  as  a  transmittance 
table,  various  profile  tables,  a  table  giving  the  irradiance 
for  a  scan  across  a  cloudy  scene,  etc.  The  format  of  each 
type  of  table  would  be  independent  of  the  program  producing 
it,  except  for  the  number  of  lines  and  number  of  digits 
displayed,  which  would  be  specified  by  an  integer  prolog  in 
the  table.  The  resulting  flexibility  would  be  valuable  in 
two  ways : 

(i)  -  A  user  ccould  develop  an  application  using 
qui ckl y- r unn i ng  low-accuracy  versions  of  each  module,  and 
then  switch  to  one  or  more  higher  accuracy  versions  for  a 
final  run,  or  as  an  inexpensive  test  of  how  much  accuracy  is 
needed  for  the  application  at  hand. 

(ii)  -  We  do  not  yet  have  a  clear  idea  of  what  types 
of  output  the  DoD  community  wants.  Moreover,  DoD  needs  will 
continue  to  evolve,  partly  in  response  to  experience  with 
whatever  cloud  modeling  capabilities  are  available.  The 
table-based  file  structure  does  not  force  us  to  wait  until 
all  desired  capabilities  are  fully  specified,  nor  does  it 
force  us  to  freeze  these  capabilities  prematurely.  We  can 
start  immediately  to  design  modules  to  produce  some  initial 
products  -  which  will  help  in  obtaining  funding  -  and  use 
pre-  and  pos t- pr ocess i ng  programs  to  adapt  this  early- 
designed  output  as  future  needs  define  themselves. 


What  about  the  non-DoD  user  community? 


Dr,  Jacobowitz  expressed  NOAA's  interest  in  being 
kept  informed  of  DoD  needs  for  atmospheric  data.  I  know 
several  NASA  scientists  and  contractors  who  would  have  been 
very  interested  in  attending  the  Tri-Service  Workshop.  They 
have  alot  to  contribute:  one  of  them  (Warren  Wiscombe)  is 
world-renowned  for  his  work  on  atmospheric  radiative 
transfer,  which  included  identifying  a  few  parameters  of  a 
particle  size  distribution  that  dominate  its  scattering 
effects;  another  has  unique  data.  There  must  be  others  at 
NCAR  and  in  universities  who  would  like  to  share  information 
with  us . 


We  could  gain  alot  from  this.  Besides  data,  there 
would  be  items  of  knowledge  that  are  difficult  to  acquire 
under  the  constraints  of  working  in  a  m i ss i on- or i ent ea 
environment.  These  people  are  also  a  market  for  our  programs 
(witness  LOWTRAN)  and  data,  and  in  the  course  of  using  them 
they  will  be  testing  and  evaluating  them.  They  may  have  the 
opportunity  to  work  on  problems  we  need  to  have  solved  but 
cannot  get  funded  within  DoD,  and  can  provide  a  constituency 
in  other  agencies  for  work  on  these  problems. 

This  raises  two  obvious  issues.  How  do  we  tell  the 
wider  scientific  community  what  our  problems  are,  without 
divulging  sensitive  information?  How  do  we  offer  a  quid  pro 
quo  of  programs  and  data  in  exchange  for  their  efforts  and 
data,  while  reserving  sensitive  information?  Both  problems 
seem  solvable,  (i)  We  need  a  forum  for  explaining  our  needs. 
Since  a  meeting  in  which  some  sessions  are  closed  and  others 
are  open  is  both  cumbersome  and  insulting,  occassional 
meetings  should  be  for  unclassified  material  only.  To 
preserve  DoD  support  and  avoid  loss  of  focus,  these  meetings 
should  be  even  more  narrowly  oriented  toward  specific  DoD 
needs  than  closed  Tri-Service  workshops,  which  can  and 
should  be  more  wide  ranging,  (ii)  As  LOWTRAN  illustrates,  we 
have  models  and  data  that  can  be  released  to  the  wider 
community  without  revealing  sensitive  aspects  of  the  purpose 
or  means  of  data  collection  or  model  development.  For  use  as 
a  quid  pro  quo,  we  may  have  to  increase  the  volume  of  this 
material,  scrutinizing  the  data  and  models  we  have  to 
identify  additional  material  that  can  safely  be  released 
(with  proper  approval,  of  course).  To  avoid  embarrassment 
and  to  ensure  that  we  have  something  on  hand  to  show 
potential  col  1 eagues  ,  i t  would  be  best  to  identify  such 
material  before  we  actively  solicit  the  participation  of  the 
wider  community.  Should  this  be  dohe  i hd i vi dual  1 y ,  or  should 
there  be  a  DoD-wide  catalog  of  data  and  other  materials  that 
can  be  shared?  The  latter  option  is  related  to  the  clearing¬ 
house  issue,  discussed  hext. 
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The  discussion  panel  broached  the  possibility  of  an 
evaluation  board  to  approve  modules  and  data  before  they  are 
accepted  as  full-fledged  ingredients  of  CLOUDS  or  related 
programs.  I  feel  that  this  should  not  mean  that  "public" 
access  to  a  candidate  module  is  restricted  until  the  module 
has  been  certified.  Much  information  on  validity,  utility 
and  ease  of  use  can  be  generated  by  early  wide  dissemination 
to  users,  accompanied  by  suitable  caveats. 


Who  is  the  user  community  within  DoD? 

Stan  Grigsby  noted  that  one  service  we  can  perform  is 
to  draw  up  a  list  of  applications,  and  for  each  application 
indicate  how  clouds  affect  it,  which  of  these  effects  depend 
on  s t i 1 1 - unr es ol ve d  cloud  issues,  and  the  plausible  impact 
their  resolution  would  have  for  that  application.  This  list 
would  have  several  uses: 

-  they  indicate  some  of  the  capabilities  that  the 
CLOUDS  package  should  provide,  helping  to  define  the  system; 

-  they  show  the  relevance  to  weapons  effectiveness. 
Some  needs  may  have  such  a  narrow  range  of  application  that 
they  will  identify  a  particular  sponsor,  who  must  either 
fund  those  parts  of  the  work  or  expect  to  see  his  needs 
unmet.  (The  application  may  be  important  despite  being 
narrow.)  This  would  be  especially  favorable  for  getting  the 
program  off  the  ground,  since  it  gets  around  the  reluctance 
of  sponsors  to  fund  capabilities  that  other  people  will  also 
use.  On  the  other  hand,  some  needs  will  be  relevant  to  so 
many  applications  that  their  cumulative  importance  will 
elicit  funding,  or  provoke  NASA  or  NCAR  (NSF)  or  NOAA 
interest.  As  several  attendees  stressed,  the  key  to  DoD 
funding  (or  NASA/NCAR/NOAA  activity  on  a  problem  relevant  to 
our  needs)  is  for  the  "cloud  community"  to  push  a  fixed  set 
of  problems  and  capabilities  that  respond  to  a  few 
identifiable  and  documented  needs. 

What  person  or  group  will  start  drawing  up  this  list? 


Another  aspect  of  "What  is  the  DoD  user  community?" 
is  "short  wave  chauvinism".  The  Workshop  stressed  IR  and 
optical  aspects,  bit  both  designers  and  users  of  passive  and 
active  microwave  systems  suffer  from  cloud  (mostly 
pr eci pi t at i on- i nduce d )  clutter  and  attenuation.  Some  of 
these  people  have  considerable  interest  and  expertise  in 
CFLOS  and  cloud  microphysics  and  dynamics.  They  should 
probably  be  invited  to  Join  us.  Besides  bringing  their 
knowledge  and  data,  they  are  a  "market"  for  our  models  and 
data,  increasing  our  constituency. 
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Several  attendees  wanted  a  central  clearing-house  for 
cloud  data,  programs  and  lore.  The  panel  mentioned  that  it 
had  neither  funds  nor  a  mandate  for  a  clearing-house  yet.  At 
present  it  can  serve  only  as  an  unofficial  contact  point. 
One  attendee  stressed  that  a  great  deal  of  time  and  effort 
could  be  saved  if  one  person  or  group  undertook  to  at  least 
know  the  location  of  all  data,  programs  and  experts.  Someane 

One  attendee  raised  the  related  issue  of  efficient 
access  for  newcomers  to  the  lore;  typical  microphysical  and 
optical  parameters  for  various  types  of  clouds;  relative 
importance  of  scattering  and  absorption  in  various  spectral 
regions;  etc.  Since  this  material  doesn't  change  rapidly,  it 
doesn't  seem  necessary  to  have  a  person  or  group  "on  call" 
with  this  material.  The  need  could  be  met  by  a  good  chapter 
devoted  solely  to  clouds  and  fog,  in  something  like  the 
Handbook  of  Geophysics  or  the  IR  Handbook;  preferably  the 
former,  since  the  microwave  community  also  needs  this 
material.  The  current  treatments  in  these  two  books  are  too 
scanty  and  out  of  date,  and  neglect  some  important  topics. 


A  problem  with  contract  studies. 


A  final  issue  of  access  concerns  DoD  access  to 
intermediate  results  of  contract  studies.  Ron  Nelson  jaised 
this  issue,  but  it  has  also  been  a  problem  for  our  group. 
Quite  often  intermediate  results  or  reasonings  would  be 
useful  in  themselves  or  for  evaluating  the  suitability  of  a 
contract  deliverable  for  a  possible  application.  Attempts  to 
obtain  these  final  but  intermediate  results  frequently 
encounter  a  stone  wall.  What  is  needed  is  not  a  formal 
report,  which  is  troublesome  and  expensive,  just  a 
recognition  that  this  material  is  not  proprietary,  and  a 
willingness  to  discuss  it.  Perhaps  the  panel  of  the  Tri- 
Service  workshop  knows  a  few  people  at  appropriate  levels 
who  could  insist  that  contract  work  sponsored  by  their 
commands  must  make  their  final  intermediate  results  and 
methods  of  reasoning  available  -  in  an  informal  way  and  only 
on  demand  -  to  qualified  requestors. 


John  Hornstein 
Code  6521 

Naval  Research  Laboratory 
Washington,  DC  20375 


:  Donald  Grantham  (AFGL) 
Edward  Stone  (NRL) 


(202)  767-3069 


DEPARTMENT  OF  THE  ARMY 

HARRY  DIAMOND  LABORATORIES 
2800  POWDER  MILL  RO..  ADELPHI.  MD  20783 


REPLY  TO 
ATTENTION  OF 

DELHD-RT-CB  28  June  1984 

SUBJECT:  Suggestions  to  Tri-Service  Cloud  Modeling  Committee 


Commander 

Air  Force  Geophysics  Laboratory 
Atmospheric  Sciences  Division 
Tropospheric  Structure  Branch 
ATTN:  LYT/Dr.  D  Grantham 
Hanscom  Air  Force  Base,  MA  01731 


1.  First,  I  would  like  to  thank  you  and  your  collaborators  for 
putting  on  such  a  fine  workshop.  I  found  it  a  very  worthwhile 
experience. 

2.  In  response  to  your  request  during  the  wrap-up  session  for 
written  comments,  I  would  like  to  endorse  the  suggestion  by 
various  people  that  a  recommended  program  be  put  together  now. 
This  program  should  consist  of  high  priority  unfunded  tasks  that 
the  committee  believes  are  a  necessary  part  of  an  overall  effort 
that  should  be  pursued  as  soon  as  funds  become  available. 

3.  One  such  task,  mentioned  by  Dr.  Ernest  Bauer  of  IDA  during 
the  wrap-up,  concerns  experimental  cloud  edge  data.  The  need  for 
more  such  data  became  quite  evident  during  the  workshop.  As  you 
know,  the  Harry  Diamond  Laboratories  have  the  measurement 
capability  for  obtaining  such  data,  and  have  already  obtained  a 
significant  data  bank  containing  cloud  edge  data.  Examples  of 
these  data  were  presented  at  the  First  Cloud  Modeling  Workshop 
and  published  in  the  collected  viewgraphs,  and  more  recently 
presented  at  the  IRIS  Targets,  Backgrounds,  and  Discrimination 
Subgroup  Meeting  and  published  in  its  proceedings  (just  out). 
Reduction  and  analysis  of  more  of  these  data,  which  are  to 
resolutions  as  fine  as  =1.5  m,  should  significantly  enhance  our 
understanding  of  cloud  edges. 

FOR  THE  DIRECTOR: 


Z.  G.  SZTANKAY  ^ 

Chief,  Near  Millimeter 
Wave  Branch 


PART  VI 


ON  THE  USE  OF  AIRBORNE  LIDARS  TO  COLLECT 
METEOROLOGICAL  DATA  (WINDS,  DENSITIES, 
AEROSOLS)  DURING  MISSILE  REENTRY 
DURING  BROAD  OCEAN  AREAS 


(Report  of  a  meeting  on  above  subject,  held  at  the 
Institute  for  Defense  Analyses  on  24  January  1984) 


SCIENCE  AND  TECHNOLOGY  DIVISION 

1801  N.  Beauregard  Street,  Alexandria,  Virginia  22311  •  Telephone  (703)  845-2000 

iNsnTUTsro* 

OIFINSI  ANALYSES 

25  January  1984 

MEMORANDUM  FOR  DISTRIBUTION 

SUBJECT:  On  the  Use  of  Airborne  Lidars  to  Collect  Meteorological 
Data  (Winds , Dens i ties , Aerosols )  During  Missile  Reentry 
over  Broad  Ocean  Areas (BOA).  Pre  liminary-f or  Review. 

FROM:  Ernest  B^ITer  and  LCDR  Sta^46y  Gr igsliy  (NAVSEA) 

On  24  January  1984  an  unclassified  meeting  on  the  referenced 
subject  was  held  at  IDA.  The  list  of  attendees  is  included  as 
Attachment  1. 

The  reason  for  the  meeting  was  a  response  to  needs  by  the  Navy 
Fleet  Ballistic  Missile  Office  (SSPO)  for  weather  measurements 
during  RV  reentry  over  BOA's  ,  principally  in  the  S. Atlantic  Ocean. 
The  problem  was  defined  by  R. Sokol  ( SP-2 7 , Reen try  engineer  for  the 
D-5  vehicle)  and  by  Jay  Berkowitz  {SP-25,  Test  engineer)  who  stated 
that  current  primary  requirements  are  for  winds  and  densities  at 
altitudes  to  50/100  Kft,  with  less  emphasis  on 

a.  large  ice  crystals  (  ^  300  ^m) 

b.  higher  altitudes  (to  250-350  Kft) . 

Earlier  concern  with  cirrus  has  been  reduced  by  aerodynamic  redesign 
leading  to  higher  altitudes  of  transition  to  turbulence  on  the 
vehicle . 

Lee  Dickinson  (LMSC)  reported  on  the  ADMSS  (-Aircraft  Deployed 
Meteorological  Sounding  System)  study,  which  was  conducted  some  two 
years  ago. by  Space  Data  Corp.  It  included  a  review  (by  Kentron  Corp) 
of  indirect  sensing  which  concluded  that  at  the  time  of  the  study 
lidar  and  other  remot-e  systems  were  not  suitable  for  gathering  the 
necessary  data,  although  the  technology  was  progressing  rapidly. 

The  ADMSS  study  recommended  using  dropsondes  and  sounding  rockets. 
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all  deployed  from  a  P-3  aircraft.  The  recommendation  has  not 
been  acted  upon,  principally  because  of  the  high  risk  and  cost 
of  ensuring  safety  for  the  aircraft  during  the  rocket  launches. 

It  was  stated  that  meteorological  data  are  actually  used  in 
computer  programs  with  1000  ft.  vertical  resolution. 

Freeman  Hall  (NOAA/WPL)  reported  on  their  ground-based  Doppler 
lidar  system  for  measuring  winds;  Geoff  Kent  (IFAORS-  substituting 
for  Pat  McCormick  of  NASA/LaRC)  reported  on  airborne  aerosol 
measurements,  and  Tom  Wilkerson  (U.  of  Maryland)  reported  on  the 
capabilities  of  DIAL  and  other  lidars  to  measure  temperature  and 
pressure  (from  which  density  is  derived). 

Capt.  George  Fisher  (USAF  ;  BMO/WE)  indicated  BMO  concerns  and 
needs . 

The  conclusions  of  the  group  are  listed  in  Attachment  2.  They 
represent  a  consensus  of  these  experts  in  the  lidar  area. 

The  following  material  is  available  from  IDA  on  request: 

Attachment  1.  Attendees 

2.  Conclusions 

3.  Requirements  (R. Sokol) 

4.  AOMSS  briefing  (L. Dickinson) 

5.  Section  4  of  ADMSS  report  on  Indirect  Sensing 

6.  Remote  wind  sensing  by  IR  Doppler  lidar  (F.Hall) 

7.  Lidar  measurements  of  aerosols  and  clouds 
(G.Kent  /  NASA) 

Atmospheric  lidar  (T. Wilkerson) 


8. 


ATTACHMENT  1 


Attendees . 


Ernie  Bauer 


703-845-2290 


Jay  B.Berkowitz 


SSPO 


202-697-0600 


Roger  Carson 


Lee  G.  Dickinson 


LCDR  Stan  Grigsby 


NSWC/Dahlgren 


LMSC 


Capt.  George  F. Fisher  BMO/WE 


703-663-8740  or  8138 


408-742-4956 


714-382-6891 


NAVSEA/PMS-405  202-692-5626 


Freeman  Hall 


NOAA/WPL 


303-497-6312 


Jerry  L. Harper 


Kaman  Sci.Corp.  303-599-1931 


Geoff  S.Kent 


IFAORS 

(or  NASA/LaRC 


804-865-0811 

804-865-2065) 


Susan  Masters 


NSWC/Dahlgren 


703-663-8346 


Richard  Sokol 


SSPO 


202-697-1352 


Tom  Wilkerson 


U,  of  MD. 


301-454-5401 
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CONCLUSIONS. 


ATTACHMENT  2.  Conclusions. 
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e  upper  altitude  limit  to  wind  measurements  is  set 
the  low  aerosol  density  above  20-25  km,  especially  in 
volcanically  quiet  times. 


